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Abstract 


The Loop-Tree Duality (LTD) is a novel perturbative method in QFT that establishes 
a relation between loop-level and tree-level scattering amplitudes. This is achieved by 
directly applying the Residue Theorem to the loop-energy-integration. The result is a 
sum over all possible single cuts of the Feynman diagram in consideration integrated 
over a modified phase-space. These single-cut integrals, called Dual contributions, are 
in fact tree-level objects and thus give rise to the opportunity of bringing loop- and 
tree-contributions together, treating them simultaneously in a common Monte Carlo 
event generator. Initially introduced for one-loop scalar integrals, the applicability of 
the LTD has been expanded ever since. In this thesis, we show how to deal with Feyn¬ 
man graphs beyond simple poles by taking advantage of Integration By Parts (IBP) 
relations. Furthermore, we investigate the cancellation of singularities among Dual con¬ 
tributions as well as between real and virtual corrections. For the first time, a numerical 
implementation of the LTD was done in the form of a computer program that calculates 
one-loop scattering diagrams. We present details on the contour deformation employed 
alongside the results for scalar integrals up to the pentagon- and tensor integrals up to 
the hexagon-level. 


Resumen 

La Dualidad Loop-Arbol (LTD) representa un nuevo metodo perturbative en Teoria 
Cuantica de Campos que establece una relacion entre amplitudes de dispersion virtuales 
y de arbol. Se logra hacer esto por aplicacion directa del Teorema de los Residues a la 
integration de la componente de energia. El resultado es la suma de todos los cortes 
simples posibles del diagrama de Feynman considerado integrada sobre un espacio fasico 
modificado. Estas integrales de corte simple, denominadas Contribuciones Duales, de 
hecho son objetos de tipo arbol y por lo tanto dan lugar a la oportunidad de combinar 
las contribuciones virtuales y de arbol con el motive de tratarlas simultaneamente en 
un generador de eventos de Monte Carlo. A pesar de ser introducido inicialmente para 
integrales escalares de un loop, la practicabilidad de la LTD fue extendida tremenda- 
mente. En esta tesis demonstramos como aplicar la LTD a diagramas con polos de 
multiplicidad elevada utilizando relaciones de Integracion Por Partes (IBP). Ademas, 
examinamos la cancelacion de singularidades entre Contribuciones Duales tanto como 
entre correcciones reales y virtuales. Por primera vez una implementacion numerica 
de la LTD fue realizada en forma de un programa de ordenador que calcula diagramas 
de dispersion. Presentamos detalles sobre la deformacion de contorno empleada y los 
resultados de integrales escalares hasta el nivel de pentagono y de integrales tensoriales 
hasta el nivel de hexagono. 
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Chapter 1 


Introduction 


The Large Hadron Collider (LHC), the largest and most complex machine ever built 
by mankind, gives particle physicists a powerful tool at hand to verify existing models 
and probe the fundamental laws at very high energies. The aim of the LHC is: First, 
to investigate whether the Standard Model (SM) is still valid at the collider’s energies. 
Second, to shed light on the electroweak symmetry breaking mechanism and to find or 
exclude a particle that fits the SM description of the Higgs boson [7-9]. Third, to search 
for Beyond Standard Model (BSM) physics like supersymmetry or extra-dimensions or 
particles that could be Dark Matter candidates. In the first run of the LHC, a Higgs- 
like particle has been found with a mass of 125 GeV. Apart from that, several other 
discoveries have been made including the first creation of a quark-gluon plasma or the 
rare B^-decay. In the second run, in which the center-of-mass energies will be increased 
even further, signals of BSM physics are hoped to be detected and the properties of the 
Higgs further explored. 

Despite the existence of observations which the SM cannot accommodate at the mo¬ 
ment, e.g. the existence Dark Matter [10] or neutrino oscillations [11], it is still the only 
established theory of particle physics to date that describes experimental data, a great 
theoretical achievement of its own. Leaving out gravity, for which there is no quantum 
theory available yet, it successfully describes (almost) all relevant particle physical ob¬ 
servables. 

The SM is a relativistic quantum field theory with an SU{3) x SU{2) x U{1) gauge 
symmetry. It describes three of the four fundamental forces of nature at microscopic 
distances, namely electromagnetism, strong and weak force. ^ 

The SU (3) symmetry accounts for the strong force and the corresponding quantum field 
theory is called QCD (Quantum Chromodynamics, [14]). It describes the interactions 
^For a full review of the Standard Model, see for example [12, 13]. 
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between quarks and gluons. Quarks are, besides leptons, the fundamental matter con¬ 
stituents and come in six different flavors. Gluons are massless spin 1 particles (bosons) 
carrying an SU{3) colour charge and act as the mediators of the strong force. Because 
the gauge group is non-Abelian [15], QCD features two remarkable properties: Confine¬ 
ment [16] and Asymptotic Freedom [17, 18]. 

Confinement accounts for the fact that physical objects are always colour-neutral at low 
energies, in particular no individual free quarks or gluons are observed. Nonetheless, an 
analytical proof of this property is still missing. 

Asymptotic Freedom is a property of the theory. It means that due to the running of 
the strong coupling the strong force becomes small at high energies (equivalentely, 
small distances), a fact that permits the employment of perturbative techniques. 

The weak force describes processes such as the decay of nuclei and the interaction of 
neutrinos with matter. In the modern context, it is better understood within the frame¬ 
work of the electroweak sector. By electroweak sector, we mean the unification of weak 
and electromagnetic forces [19-22] by Weinberg, Salam and Glashow, a major success 
which allowed to understand both forces in the common framework of the Electroweak 
Symmetry Breaking mechanism. Very similar to the strong force, its gauge group is 
the non-Abelian SU{2). Contrary to QCD, the W^- and Z^-gauge bosons of the weak 
interaction are not only massive, but with masses of 80 GeV and 91 GeV respectively, 
they are quite heavy. Hence it is a short distance interaction and appears to be ‘weak’. 
Finally, 17(1) is the gauge group of Quantum Electrodynamics (QED) [23] which is the 
quantum field theory describing the electromagnetic force. Its gauge boson, the photon, 
is massless and due to the Abelian nature of t7(l), photons do not interact with each 
other and the force is of unlimited range. 

What we today know as the Standard Model of particle physics is QCD together with 
the electroweak sector. It has undergone countless checks and investigations over many 
different aspects and it has been exceptionally successful in making correct and accurate 
predictions for a wide range of physical observables. 

However, the SM is believed to be a mere low-energy approximation to a still to be con¬ 
structed unified field theory which would describe all four forces at all energy regimes. 
The SM-Lagrangian encodes all previosly known symmetry principles, for example con¬ 
servation laws. Although its mathematical formulation is very simple, we cannot an¬ 
alytically solve the equations of motion. Only in certain parts of phase-space we can 
perform reliable theoretical calculations which are all based on perturbation theory. The 
underlying idea is that when the coupling of the interaction term is small, one can do a 
series expansion in which every term can be represented in a pictorial form. This is done 
by the so-called Eeynman diagrams [24, 25]. Due to its clarity and predictive power, 
the diagrammatical approach is the most popular for theoretical calculations. The first 
term of the expansion usually gives an estimate of the order of magnitude whereas in 
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order to obtain the first proper estimate, one has formally to go to next-to-leading order 
(NLO) precision. More and more processes processes demand next-to-next-to-leading 
order (NNLO) precision to match the precision of the experimental data. 

The discovery of the Higgs-like boson in 2012 has been a huge success [26, 27]. At 
the moment of writing this thesis, the LHC is warming up for its second phase with a 
center-of-mass energy of up to 13 TeV. In this second phase, as previously mentioned, it 
will be of great importance to measure as many properties of the discovered particle as 
possible and to continue the hunt for physics beyond the Standard Model. So far, the 
absence of signals hinting to BSM physics is a little disappointing as it does not dehne 
a clear direction for the theorists to follow. 

The high quality of LHC data raises the need for high-precision theoretical predictions. 
The processes at the LHC are rather challenging to calculate because they typically 
involve many particles and because QCD plays the dominant role at the LHC. Fur¬ 
thermore, higher orders of the perturbation expansion have to be calculated in order to 
match the experimental precision. This has led to considerable progress in the analytical 
and numerical techniques for the calculation of Standard Model cross-sections. Apart 
from the usual diagrammatic approach, there are other methods, some of the most pop¬ 
ular ones being Unitarity Methods [28-30], the OPP-Method [31, 32], Mellin-Barnes 
Representations [33, 34] and Sector Decomposition [35-38]. Thanks to these techniques, 
2—7-4 processes at NLO are the standard nowadays, and even higher multiplicities are 
becoming more accessible [39-42]. They have achieved an incredible feat: The com¬ 
putation of Feynman graphs up to NNLO-level and in some cases even beyond. Still, 
many important issues remain. When calculating cross-sections one needs to consider 
tree- and loop-contributions separately. Thus, a lot of effort has to be put into cancelling 
infrared singularities between real and virtual corrections [43-47] . Additional difficulties 
arise from threshold singularities that lead to numerical instabilities. 

Recently, a new method called the Loop-Tree Duality (LTD) [48] has been developed, 
which is designed to attack these problems. The basic concept at one-loop is to directly 
apply the Cauchy Residue Theorem to the Feynman integrals. The outcome is a sum 
of tree level objects in which each represents all possible single cuts of the considered 
diagram. This form is called the ‘dual integral’, which closely resembles the real correc¬ 
tions. The idea is to then combine the dual integral with the tree-level contributions in 
order to treat them simultaneously in a common Monte Carlo event generator. While 
initially the technique was limited to one-loop graphs, it has been greatly expanded since 
then. In [49] it has been shown how to extend it to diagrams with an arbitrary number 
of loops and in [1] how to deal with graphs which involve propagators that are raised to 
higher powers (higher order poles). 
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1.1 Outline 

The remainder of this thesis is organised as follows: In Chapter 2, we establish the 
fundamentals of this work. In Chapter 3 we introduce the Loop-Tree Duality method 
alongside some illustrative examples. In Chapter 4, we formalize the notation and 
extend the Loop-Tree Duality to double and higher order loop graphs, in Chapter 5 
we show how to deal with poles of higher multiplicities. In Chapter 6 we report on 
the cancellation of singularities among dual contributions as well as between real and 
virtual corrections for massless internal lines. In Chapter 7, we present details on the 
numerical implementation of the Loop-Tree Duality for scalar one-loop integrals. This 
is the first time that the LTD has been applied to explicitly calculate Feynman diagrams 
and constitutes the main result of this thesis. In Chapter 8, we demonstrate that the 
computer program used in Chapter 7 is also able to deal with tensor integrals. We 
conclude the thesis with Chapter 9 in which we give a summary and our future plans. 
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Standard Model Phenomenology 


The Standard Model is a theory of particle physics that has been refined continuously 
over the past, so that it is nowadays able to describe almost all particle reaction processes 
that we observe in the laboratory (particle colliders) as well as in space. Its tremendous 
success is threefold: 

• The SM has the ability to explain a wide variety of experimental results. 

• The SM repeatedly predicted the existence of particles before their experimental 
discovery. This has been the case for the and bosons [50, 51], the gluon 
[52-55] and the charm [56, 57] and top quarks [58, 59] and very recently the 
Higgs boson [26, 27]. For each of these particles, experiments later confirmed the 
predicted properties with good precision. 

• The SM has passed a huge number of precision tests with flying colours, the most 
famous one being the anomalous magnetic dipole moment of the electron [60]: 
The agreement between theory and measurement is up to the 13th digit, which 
is a unique achievement not only in particle physics, but in all science. Other 
examples of precisely predicted quantities are the Neutron Compton Wavelength 
[61] or the mass of the boson [62]. 

All of these accomplishments have strengthend the conhdence in the SM as the proper 
theory for the description of the behaviour of elementary particles excluding gravity. 
We dedicate the next section to an instructive review of the basic principles of the SM. 
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2.1 The Standard Model of Particle Physics 

The Standard Model of particle physics is a quantum field theory which nowadays allows 
us to correctly describe almost all physical processes with high precision on a fundamen¬ 
tal level. It is based on the gauge group 

SU{?,) X SU{2) X U{1) (2.1) 

where SU (3) is the gauge group of the strong interaction, SU (2) the gauge group of the 
weak isospin, and 17(1) the gauge group of the hypercharge which is not identical with 
the gauge group of QED. To avoid confusion, the notations U{1)y and U{1)em with 
subscripts hinting to either the hypercharge (Y) or QED (EM) will be used. 

Despite the interacting constituents being fields, the SM’s perturbative description is 
carried out almost entirely in the particle picture. We imagine the force between two 
matter particles to be mediated via exchange particles. All interactions have in common 
that the matter fields are spin ^ fermions and the exchange particles spin 1 bosons, with 
the exception of the Higgs boson which has spin 0. 

2.1.1 QCD 

The QCD, Quantum ChromoDynamics, is a non-Abelian gauge theory with the action 

S = y - m)V’- |trF^,yF^^| , (2.2) 

and the symmetry group SU{3). The fermion fields V’ transform under its fundamental 
representation. They are, apart from the gluons, the dynamical degrees of freedom of 
the theory. There are six flavors of quarks, each with a different mass. D is the so-called 
covariant derivative, defined as 


D^ = d^,- igsAf, , (2.3) 

with = A‘IjT°' being the gluon field, which lives in the adjoint representation of 
the gauge group. The group SU{N) has N'^ — 1 generators, hence SU{3) has eight 
of them. They span the Lie algebra su(3) and are normalised by tr(T“r^) = and 
[T“,r^] = is the field strength tensor. It is an algebra-valued two-form 

defined by 

F^Jiv = = d^Ay - dyA^ + igs[A^,Ay]. (2.4) 
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The first two terms are a four-dimensional rotation known from QED; the third term is 
typical for non-Abelian symmetry groups and the origin of the gluon self-interaction. 
Considering QCD in the absence of quarks, the only thing left is the gluon self-interaction. 
This branch is called pure Yang-Mills theory. As we mentined in the Introduction, QCD 
has two special properties, Confinement and Asymptotic Freedom. They mainly have 
to do with the running coupling as{Q) of QCD. By running coupling we mean that the 
coupling is dependent on the energy at which the process happens [12]: 

27r 

= (ll-2n^/3)ln(Q/AQCD) ' 

Confinement takes place at low energies at which the coupling is strong. Physical states 
are always color singlets: As a consequence of the non-Abelian gauge symmetry, the 
energy cost grows proportionally to the distance, if one tries to separate the particles 
of a color singlet. At some point it is energetically favorable for a new quark-antiquark 
pair to appear. Hence, neither free quarks nor free gluons are observed in nature. 
Asymptotic freedom describes the phenomenon that the coupling of the theory becomes 
weaker as the energy increases. At very high energies, quarks and gluons are almost 
able to move as if they were free particles. In this energy regime, perturbation theory 
is applicable. However, one has to quantise the action first, a good method to do so 
is the one by Faddeev and Popov [63]. An unwanted side effect of the method is the 
introduction of ghost fields into the theory. They are called ghosts because of their 
incorrect spin-statistic relation and therefore are unphysical. To get rid of them, one 
can choose the axial gauge, in which ghosts and gluons decouple and the gauge of the 
action is fixed. From the gauge fixed action one can extract the propagators and vertices 
which allow to calculate the contributions to the perturbation series via the Feynman 
diagram approach. These diagarams are pictorial representations of the terms of the 
perturbation series. To calculate an amplitude of arbitrary order in ccs, one draws all 
the Feynman diagrams that contribute to the process, contracts them with the outer 
polarisation vectors and spinors and finally sums them up. A downside of this approach 

P 

/■-*- 7 

P ^ 

p, a v, b 

Figure 2.1: QCD propagators in Feynman gauge, without ghosts. 

is that it consumes a lot of computational power. For example, the number of Feynman 
diagrams that contribute to the n-gluon amplitude grows faster than n![64]. Techniques 
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P, a 



y, b 


gfabc[g,u^p _ 

+g‘'^{q — pY 

+gPYk - qf] 



—ig‘i[jabejcde^gtipgua _ gpagUp'^ 
+ fbCe^gPUgPa _ gPPgUa^ 

+ fbde ^gPUgpa _ gpagUpy^ 


Figure 2.2: QCD vertices in Feynman gauge, without ghosts. All momenta outgoing. 

like color-ordered amplitudes [65] and (tree-level) recurrence relations [66-68] have been 
developed to deal with that problem. 

2.1.2 Electroweak Interaction 

Until 1967, QED and weak interaction were two separate theories. In that year, Salam, 
Weinberg and Glashow succeeded in understanding both interactions as special cases of 
one unified theory called Electroweak theory [19-22]. The reason why both interactions 
appear to be so different is the spontaneously broken symmetry of the Electroweak 
theory. 

The corresponding action resembles the one of the QCD: 

5 = j d^x |v;(0 - mY - • 

The matter fields of the electroweak interaction are the six leptons, bnt quarks do 
carry a weak charge as well. They can be organised into three generations. Electron, 
muon and tau participate in both the weak and the electromagnetic interactions. The 
corresponding neutrinos interact exclusively weakly. Within the SM, they are supposed 
to be massless, though recent experiments show that they actually have a small nonzero 
mass. This allows to solve the solar nentrino problem using the concept of nentrino 
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oscillations [11]. The operator D is defined as: 



( 2 . 6 ) 


The underlying symmetry group of the electroweak interaction is the SU{2) x U{\)y- 
In Eq. (2.6), the (t“ are the generators of the SU{2) and the hypercharge Y is the 
generator of the U{1)y- Because of the direct product of the two groups, they carry 
different coupling constants g and g'. 

The fields W and B are the gauge fields of the theory. IT is a field similar gluons 
of QCD with respect to their non-Abelian character. The corresponding field strength 
tensor is obtained by replacing A with W in equation (2.4). The field B transforms 
under the Abelian symmetry group U{1)y- Thus its field strength tensor does not have 
a commutator. 

According to the way B and W appear in the Lagrangian, they are massless. However, 
experiments show that the W and Z bosons indeed have a mass. To get rid of this flaw 
of the theory, one introduces an additional field 4>, the so-called Higgs-held, that couples 
to the gauge bosons by the term The a priori unphysical and massless 

helds B and W now obtain their masses via the Higgs-mechanism. By transforming the 
unphysical helds B and W in physical ones, the SU{2) x C/(l)y-symmetry is broken 
down to a 17(l)£;A^-symmetry. Simultaneously the Lagrangian picks up the correct mass 
terms for the IT and Z bosons as well as the photon. The other terms of the Langrangian 
represent the potential of the Higgs held. This is illustrated in greater detail in the next 
subsection. 

Despite the similarities between QCD and electroweak interaction, the latter one does 
not display the phenomenon “Conhnement”. The reason for this is that the exchange 
bosons of the electroweak interaction are very heavy, namely 80 GeV for the IT^ and 
91 GeV for the Z^, which does not allow for bound states. “Asymptotic Freedom”, 
however, is observed for the weak interaction just as in QCD. 

2.1.3 Higgs Boson and Electroweak Symmetry Breaking 

In the Standard Model, the gauge group SUi2) x U{1)y is spontaneously broken. The 
left-over symmetry is U{1)em- One assumes the existence of an additional complex 
scalar held which transforms under SU(2) and has hypercharge T = 1. In the space 
of weak isospin, it can be represented by a two-component vector with complex entries 


[69]. 



(2.7) 
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(/>■*“ (x) is a complex field, i.e. it is composed of two components. The three components 
(/>■'“(x) and x{x) are absorbed into the longitudinal modes of and Z^. H{x) is the 
actual physical Higgs-held. The Lagrangian of the Higgs-sector then reads 

^Cffiggs = Yukawa , (2.8) 

'-v-" 

V{4>) 

For the theory to be bound from below, A has to be greater than zero, A > 0, and 
furthermore vn? < 0 for the symmetry to be spontaneously broken. The covariant 
derivative is given by Eq. (2.6). The Higgs-dublet has hypercharge T = 1; from the 
minimization of the potential V{(j)) follows 


\4>r, 



(2.9) 


In the next step, the physical contributions to {Dcj) will be calculated. 


(l)^</.)t(Z)^</))|phys 


= + 


5-'^ -99'\ ( 

-99' 9^ J 


( 2 . 10 ) 


H is already a physical held, but B and W are not. They are transformed to physical 
helds by virtue of the transformations 




( 2 . 11 ) 


and 




cos^w sin0ty\ / 
-sin% cosew) Wu, 


Here, sin Ow and cos 9w are given by 

9 


cos 9w = 


\/pTff 


and sin0pi/ = 


This way, one can rewrite the part of the Lagrangian of equation (2.10) as 

(D^</-)t(Z^/^,/,)|phys = + 2 (f)' (W;*VF+^ + W-*W-f^) 


( 2 . 12 ) 


(2.13) 


(2.14) 
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Now, the masses of the W- and Z-bosons can be read off easily: 


= ^9 and mz = + 


,/2 


2 -^ "2 
The couplings g and g' are related to the elementary charge e by 


(2.15) 


e = 


99 




,/2 


or g = - 


sin 6w 


and g' = 


cos 9w 


respectively. 


(2.16) 


Within this framework, it is possible to have massive gauge bosons which otherwise are 
forbidden by unbroken gauge symmetries. 


2.2 Particle Phenomenology 

Phenomenology is the branch of particle physics that deals with the calculation of phys¬ 
ical observables that can then be compared to experimental measurements. In order to 
test a theory, solutions have to be compared to experimental data. For the SM, there are 
two established ways to calculate them: Lattice calculations and perturbation theory. 
Since lattice techniques lie beyond the scope of this thesis, we focus on the latter. As 
mentioned in Section 2.1.1, perturbation theory (in the coupling) is employed at high 
energies when the coupling of the theory becomes small. Usually, Feynman diagrams 
are used to calculate the individual contributions to the perturbation series. We distin¬ 
guish betweeen two main types of diagrams: Tree and loop diagrams. The the leading 
order term is a tree-level contribution. Beyond leading order, a Feynman diagram gets 
dressed with more lines. These can connect two existing lines, thus leading to a loop 
or they connect to only one existing line resulting in an additional external leg. We 
use the terms next-to-leading order (NLO, one-loop) or next-to-next-to-leading order 
(NNLO, two-loops) to indicate the accuracy up to which a certain calculation has been 
performed. Today, the quality of data collected at colliders like the LHC has reached 
a level that makes the inclusion of terms beyond the leading order of perturbation the¬ 
ory necessary in order to provide theoretical estimates of similar precision. Apart from 
the cancellation of singularities between tree and loop corrections and phase-space in¬ 
tegrations with many external legs, loop diagrams are the bottleneck when calculating 
scattering amplitudes. Hence, theorists are looking for ways to efficiently calculate NLO 
and NNLO corrections. The last decade in particular has seen a huge progress in the 
development of calculational techniques to tackle computations at NLO and also some 
progress at NNLO accuracy. In the LHC era, this is important, because of the avail¬ 
ability of large amounts of data to which theoretical predictions can be compared to. 
Still, the task of obtaining theoretical estimates that can be compared with LHC-data 
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is not easy. Tree and loop contributions have to be evaluated independently before 
bringing them together. The Loop-Tree Duality method, which will be introduced in 
Chapter 3 and which is at the heart of this thesis, reexpresses loop-level objects in terms 
of tree-level objects. Thus it allows to directly combine the two, in order to evaluate 
them simultaneously. This evaluation is carried out using numerical integration. In the 
latter chapters of the thesis, we present results of the of the first implementation of the 
Loop-Tree Duality method as a computer progam. This program calculates one-loop 
integrals using a numerical integrator (Cuhre or VEGAS). Since the numerical integra¬ 
tor is employed as a black box, we recapitulate the basics of numerical integration in 
Section 2.3. 


2.3 Numerical Integration Techniques 

2.3.1 Motivation for employing numerical techniques 

In Chapters 7 and 8, we are going to calculate a certain type of three-dimensional 
integral. 

Usually one is interested in knowing the value of an integral up to a certain accuracy 
and wants to obtain the result within an acceptable amount of time. Hence, a numerical 
integrator in the form of a computer program is the proper tool for this task. 

Our long-term goal is to set up a fully automated program which is able to calculate 
(N)NLO cross-sections for arbitrary kinematics. In this thesis, we do the first step by 
addressing the virtual part. 

2.3.2 Classical numerical integration 

For the rest of this chapter, we closely follow the paper [70] by Weinzierl, which is a nice 
introduction to the subject and is recommended to readers with further interest in the 
matter. 

Quadrature rules have been known for a long time. The numerical integrator “Cuhre”, 
part of the Cuba-library [71], uses quadrature rules to estimate the integral, therefore 
their idea shall be illustrated here. One distinguishes between formulae that evalu¬ 
ate the integrand at equally spaced abscissas (Newton-Cotes type) and formulae which 
evaluate the integrand at carefully selected, but non-equally spaced abscissas (Gaussian 
quadrature rules). The simplest example of a Newton-Cotes type rule is the so-called 
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trapezoidal rule: 

a:o+Aa: 

J dxf{x) = ^[/(xo) + /(xo + Ax)] - > (2-17) 

where xq < ^ < xq + Ax. To approximate an integral over a finite interval [xq, Xr ^ with 
the help of this formula, one divides the interval into n sub-intervals of length Ax and 
applies the trapezoidal rule to each sub-interval. With the notation Xj = xq + j ■ Ax, 
one arrives at the compound formula 


Xn 

Jdxf{x) 

XQ 


Xn Xq 

n 


j=0 


1 {Xn - Xq)^ 
12 n2 ^ 


(2.18) 


with wq = Wn = 1/2 and Wj = 1 for 1 < j < n — 1. Further 


/" 


1 


(2.19) 


where is somewhere in the interval [xj_i,Xj]. Since the position of the cannot 
be known without knowing the integral exactly, the last term in Eq. (2.18) is usually 
neglected and introduces an error in the numerical evaluation. This error is proportional 
to 1/n^ and one has to evaluate the function /(x) roughly re-times. 

An improvement is given by Simpson’s rule, which evaluates the function at three points: 

X2 g 

J dxf{x) = ^[/(xo) -b 4/(xi) -b /(X 2 )] - (2.20) 

This yields the compound formula 


312 

Jdxf{x) 

XQ 


Xn Xq 
n 


'^WjfiXj) 

j=0 


1 (Xn - Xq)^ -^ 4 ^ 
180 re^ ^ 


( 2 . 21 ) 


where re is an even number, wo = Wn = 1/3, and for 1 < j < re we have Wj = 4/3 if j is 
odd and Wj if j = 2/3 is even. The error estimate scales now as l/re"^. 


As mentioned earlier, there are also rules involving non-equally spaced abscissas. A 
well known representative is the main formula of Gaussian quadrature. 

If w{x) is a weight function on [a,b], then there exist weights Wj and abscissas Xj for 
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1 < j < JT' such that 


with 


dxw{x) f{x) = '^Wjf{xj) + 




(2n)! 


b 

a 


( 2 . 22 ) 


n(x) = {x - Xl)(x - X2) . . . (x - Xn), 

a < xi < X 2 < ■ ■ ■ < Xn < b, a < ^ < b (2.23) 


The abscissas are given by the zeros of the orthogonal polynomial of degree n associated 
to the weight function w{x). In order to find them numerically, it is useful to know that 
they all lie in the interval [a, 6]. The weights are given by the (weighted) integral over 
the Lagrange polynomials: 

b 

Wj = j dxw{x)Vj{x) (2.24) 

a 


where the fundamental Lagrange polynomials are given by 

X _ (x - Xq) . . . (x - Xi_i)(x - Xj+i) . . . (x - Xn) 

* (Xj -Xq) ... (Xj - Xi_i)(Xj - Xj+i) . . . (Xj - Xn) 


(2.25) 


2.3.3 Monte Carlo techniques 

Monte Carlo integration is one of the two methods to perform the integrations in Chap¬ 
ters 7 and 8, although we also present only the results of CL This is common for the 
evaluation of multi-dimensional integrals. Generally speaking one is looking for an al¬ 
gorithm that meets a certain set of properties. It should 


• give a numerical estimate of the integral together with an estimate of the error, 

• yield the result in a reasonable amount of time, e.g. at low computational cost, 

• be able to handle multidimensional integrals. 


Monte Carlo integration delivers on all points of the list. In particular, as will be shown 
later, its error scales like 1/y/N, independent of the number of dimensions. This makes 
Monte Carlo integration the preferred method for integrals in high dimensions. 

Another integrator provided by the Cuba-library [71] is “VEGAS”. It is a Monte Carlo 
integrator paired with variance reducing techniques. These two concepts shall be shown 
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in the following subsection. 

Quadrature rules are inefficient for multidimensional integrals. Therefore one resorts to 
Monte Carlo integration. On the plus side, its error scales like l/^/N independent of the 
number of dimensions. But, on the other hand, a convergence by a rate of 1/\/iV is pretty 
slow. To improve efficiency, variance reducing techniques like importance sampling are 
employed. 

Consider the integral of a function /(ui,..., Ud)-, depending on d variables ui,... ,Ud over 
the unit hypercube [0,1]*^. Furthermore, / is assumed to be square-integrable. From 
now on, the short-hand notation x = {ui, ..., Ud) and for the function evaluated at that 
point f{x) = f{ui, ..., Ud) will be used. The Monte Carlo estimate for the integral 

I = j dxf{x) = j d^uf{ui, ...,Ud) (2.26) 

is given by 

1 ^ 

^ = (2.27) 

n=l 

The law of large numbers ensures that the Monte Carlo estimate converges to the true 
value of the integral: 

1 ^ 

j™ = P' 28 ) 

n=l 

The corresponding error S reads 

1 ^ 

= • (2-29) 

n=l 


2.3.3.1 Stratified sampling 


This technique consists of dividing the full integration space into subspaces, performing 
a Monte Carlo integration in each subspace, and adding up partial results in the end. 
Mathematically, this is based on the fundamental property of the Riemann integral 


J dxf{x) 
0 


a 1 

J dxf{x) + J dxf{x), 

0 a 


0 < a < 1. 


(2.30) 


More generally, one splits the integration region M = [0,1]'^ into k regions nij where 
j = 1,... ,k. In each region one performs a Monte Carlo integration with nj points. For 
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the integral /, one obtains the estimate 

P-31) 

j=l ^ n=l 

and instead of one has now the expression 


Y . AT. 52 

^ IV,- 

j=l 


(2.32) 


If the subspaces and the number of points in each subspace are chosen carefully, this 
can lead to a dramatic reduction in the error compared with crude Monte Carlo, but it 
should be noted, that it can also lead to a larger error if the choice is not appropriate. In 
general, the total variance is minimized when the number of points in each sub-volume 
is proportional to N ■ 8“^ . 


2.3.3.2 Importance sampling 


Mathematically, importance sampling corresponds to a change of integration variables: 

J dxf{x) = j ^p{x)dx = j ^dP{x) (2.33) 

with 

If one restricts p{x) to be a positive-valued function p{x) > 0 and to be normalized to 
unity 

J dxp{x) = 1 (2.35) 

one may interpret p{x) as a probability density function. If one has a random number 
generator corresponding to the distribution P{x) at his disposal, one may estimate the 
integral from a sample xi,..., xn of random numbers distributed according to P{x): 



n=l 


fix) 

p{x) 


The statistical error S of the Monte Carlo integration is given by 


52 



(2.36) 


(2.37) 
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It becomes evident that the relevant quantity is now f{x)/p{x) and it will be advan¬ 
tageous to choose p{x) as close in shape to f{x) as possible. In practice, one chooses 
p{x) such that it approximates |/(x)| reasonably well in shape and such that one can 
generate random numbers distributed according to P{x). 

One disadvantage of importance sampling is the fact, that it is dangerous to choose 
functions p{x), which become zero, or which approach zero quickly. If p goes to zero 
where / is not zero, S‘^{f/p) may be infinite and the usual technique of estimating the 
variance from the sample points may not detect this fact if the region where p = 0 is 
small. 


2.3.3.3 The VEGAS-algorithm 

The techniques described before require some advanced knowledge of the behaviour of 
the function to be integrated. In many cases this information is not available and one 
prefers adaptive techniques, e.g. an algorithm which learns about the functions as it 
proceeds. In the following, the VEGAS-algorithm will be presented. It combines the 
basic ideas of importance sampling and stratified sampling into an iterative algorithm, 
which automatically concentrates evaluations of the integrand in those regions where the 
integrand is largest in magnitude. VEGAS starts by subdividing the integration space 
into a rectangular grid and performs an integration in each subspace. These results 
are then used to adjust the grid for the next iteration according to where the integral 
receives dominant contributions. In this way VEGAS uses importance sampling and 
tries to approximate the optimal probability density function 

Poptimal(x) = (2-38) 

by a step function. Due to storage requirements one has to use a separable probability 
density function in d dimensions: 


p{ui,...,Ud) =Pi{ui) ■P2{u2) . Pd{ud)- (2.39) 

Eventually after a few iterations the optimal grid is found. In order to avoid rapid 
destabilizing changes in the grid, the adjustment of the grid includes usually a damping 
term. After this initial exploratory phase, the grid may be frozen and in a second 
evaluation phase the integral may be evaluated with high precision according to the 
optimized grid. The separation in an exploratory phase and an evaluation phase allows 
one to use less integrand evaluations in the first phase and to ignore the numerical 
estimates from this phase (which will in general have a larger variance). Each iteration 
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yields an estimate Ej together with an estimate for the variance 



(2.40) 


Here Nj denotes the number of integrand evaluations on iteration j. The results of each 
iteration on the evaluation phase are combined into a cumulative estimate, weighted by 
the number of calls Nj and their variances: 


E = 




SI j- 


(2.41) 


If the error estimates S’? become unreliable (for example if the function is not square 
integrable), it is more appropriate to weight the partial results by the number Nj of 
integrand evaluations alone. In addition VEGAS returns the per degree of freedom: 


xVdof 


1 ^ {E, - Ef 

m — I ^ 5^ 

j=i j 


(2.42) 


This allows a check whether the various estimates are consistent. One expects a x^/ dof 
not much greater than one. 
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Loop—Tree Duality at One—Loop 


This chapter serves to introduce the Loop-Tree Duality method (LTD) [72, 73] which 
we will be using subsequently throughout this thesis. It is a technique to numerically 
calculate multi-leg one-loop cross-sections in perturbative field theories. In fact, it is 
applicable in any quantum field theory in Minkowsky space with an arbitrary number of 
space-time dimensions. At its core, it establishes a relation between loop-level and tree- 
level amplitudes similar to the Feynman Tree Theorem (FTT) [74, 75]. Both methods 
allow to write basisc loop Feynman diagrams in terms of tree-level phase-space integrals 
which are obtained by cutting the original loop-integral. This is achieved by directly 
applying the Residue Theorem^, to the loop integrand. However, there are also some 
inportant differences between them: While the LTD produces only single cuts, the FTT 
also involves higher order cuts (double, triple, and so forth). 

In this chapter, we closely follow Refs. [48, 49] ^ to illustrate and derive the Loop-Tree 
Duality. 


^ Within the context of loop integrals, the use of the residue theorem has been considered many times 
in textbooks and in the literature. 

^ Ref. [48] is a good introduction to the Duality relation. It gives further insight into its nature and 
its comparison to the FTT. We refer the interested reader to that paper. For the purpose of this thesis 
we will follow it (and [49]) only as far as necessary to understand the later chapters. 
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3.1 Notation 


The FTT and the LTD can be illustrated with no loss of generality by considering 
their application to the basic ingredient of any one-loop Feynman diagrams, namely a 
generic one-loop scalar integral L^^'^ with N {N > 2) external legs. In the following and 
throughout this thesis, we assume the considered diagrams to be free from UV and IR 
divergencies. 

P2 


Figure 3.1: Momentum configuration of the one-loop iV-point scalar integral. 



The momenta of the external legs are denoted hy p^,P 2 , ■ ■ ■, p^ and are clockwise ordered 
(Fig. 3.1). All are taken as outgoing. To simplify the notation and the presentation, we 
also limit ourselves in the beginning to considering massless internal lines only. Thus, 
the one-loop integral can in general be expressed as; 


L^^Kpi,V2, ■ ■ ■,Pn) 



(3.1) 


where is the loop momentum (which flows anti-clockwise). The momenta of the 
internal lines are denoted by ; they are given by 


Qi = (! + '^Pk , (3.2) 

k=l 

and momentum conservation results in the constraint 

N 

= 0 . (3.3) 

i=l 

The value of the label i of the external momenta is defined modulo N, i.e. pN+i = Pi- 

The number of space-time dimensions is denoted by d (the convention for the Lorentz- 

indices adopted here is /r = 0,1,..., d—1) with metric tensor = diag(-|-l, —1,..., —1). 
The space-time coordinates of any momentum are denoted as = (fco,k), where 
ko is the energy (time component) of k^. It is also convenient to introduce light-cone 
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coordinates = (A;+, k_i_, k_), where k± = (/cq i kd-i)lV^- Throughout the chapter 


we consider loop integrals and phase-space integrals. If the integrals are ultraviolet or 
infrared divergent, we always assume that they are regularized by using analytic contin¬ 
uation in the number of space-time dimensions (dimensional regularization). Therefore, 
d is not fixed and does not necessarily have integer value. 

We introduce the following shorthand notation: 



(3.4) 


When we factorize off in a loop integral the integration over the momentum coordinate 
£o or £+, we write 



(3.5) 


and 



(3.6) 


respectively. The customary phase-space integral of a physical particle with momentum 
£ (i.e. an on-shell particle with positive-definite energy: = 0, £o > 0) reads 



(3.7) 


where we have defined 


6ii) = 27rie{£o) 5{£^) = 2Tri 5+i£^) . 


(3.8) 


Using this shorthand notation, the one-loop integral in Eq. (3.1) can be cast into 


N 




Pn) = / U Gpiqi) , 

^ 7 — 1 


(3.9) 


2 = 1 


where Gpiq) denotes the customary Feynman propagator, 



(3.10) 


We also introduce the advanced propagator GA{q) 



(3.11) 


We recall that the Feynman and advanced propagators only differ in the position of the 
particle poles in the complex plane (Fig. 3.2). Using q^ = q^ — = ‘2q+q- ~ ql) we 
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therefore have 


[G^(g)]-1=0 = 

/—z - q 1 — iO 

^ go = ±V*r lO ,OT q± = 

2g=F 

(3.12) 

[GA{q)]-^=0 = 

q2 

go — 5 or g± ~ + zO . 

zgzp 

(3.13) 


Thus, in the complex plane of the variable qo (or, equivalently^, q±), the pole with 
positive (negative) energy of the Feynman propagator is slightly displaced below (above) 
the real axis, while both poles (independently of the sign of the energy) of the advanced 
propagator are slightly displaced above the real axis. 




Figure 3.2; Location of the particle poles of the Feynman (left) and advanced (right) 
propagators, Gpiq) and Ga(<i), in the complex plane of the variable go or q±. 


3.2 The Feynman Tree Theorem 

In this section we briefly recall the FTT [74, 75]. To this end, we first introduce the 
advanced one-loop integral which is obtained from in Eq. (3.9) by replacing 
the Feynman propagators Gpiqi) with the corresponding advanced propagators GAiqi)'- 

r ^ 

^A^(Pi,P2,---,7'iv) = / Ga(%) . (3.14) 

i=i 

Then, we note that 

-^A^(pi, 7'2, •.. ,PAr) = 0 . (3.15) 

^ To be precise, each propagator leads to two poles in the plane qo and to only one pole in the plane 
q+ (or q-). 
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The proof of Eq. (3.15) can be carried out in an elementary way by using the Cauchy 
Residue Theorem and choosing a suitable integration path Cl- We have 


= / f dto ^ GA(gi) 

-h -) 

N 

= f f Mo =-2vri / ^ Res 

J J C T. ^ 1 ^ 


(3.16) 


{Im £o<0} 


N 


n GA{q^ 


. ^=1 


= 0 . 


The loop integral is evaluated by integrating first over the energy component Iq. Since 
the integrand is convergent when £o oo, the Iq integration can be performed along 
the contour Cl, which is closed at oo in the lower half-plane of the complex variable io 
(Fig. 3.3-left). The only singularities of the integrand with respect to the variable £o are 
the poles of the advanced propagators GA{qi), which are located in the upper half-plane. 
The integral along Cl is then equal to the sum of the residues at the poles in the lower 
half-plane and therefore it vanishes. 




Figure 3.3: Location of poles and integration contour Cl in the complex £o-plane 
for the advanced (left) and Feynman (right) one-loop integrals, and Lb). 


The advanced and Feynman propagators are related by 

CA{q)=GF{q) + S{q) , (3.17) 


which can straightforwardly be obtained by using the elementary identity 

= TiTrSix) , (3.18) 

X ± zU \x J 

where PV denotes the principal-value prescription. Inserting Eq. (3.17) into the right- 
hand side of Eq. (3.14) and collecting the contributions with an equal number of factors 
Gpiqi) and h{qj), we obtain a relation between and the one-loop integral 


r ^ 

La\pi,P2,. .. ,Pn) = / n [Cpiqi) + 6{qi) 

J £ .•_ 1 


2=1 

( 1 ) 


(3.19) 


= L'^^\pi,P2, . . . ,Pn) + L^ilc^.^{pi,P2, ...,Pn)^ -h ■ ■ • ,Pn) 
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Here, the single-cut contribution is given by 

-bS-cut(Pl>P2,...,Piv) = / S{qi) Gir(gj) . (3.20) 

i=l 3=1 

j¥=^ 

In general, the m-cut terms (m < N) are the contributions with precisely m 

delta functions S{qi): 

-^mLut(Pi)P 2 , ■ • ■,Pn) = j |5(gi). ..Siqm) Gpiqm+i) ■ ■ ■ GpiqN) + uneq. perms.j , 

(3.21) 

where the sum in the curly brackets includes all the permutations of gi,..., gw that give 
unequal terms in the integrand. 

Recalling that vanishes, cf. Eq. (3.15), Eq. (3.19) results in: 


LilcutiPi^P2 ,... ,Pw) + • • • + iyicut(^’i>?52, ■■■,Pn) ■ 

(3.22) 

This equation is the ETT in the specific case of the one-loop integral L^^\ The ETT 
relates the one-loop integral L^^l to the multiple-cut^ integrals Each delta func¬ 
tion S{qi) in replaces the corresponding Eeynman propagator in L^^^ by cutting 

the internal line with momentum g^. This is synonymous to setting the respective parti¬ 
cle on shell. An m-particle cut decomposes the one-loop diagram in m tree diagrams: in 
this sense, the ETT allows us to calculate loop-level diagrams from tree-level diagrams. 


L^^Hpi,P2, ■ ■ ■ ,Pn) = - 



Figure 3.4: The single-cut contribution of the Feynman Tree Theorem to the one- 
loop A-point scalar integral. Graphical representation as a sum of N basic single-cut 

phase-space integrals. 


In view of the discussion in the following sections, it is useful to consider the single-cut 

^ If the number of space-time dimensions is d, the right-hand side of Eq. (3.22) receives contributions 
only from the terms with m < d\ the terms with larger values of m vanish, since the corresponding 
number of delta functions in the integrand is larger than the number of integration variables. 
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contribution right-hand side of Eq. (3.22). In the case of single-cut contri¬ 

butions, the FTT replaces the one-loop integral with the customary one-particle phase- 
space integral, see Eqs. (3.7) and (3.20). Using the invariance of the loop-integration 
measure under translations of the loop momentum we can perform the momentum shift 
£—>•£ — Yl\=iPk ill the term proportional to 5[qi) on the right-hand side of Eq. (3.20). 
Thus, cf. Eq. (3.2), we have qi ^ i and qj ^ i + {Pi+i + Pi +2 + • • • + Pi+j), with i / j. 
We can repeat the same shift for each of the terms {i = 1,2,, N) in the sum on the 
right-hand side of Eq. (3.20), and we can rewrite as a sum of N basic phase-space 

integrals (Eig. 3.4): 

^i-cut(7'i>7’2, • • ■ ,Pn) = ...,pi+P 2 -\ -hpAT-i) -h cyclic perms. 

N 

= ^ l[^~,^^{pi,pi +Pi+i,... ,pi+pi+i H- \-pi+N-2) (-3.23) 

i=l 

We denote the basic one-particle phase-space integrals with n Feynman propagators by 
They are defined as follows: 

= YiaHi+k,) = lm 

j — 1 j — 1 ^ 


The extension of the FTT from the one-loop integrals to one-loop scattering am¬ 
plitudes (or Green’s functions) in perturbative field theories is straightfor¬ 

ward, provided the corresponding field theory is unitary and loeal. The generalization 
of Eq. (3.22) to arbitrary scattering amplitudes is [74, 75]: 


^(l-loop) 


,(l-loop) .(l-loop) 
•^l-cut ' “^2-cut ' 


(3.25) 


where obtained in the same way as i.e. by starting from 

and considering all possible replacements of m Feynman propagators GF{qi) of its loop 

internal lines with the ‘cut propagators’ 6{qi). 

The proof of Eq. (3.25) directly follows from Eq. (3.22): ^(1“1 °°p) is a linear combination 
of one-loop integrals that differ from only by the inclusion of interaction vertices. 
As briefly recalled below, this difference has harmless consequences on the derivation of 
the FTT. 


Including particle masses in the advanced and Feynman propagators has an effect on 
the location of the poles produced by the internal lines in the loop. However, as long as 
the masses are real, as in the case of unitary theories, the position of the poles in the 
complex plane of the variable go is affected only by a translation parallel to the real axis. 
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with no effect on the imaginary part of the poles. This translation does not interfere 
with the proof of the FTT as given in Eqs. (3.14)-(3.22). Therefore, the effect of a 
particle mass Mj in a loop internal line with momentum qi simply amounts to modifying 
the corresponding on-shell delta function 5{qi) when this line is cut to obtain 
This modification then leads to the obvious replacement: 

= 2vri6l(gi,o) - m^) = 27rz - m-) • (3.26) 

Including interaction vertices has the effect of introducing numerator factors in the 
integrand of the one-loop integrals. As long as the theory is local, these numerator factors 
are at worst polynomials of the integration momentum I ^ . In the complex plane of 
the variable (.q, this polynomial behavior does not lead to additional singularities at any 
finite values of ^o- The only danger, when using the Cauchy theorem as in Eq. (3.16) to 
prove the ETT, stems from polynomials of high degree that can spoil the convergence of 
the £o-iiit6grcition at infinity. Nonetheless, if the field theory is unitary, these singularities 
at infinity never occur since the degree of the polynomials in the various integrands is 
always sufficiently limited by the unitarity constraint. 


3.3 The Loop—Tree Duality 


In this Section we derive and illustrate the Loop-Tree Duality relation between one-loop 
integrals and single-cut phase-space integrals. This relation is the main general result 
of this chapter. 

Rather than starting from we directly apply the Residue Theorem to the compu¬ 
tation of lT). We proceed exactly as in Eq. (3.16), and obtain 


N 


L^^Hpi,P2,---,Pn)= I dio G'f’(gi) 


2=1 



N 


dio Tf GF(gi) =/ V Res{im4<o} 
I JCl Jt 


N 


n GF{q^) 


_ ^=1 


.(3.27) 


At variance with G^(gj), each of the Eeynman propagators Gpigi) has single poles in 
both the upper and lower half-planes of the complex variable Iq (see Eig. 3.3-right) and 
therefore the integral does not vanish as in the case of the advanced propagators. In 
contrast, here, the N poles in the lower half-plane contribute to the residues in Eq. (3.27). 

®This statement is not completely true in the case of gauge theories and, in particular, in the case of 
gauge-dependent quantities. For the discussion of the additional issues that arise in gauge theories see 
Ref. [48] Section 9. 
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The calculation of these residues is elementary, but it involves several subtleties. The 
detailed calculation, including a discussion of its subtle points, is presented in the fol¬ 
lowing. 

The sum over residues in Eq. (3.27) receives contributions from N terms, namely the N 
residues at the poles with negative imaginary part of each of the propagators Gpiqi), 
with i = 1,... ,N, see Eq. (3.12). Considering the residue at the z-th pole we write 



N 


N 

R,6S{2_th pole} 

n Griqj) 

. 

[Res|j_^,jj pole} Gp{qi)~^ 

1 

1 _ 


(3.28) 

where we have used the fact that the propagators Gp^qj)-, with j ^ i, are not singular at 
the value of the pole of Gpiqi)- Therefore, they can be directly evaluated at this value. 

Applying the Residue Theorem in the complex plane of the variable io, the computa¬ 
tion of the one-loop integral reduces to the evaluation of the residues at N poles, 
according to Eqs. (3.27) and (3.28). 

The evaluation of the residues in Eq. (3.28) is a key point in the derivation of the Loop- 
Tree Duality relation. To make this point as clear as possible, we first introduce the 
notation ( 7 ^^^ to explicitly denote the location of the i-th pole, i.e. the location of the 
pole with negative imaginary part (see Eq. (3.12)) that is produced by the propagator 
Gpiqi)- We further simplify our notation with respect to the explicit dependence on 
the subscripts that label the momenta. We write Gpiqj) = Gpiqi + {qj — qi)), where 
qi depends on the loop momentum while {qj — qi) = kji is a linear combination of the 
external momenta (see Eq. (3.2)). Therefore, to carry out the explicit computation of 
the Tth residue in Eq. (3.28), we re-label the momenta hy qi ^ q and qj ^ q + kj, and 
we simply evaluate the term 


where (see Eq. (3.12)) 

= Vq2 - iO . (3.30) 

In the next paragraphs, separately compute the residue of Gpiq) and its pre-factor - 
the associated factor arising from the propagators GF{q + kj). 


Gpiq + kj 


(3.29) 
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The computation of the residue of Gpiq) gives 






lim 
—>■ ( 

1 




— + iO 


2 9o 


{+) 


= / dqQ5+{q^) , 


(3.31) 


thus leading to the result in Eq. (3.32). Note that the first equality in the second line of 
Eq. (3.31) is obtained by removing the iO prescription from the previous expression. This 
is fully justified. The term (^q"*^^)”^ = {y /— f0)“^ becomes singular when q^ —)■ 0, 
and this corresponds to an end-point singularity in the integration over q: therefore the 
iO prescription has no regularization effect on such end-point singularity. The second 
equality in the second line of Eq. (3.31) simply follows from the definition of the on-shell 
delta function S+{q‘^), which we will later call the dual delta function. 


Hence, the calculation of the residue of Gpiqi) gives 

= I dio S+{q^) . (3.32) 

This result shows that considering the residue of the Eeynman propagator of the internal 
line with momentum qi is equivalent to cutting that line by including the corresponding 
on-shell propagator 6+{qf). The subscript -|- of 5+ refers to the on-shell mode with 
positive definite energy, qi^ = the positive-energy mode is selected by the Feynman 

iO prescription of the propagator Gpiqi)- The insertion of Eq. (3.32) in Eq. (3.27) directly 
leads to a representation of the one-loop integral as a linear combination of N single-cut 
phase-space integrals. 


[l^es|j_^jj pole} Gpi^qi')^ 


Resr,_ 


{i—thpolej 2 I -rv 

^ qf +l0 


We now consider the evaluation of the residue pre-factor (the second square-bracket 
factor in Eq. (3.29)). We first recall that the iO prescription of the Feynman propagators 
has played an important role in the application (see Eqs. (3.27) and (3.29)) of the Residue 
Theorem to the computation of the loop integral: having selected the pole with negative 
imaginary part, q^ = qQ^\ the prescription eventually singled out the on-shell mode with 
positive definite energy, go = (see Eq. (3.31)). However, the evaluation of the one- 
loop integrals by the direct application of the Residue Theorem (as in Eq. (3.27)) involves 
some subtleties. The subtleties mainly concern the correct treatment of the Feynman iO 
prescription in the calculation of the residue pre-factors. A consistent treatment requires 
the strict computation of the residue pre-factor in Eq. (3.29): the iO prescription in both 
Gf{q + kj) and has to be dealt with by considering the imaginary part iO as a 
finite, though possibly small, quantity; the limit of infinitesimal values of iO has to be 
taken only at the very end of the computation, thus leading to the interpretation of the 
ensuing iO prescription as mathematical distribution. Applying this strict procedure, we 
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obtain 


Gpiq + kj) 
j 


9o='!'o^^ 



1 


{q + kjY + iO 


90 = 90 "^’ 


n 


1 

2 gQ^^/cjo - 2q • kj + /c? 


n 


1 

2 |q|A:jo - 2q • kj + kj - iOA;jo/|q| 



1 

2qkj + kj - iOkjo/qo 


9o=|q| 


.(3.33) 


The last equality on the first line of Eq. (3.33) simply follows from setting Qq = 
in the expression on the square-bracket (note, in particular, that q^ = —iO). The first 
equality on the second line follows from ~ 2|q| — fO/lql (i.e. from expanding 
at small values of iO). 


The result in Eq. (3.33) for the residue pre-factor is well-defined and leads to a well- 
defined (i.e. non singular) expression once it is inserted in Eq. (3.27). The possible 
singularities from each of the propagators l/{q + kj)'^ are regularized by the displace¬ 
ment produced by the associated imaginary amount iOkjo/qQ. Performing the limit of 
infinitesimal values of iO, only the sign of the iO prescription (and not its actual mag¬ 
nitude) is relevant. Therefore, since qo is positive, in Eq. (3.33) we can perform the 
replacement iOkjo/qo —)• iOrjkj, where is a future-like vector with 


Vu = iVo,r]) , % > 0, 7/^ = > 0 , (3.34) 


i.e. a d-dimensional vector that can be either light-like (r/^ = 0) or time-like ( 77 ^ > 0) 
with positive definite energy rjQ. Hence, we finally obtain 

YlGpiq + k,) 

j 

After reintroducing the original labels of the momenta of the loop integral according to 
the replacements q —)• qi, kj —)> qj — qi, see the discussion above Eq. (3.29), the calculation 
of the residue pre-factor on the r.h.s. of Eq. (3.28) yields 


n 


{q -\- kjY — iO'qkj 


(3.35) 


n = lnl 


n 


{i—th pole} 


n 


. q'j +i0 


{2—th pole} 


n 


1 

iOrjiqj - qi) 


(3.36) 


Note that the calculation of the residue at the pole of the internal line with momen¬ 
tum qi changes the propagators of the other lines in the loop integral. Although the 
propagator of the j-th internal line still has the customary form l/^j, its singularity at 
qj = 0 is regularized by a different iO prescription: the original Eeynman prescription 
qj + iO is modified in the new prescription qj — i0r]{qj — qi), which we name the ‘dual’ 
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iO prescription or, briefly, the rj prescription. The dual iO prescription arises from the 
fact that the original Feynman propagator 1/(^| + i0) is evaluated at the comp/ex value 
of the loop momentum q, which is determined by the location of the pole at qf + i0 = 0. 
The iO dependence from the pole has to be combined with the iO dependence in the 
Feynman propagator to obtain the total dependence as given by the dual iO prescrip¬ 
tion. The presence of the vector rj^ is a consequence of using the Residue Theorem. To 
apply it to the calculation of the d dimensional loop integral, we have to specify a system 
of coordinates (e.g. space-time or light-cone coordinates) and select one of them to be 
integrated over at fixed values of the remaining d—\ coordinates. Introducing the aux¬ 
iliary vector rj^ with space-time coordinates = {r]o,0±,r]d-i), the selected system of 
coordinates can be denoted in a Lorentz-invariant form. Applying the Residue Theorem 
in the complex plane of the variable qo at fixed (and real) values of the coordinates q_L 
and q'^_^ = qu-i — qor]d-i/vo (to be precise, in Eq. (3.27) we actually used = (1,0)), 
we obtain the result in Eq. (3.36). 


The rj dependence of the ensuing iO prescription is thus a consequence of the fact that the 
residues at each of the poles are not Lorentz-invariant quantities. The Lorentz-invariance 
of the loop integral is recovered only after summing over all the residues. 



Pi-i Pi 



Figure 3.5: The Duality relation for the one-loop A-point scalar integral. Graphical 
representation as a sum of N basic dual integrals. 


Inserting the results of Eq. (3.28)-(3.36) in Eq. (3.27) we directly obtain the Duality 
relation between one-loop integrals and phase-space integrals: 


L^^Kpi,P2, ...,Pn) = - L^^\pi,P2, ... ,pn) , 


where the explicit expression of the phase-space integral is (Fig. 3.5) 
L^^Hpi,P 2 ,---,Pn) = / V S{qi) n Goiqpqj) , 

i=i 


t=i 


I 


with Gniqi] qj) = —tetA 


— iO pk 




(3.37) 


(3.38) 

(3.39) 
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where r] is the auxiliary vector defined in Eq. (3.34). Each of the N — 1 propagators 
in the integrand is regularized by the dual iO prescription and, thus, it is named ‘dual’ 
propagator. Note that the momentum difference qj — qi is independent of the integration 
momentum l\ it only depends on the momenta of the external legs (see Eq. (3.2)). 

Using the invariance of the integration measure under translations of the momentum 
we can perform the same momentum shifts as described in Section 3.2. In analogy to 
Eq. (3.23), we can rewrite Eq. (3.38) as a sum of iV basic phase-space integrals (Fig. 3.5): 


L^^\pi,P2,. .. ,pn) = ^\pi ,Pi + P 2 , ■ ■ ■ ,Pi + P 2 -\ - \-PN-i) + cyclic perms. 

N 

= -bPt+i,... -Fpi+I H-hpj+Ar- 2 ) .(3.40) 

i=l 

The basic one-particle phase-space integrals with n dual propagators are denoted by 
and are defined as follows: 

(3.41) 

We now comment on the comparison between the FTT (Eqs. (3.20)-(3.24)) and the 
Loop-Tree Duality (Eqs. (3.37)-(3.41)). The multiple-cut contributions L^^cut) with 
m > 2, of the FTT are completely absent from the Loop-Tree Duality relation, which 
only involves single-cut contributions similar to those in However, the Feynman 

propagators present in are replaced by dual propagators in This compen¬ 

sates for the absence of multiple-cut contributions in the Loop-Tree Duality. 

The iO prescription of the dual propagator depends on the auxiliary vector r]. The 
basic dual contributions 1*'”^ are well defined for arbitrary values of rj. However, when 
computing the future-like vector rj has to be the same across all dual contributions 
(propagators): only then does not depend on r]. 

In our derivation of the Loop-Tree Duality relation, the auxiliary vector r/ originates from 
the use of the Residue Theorem. Independently of its origin, we can comment on the 
role of rj in the Duality relation. The one-loop integral {pi,p2, ■ ■ ■ ,Pn) is a function 
of the Lorentz-invariants (piPj)- This function has a complicated analytic structure, 
with pole and branch-cut singularities (scattering singularities), in the multidimensional 
space of the complex variables (piPj)- The iO prescription of the Feynman propagators 
selects a Riemann sheet in this multidimensional space and, thus, it unambiguously 
defines L^^'){pi,p2, ■ ■. ,Pn) as a single-valued function. Each single-cut contribution to 
has additional (unphysical) singularities in the multidimensional complex space. 
The dual iO prescription fixes the position of these singularities. The auxiliary vector 
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r] correlates the various single-cut contributions in L^^\ so that they are evaluated on 
the same Riemann sheet: this leads to the cancellation of the unphysical single-cut 
singularities. In contrast, in the FTT, this cancellation is produced by the introduction 
of the multiple-cut contributions R^^cut- 

We remark that the expression (3.40) of L^^^ as a sum of dual contributions is just a 
matter of notation: for massless internal particles is actually a single phase-space 
integral whose integrand is the sum of the terms obtained by cutting each of the internal 
lines of the loop. In explicit form, we can write: 

L^^\p1,P2,---,Pn) = ^'>{e]Pi,Pi+Pi+l,...,Pi+Pi+l-\ - \-pi+N-2) , 

i=l 

(3.42) 

where the function is the integrand of the dual contribution in Eq. (3.41). There¬ 
fore, the Loop-Tree Duality relation (3.37) directly expresses the one-loop integral as the 
phase-space integral of a tree-level quantity. To name Eq. (3.37), we have introduced the 
term ‘duality’ precisely to point out this direct relation^ between the d-dimensional inte¬ 
gral over the loop momentum and the (d — l)-dimensional integral over the one-particle 
phase-space. Eor the ETT, the relation between loop-level and tree-level quantities is 
more involved, since the multiple-cut contributions (with m > 2) contain inte¬ 

grals of expressions that correspond to the product of m tree-level diagrams over the 
phase-space for different number of particles. 

The simpler correspondence between loops and trees in the context of the Loop-Tree 
Duality relation is further exploited in Ref. [73] Sect. 10, where the Green’s functions 
and scattering amplitudes are discussed. 


3.4 Explicit example: The scalar two-point function 

In this Section we illustrate the application of the Loop-Tree Duality relation to the 
evaluation of the one-loop two-point function L^‘^\ A detailed discussion (including a 
comparison between ETT and Loop-Tree Duality as well as detailed results in analytic 
form and numerical results) of higher-point functions is found in [48] (see also Ref. [72]). 

The two-point function (Eig. 3.6), also known as bubble function Bub, is the simplest 
non-trivial one-loop integral with massless internal lines: 

Bub(rf) ^ --i [P + iO]|(/+KP + ffll ■ 

^The word duality also suggests a stronger (possibly one-to-one) correspondence between dual inte¬ 
grals and loop integrals, which is further discussed in Ref. [73] Sect. 7. 
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i + Pi 



Figure 3.6: The one-loop two-point scalar integral L^^^(pi,p 2 ). 


Here, we have visibly implemented momentum conservation {pi -\-p 2 = 0) and exploited 
Lorentz invariance {L^‘^\pi,p 2 ) can only depend on pf, which is the sole available invari¬ 
ant). Since most of the one-loop calculations have been carried out in four-dimensional 
field theories (or in their dimensionally-regularized versions), we set d = 4 — 2e. Note, 
however, that we present results for arbitrary values of e or, equivalently, for any value 
d of space-time dimensions. 


The result of the one-loop integral in Eq. (3.43) is well known: 

Bub(p2) = cr {-P^ - iO)"" 


(3.44) 


where cp is the customary d-dimensional volume factor that appears from the calculation 
of one-loop integrals: 


r(l + e) r^(i-e) 
(47r)2""r(l-2e) 


(3.45) 


We recall that the it) prescription in Eq. (3.44) follows from the corresponding prescrip¬ 
tion of the Feynman propagators in the integrand of Eq. (3.43). The iO prescription 
defines Bub(p^) as a single-value function of the real variable p^. In particular, it gives 
Bub(p^) an imaginary part with an unambiguous value when p^ > 0: 

Bub(p^) = cr (Ip^I)"" [6i(-p^)+ 6»(/) . (3.46) 


3.4.1 General form of single-cut integrals 


To apply the Loop-Tree Duality relation, we have to compute the single-cut integrals 
-^i-cut respectively. Since these integrals only differ because of their it) prescrip¬ 

tion, we introduce a more general regularized version, dreg, of the single-cut integral. We 
define: 



1 

2dA; -\-hf + i0c(/c) 


d^d 

(27r)'^“^ 




1 

2(.k + kf + i0c{k) 

(3.47) 
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Although c{k) is an arbitrary function of k, L^e^ only depends on the sign of the iO 
prescription, i.e. on the sign of the function c{k): setting c{k) = +1 we recover 
cf. Eq. (3.24), while setting c{k) = —rjk we recover (see Eq. (3.41)). 

The calculation of the integral in Eq. (3.47) is elementary, and the result is 


lgl{k-c{k)) = -, 


cr 


1 


^ - iQk^ c{k) 

TO 


[kQ-it)k^c{k)] " . (3.48) 


2cos(7re) e(l — 2e) 

Note that the typical volume factor, cr, of the d-dimensional phase-space integral is 

r(l-e) r(l + 2e) 


cr = 


(4vr) 


2-e 


(3.49) 


The factor cos(7re) in Eq. (3.48) originates from the difference between cr and the volume 
factor Cr of the loop integral: 


cr _ r(l-h2e)r(l-2e) _ 1 

Cr r(l-|-e) r(l — e) cos(7re) 


We also note that the result in Eq. (3.48) depends on the sign of the energy /cq. This 
follows from the fact that the integration measure in Eq. (3.47) has support on the 
future light-cone, which is selected by the positive-energy requirement of the on-shell 
constraint (5+(^^). 

The denominator contribution {2lk+k‘^) in the integrand of Eq. (3.47) is positive definite 
in the kinematical region where > 0 and k^ > 0. In this region the iO prescription 
is inconsequential, and Ireg has no imaginary part. Outside this kinematical region, 
{2ik + k"^) can vanish, leading to a singularity of the integrand. The singularity is 
regularized by the it) prescription, which also produces a non-vanishing imaginary part. 
The result in Eq. (3.48) explicitly shows these expected features, since it can be rewritten 
as 


lgl{k;c{k)) = 


+ 


[cos(7re) - i sin(7re) sign(c(fc))] 

2cos(7rej e(i — 2ej 

6{k^) [0{ko) + 9{—ko) (cos(27re) -|- i sin(27re) sign(c(fc)))]} . (3.51) 


We note that the functions Bub(/c^) and /reg(/c; c(/c)) have different analyticity proper¬ 
ties in the complex k'^ plane. The bubble function has a branch-cut singularity along 
the positive real axis, k'^ > 0. The phase-space integral Lreg{k] c{k)) has a branch-cut 
singularity along the entire real axis if fcg < 0, while the branch-cut singularity is placed 
along the negative real axis if A:o >0. 
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3.4.2 Duality relation for the two-point function 



Figure 3.7: 


1 






One-loop two-point function: the Duality relation. 


We now consider the Loop-Tree Duality (Fig. 3.7) in the context of this example. The 
dual representation of the one-loop two-point function is given by 


L^‘^\piiP 2) = I^^\pi) + {pi -pi) , (3.52) 

cf. Eqs. (3.40) and (3.41). The dual contribution I^^\k) is obtained by setting c{k) = 
—rjk in Eq. (3.48). Since is a future-like vector, c{k) has the following important 
property: 

sign(r//c) = sign(/co) , if > 0 . (3.53) 

Using this property, the result in Eq. (3.48) can be written as 


/W(A;) 


cr {-k"^ - ^O) 
y e(l-2e) 


1 - i sign{k^rik) 

cos(7re) 


(3.54) 


Comparing this expression with Eq. (3.44), we see that the imaginary contribution in the 
square bracket is responsible for the difference with the two-point function. However, 
since sign(—r/fc) = —sign(? 7 fc), this contribution is odd under the exchange k —)• —k and, 
therefore, it cancels when Eq. (3.54) is inserted in Eq. (3.52). Taken together, 

L^‘^\pi,P2) = I^^\pi) + (pi ^ -Pi) = - cr ^ j ’ (3-55) 


which fully agrees with the Duality relation L^‘^\pi,p 2 ) 


Bub(pf). 


3.5 Loop—Tree Duality with generic masses 


This section serves to recapitulate the chapter and give the explicit formulae for generic 
masses as well as some explicit example calculations to illustrate how the Loop-Tree 
Duality works in practice. 
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3.5.1 Generic masses 


The Feynman propagators Gpiqi) in Eq. (3.1) with real internal masses mi read: 


Gpiqi) 


1 

qf - m? + iO 


(3.56) 


The derivation of the Loop-Tree Duality Theorem is exactly the same regardless of 
the internal lines being massive or massless (m* = 0), as long as the masses are real. 
Non-vanishing internal real masses only account for a displacement of the poles of the 
propagators along the real axis, which does not change the derivation of the Duality 
Theorem, as will become obvious in the following. Moreover, they do not alter the rela¬ 
tionship between Feynman, advanced, retarded and dual propagators, which is the basis 
of both, the Duality Theorem as a duality to the FTT. 

Besides the customary Feynman propagators Gpiqt), we also encounter advanced, GAiqt), 
and retarded, GR{qi), propagators, defined by: 


GA{q^) 



Gdqi) 


1 

qf - m‘f + iO qifl 


(3.57) 


The Feynman, advanced, and retarded propagators only differ in the position of the 
particle poles in the complex plane. Using g? = g?Q — q?, we therefore find the poles of 
the Feynman and advanced propagators in the complex plane of the variable qi^ at: 


[GF{qi)]~^ = 0 qifl = - mf - iO and [GA{qi)]~^ = 0 qifl ~ i-^/qf - mf+iO . 

(3.58) 

Thus, the pole with positive/negative energy of the Feynman propagator is slightly dis¬ 
placed below/above the real axis, while both poles of the advanced/retarded propagator, 
independently of the sign of the energy, are slightly displaced above/below the real axis 
(cf. Fig. 3.2). Similarly to the massless case, we further define 


6{qi) = 2-71 iB{qi^o) 6{qf-m‘f) = 2TTi6+{qf-mj) , (3.59) 


where again the subscript -|- of 5+ refers to the on-shell mode with positive definite 
energy, qi^ > 0. Hence, the phase-space integral of a physical particle with momentum 
qi, i.e., an on-shell particle with positive-definite energy, qf = mf, qifl > 0, reads: 

/ (2^)d-i ~ ••• • (3.60) 

In order to derive the Duality Theorem, one directly applies the Residue Theorem to 
the computation of L^^\pi,p 2 , ■.. ,Pn) in Eq. (3.1): Each of the Feynman propagators 
Gpiqi) has single poles in both the upper and lower half-planes of the complex variable 
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3 . 5 . Loop-Tree Duality with generic masses 


io- Since the integrand is convergent when oo, by closing the contour at oo in the 
lower half-plane and applying the Cauchy theorem, the one-loop integral becomes the 
sum of N contributions, each of them obtained by evaluating the loop integral at the 
residues of the poles with negative imaginary part belonging to the propagators Gpiqi)- 
The calculation of the residue of Gpiqi) gives 

Res[GF(gi)]im(g,,o)<o = j dio S+igf - mf) , (3.61) 

with — m?) defined in Eq. (3.59). This result shows that considering the residue 

of the Feynman propagator of the internal line with momentum qi is equivalent to 
cutting that line by including the corresponding on-shell propagator 5+(q'f — rnf). The 
propagators Gpiqj), with j ^ i, are not singular at the value of the pole of Gpiqi) and 
can therefore be directly evaluated at this point, yielding to 


n ^Fiqj) 


Gp{qi)-^=0 


n Gniq^qj) , 


(3.62) 


where 


GD{qi;qj) 


1 

qj -mj -iOr] ■ kji 


(3.63) 


is (massive) dual propagator with r/ a future-like vector, defined as in (3.34). i.e., a 
d-dimensional vector that can be either light-like (r/^ = 0) or time-like (t/^ > 0) with 
positive dehnite energy r/o and kji = qj — qi. Collecting the results from Eq. (3.61) and 
Eq. (3.62), the Loop-Tree Duality Theorem at one-loop takes the hnal form 


L^^\pi,P2,... ,pn) 



N 

d{qi) n Gniq^qj) ■ 

i=i 


(3.64) 


3.5.2 Explicit example 

Here we present a very simple yet fully featured example of a triangle graph with generic 
internal masses. The integral to calculate is 


L^^\pi,P2,P3) = I GF{qi)GF{q2)GF{q3) , 


(3.65) 


with the three Feynman propagators 


GF{qi) = 


1 


gj — mf + iO 


Tx > GF{q2) = 


1 


q2 — m2+ it) 


and GF{q 3 ) = 


1 


§3 — TTlg + iO 


(3.66) 
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where qi = i+pi, q 2 = ^+Pi+P 2 = ^ and qs = i according to the conventions established 
in Section 3.1. Applying Duality means using Eq. (3.38) to rewrite the integral I. A 
triangle has three internal lines, thus N = 3. Integral I takes the form: 


,P2,P3)= / S{Qi)GDiqi]q2)GD{qi]q3) 

Je 

+ / GD{q2]qi)Siq2)GD{q2]q3) 


+ GD{q3]qi)GD{q3]q2)Hq3) 


first contribution 
second contribution 


third contribution (3.67) 


Each summand of Eq. (3.38) is a dual contribution. In the first dual contribution the line 
carrying qi gets cut, i.e. it becomes a dual delta function, i.e. 6{qi) = 27ri 9{qifi) S{qi). 
The Eeynman propagators that correspond to the other internal line get promoted to a 
dual propagators as in Eq. (3.39). The first argument of the dual propagator indicates 
the cut, the second one the momentum of the internal line it is assigned to. 

In order to produce dual contribution two, the next internal line (the one associated 
with 52 ) gets cut, and all of the other lines (= lines carrying qi and ^ 3 ) converted to 
dual propagators, similarly to the first dual contribution. 

To evaluate the dual delta functions we take advantage of 

= E . (3.68) 


where Xi are the zeros of g{x). Hence the dual deltas yield: 

~ _ (5(4 - i-Pi,o + \/(€ + Pi)2 + m^)) 

2Vie + Pir + mj) 

^ <^(4 - j-Pifi - P2,o + Vii + Pi+ P2)^ + ml)) 

2 y^(£ + Pi + P 2)2 + ml) 

<5(4 - + ml) 

< 5 ( 93 ) =- , - • 

2w4 + ml 


(3.69) 


There is a crucial difference to note: The dual delta functions <5(qi), ( 5 (^ 2 ) and < 5 (^ 3 ) 
force the zero component of the loop integration to different values. For example, in 
contribution one we have q[~\l = -\Aii ^0 = —Pi,o + \/(-^ + Pi)^ + mf whereas in 
contribution three the zero component is fixed to => 4 = \l+ ml- 

This will heavily affect the structure of the dual contributions. Now we can insert for 
the dual propagators and apply the dual delta functions. To give the reader a better 
idea of how the outcome looks we write down contribution three, which we will call I 3 , 
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explicitly: 

h = - [ - , - -^- ^ - • (3.70) 

2pifiJ(? + + 2£ • Pi — m\ + + Pi — iOv^is 2Jl^+ rri^ 

1 

2(pi,o +P 2 ,o)-\/^^ + "I 3 + 2^ • (Pi + P 2 ) + (pi +P 2 )^ - - i0r]k23 

The other two dual contributions look very similar. This is the final result that we 
obtain from using the Loop-Tree Duality. We see that it is a phase-space integration 
which runs only over the loop three momenta. We could now put numbers for pi and p 2 
and perform an actual calculation. Therefore we would pass the integral to a numerical 
integrator as we will do in Chapter 7. In fact, Section 7.7.2 shows the results for this 
exact triangle that we have discussed here. 
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Chapter 4 


Loop—Tree Duality beyond 
One—Loop 


In the previous chapter, the Loop-Tree Duality method was introduced for the one-loop 
case. It has the appealing property of recasting the virtual corrections in a from which 
is very similar to the real ones, thus giving rise to the idea of directly combining the 
two. In the introduction, this Duality relation has been derived for the one-loop case. 
The purpose of this chapter is twofold. 

First we want to establish a better-suited notation which then in turn allows the general¬ 
ization of the Duality Theorem to situations involving two-loops. Therefore, a systematic 
procedure is presented which can be employed repeatedly to cover diagrams with even 
more loops. For the scope of this thesis it is sufficient to illustrate the technique at the 
two-loop level. Readers who are interested in applying it beyond two-lops are advised 
to check Ref. [49]. 

Second, this chapter helps to prepare to following one, in which we will be dealing with 
higher order poles. At the one-loop level, these can be avoided through an adequate 
choice of gauge, however, at two-loops this no longer the case. 
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4.1 Duality relation at one—loop 


It was shown in Chapter 3 that using the Cauchy residue theorem the one-loop integral 
can be written in the form: 


N 


L’'^\pi,P2,---,Pn) = / S{qi) Yl Gniquqj) 


(4.1) 


j=i 


where 


GD{qi]qj) = 


~ - qi) ’ 

is the so-called dual propagator, as defined in Ref. [48] , with rj a future-like vector 


(4.2) 


Vo >0,V^ = VuP^ > 0 , 


(4.3) 


The extension of the Duality Theorem to two loops has been discussed in detail in 
[49, 76]. Here we recall the basic points. We extend the definition of propagators of 
single momenta to combinations of propagators of sets of internal momenta. Let ak be 
any set of internal momenta qi,qj with i,j G ak- We then define Feynman and dual 
propagator functions of this set ak in the following way: 

Gpiak) = n Gpiqi) , Goiak) = 5iqi) Goimiqj) - (4.4) 

ieak i&Q-k j&cik 

By definition, Goiak) = d{qi), when ak = {*} and thus consists of a single four momen¬ 
tum. At one-loop order, ak is naturally given by all internal momenta of the diagram 
which depend on the single integration loop momentum I'l, ak = {1,2,..., A^}. How¬ 
ever, let us stress that ak can in principle be any set of internal momenta. At higher 
order loops, e.g., several integration loop momenta are needed, and we can define several 
loop lines ak to label all the internal momenta (cf. Eq. (4.14)) where Eq. (4.4) will be 
used for these loop lines or unifications of these. We also define: 

Gpi-ak) = ^(“5*) n Goi-qu-qj) , (4.5) 

i&ctk jeoik 

where the sign in front of ak indicates that we have reversed the momentum flow of 
all the internal lines in ak- For Feynman propagators, moreover, Gpi—ak) = Gpio^k)- 
Using this notation the following relation holds for any set of internal momenta a^: 

GA{o^k) = Gpiak) + Gpiak) , (4.6) 
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where GA{qi) is the advanced propagator: 


GA{qi) 


1 

-rnf-it) qifl ’ 


(4.7) 


and 

GA(afe) = n GA{qi) . (4.8) 

i&ak 

The proof of Eq. (4.6) can be found in Ref. [49]. Note that individual terms in Goiotk) 
depend on the dual vector p, but the sum over all terms contributing to Goictk) is 
independent of it. Another crucial relation for the following is given by a formula that 
allows to express the dual function of a set of momenta in terms of chosen subsets. 
Considering the following set Pn = ai U ... U a^, where f3]sf is the unification of various 
subsets Oj, we can obtain the relation: 


Gd{(^1 U 0(2 U ... U = E n GD{a^,) n GHai,). (4.9) 

The sum runs over all partitions oi Pn into exactly two blocks and with elements 

Ui, i G {1,..., N}, where, contrary to the usual definition, we include the case: = Pn, 

(‘ 2 '\ 

(3)^' = 0. For the case of A" = 2, e.g., where /32 = ai U 0 : 2 , we have: 


Goiai U 02 ) = Goioii) GD{oi2) + Gd{oii) Gi?(a2) + GF{ai) GD{oi2) ■ (4.10) 


Naturally it holds that: 


N 

GplyOti U 0:2 U ... U cTjv) = Gp{a.i) . (4.11) 

i=l 


Since in general relation (4.9) holds for any set of basic elements a* which are sets of 
internal momenta, one can look at these expressions in different ways, depending on the 
given sets and subsets considered. If we define, for example, the basic subsets a* to 
be given by single momenta g*, and since in that case Gp){qi) = 5{qi), Eq. (4.9) then 
denotes a sum over all possible differing m-tuple cuts for the momenta in the set Pn, 
while the uncut propagators are Feynman propagators. These cuts start from single 
cuts up to the maximal number of cuts given by the term where all the propagators of 
the considered set are cut. Using this notation, the Duality Theorem at one-loop can 
be written in the compact form: 

L^^Hpi,P2, ■ ■ ■ ,Pn) = - [ Gpiai) , (4.12) 

Jii 

where ai as in Eq. (3.2) labels all internal momenta qi. In this way, we directly obtain 
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Pi 


Pi- 



Pn 

Figure 4.1: Momentum configuration of the two-loop A^-point scalar integral. 

the Duality relation between one-loop integrals and single-cut phase-space integrals 
and hence Eq. (4.12) can also be interpreted as the application of the Duality Theorem 
to the given set of momenta ai. It obviously agrees, at one loop, with Eq. (4.1). 

4.2 Duality relation at two—loops 

We now turn to the general two-loop master diagram, as presented in Figure 4.1. Again, 
all external momenta pi are taken as outgoing, and we have pi^j = Pi+ Pi+i + ■ ■ ■ + Pj, 
with momentum conservation pi^N = 0. The label i of the external momenta is defined 
modulo N, i.e., PN+i ^ Pi- In the two-loop case, unlike at the one-loop order, the 
number of external momenta might differ from the number of internal momenta. The 
loop momenta are ii and £ 2 , which flow anti-clockwise and clockwise respectively. The 
momenta of the internal lines are denoted by qi and are explicitly given by 


/ 


4 +pi,i 

Qi = ' 4 + Pi,1-1 

£±+£ 2 +Pi,1-1 




(4.13) 


i G 02 


f G as 


where a^, with k = 1,2,3, are defined as the set of lines, propagators respectively, 
related to the momenta qi, for the following ranges of i: 
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4.2. Duality relation at two-loops 


In the following, we will use for denoting a set of indices or the set of the corresponding 
internal momenta synonymously. Furthermore, we will refer to these lines often simply 
as the “loop lines”. 

We shall now extend the Duality theorem to the two-loop case, by applying Eq. (4.12) 
iteratively. We consider hrst, in the most general form, a set of several loop lines ai to 
um depending on the same integration momentum £j, and find 

/ Gf{oii U 012 U ... U OfAf) = ~ / Gd(q^1 U 0.2 U ... U CTAf) j 

J J 

which states the application of the Duality Theorem, Eq. (4.12), to the set of loop lines 
belonging to the same loop. Eq. (4.15) is the generalization of the Duality Theorem 
found at one-loop to a single loop of a multi-loop diagram. Each subsequent application 
of the Duality Theorem to another loop of the same diagram will introduce an extra 
single cut, and by applying the Duality Theorem as many times as the number of 
loops, a given multi-loop diagram will be opened to a tree-level diagram. The Duality 
Theorem, Eq. (4.15), however, applies only to Feynman propagators, and a subset of the 
loop lines whose propagators are transformed into dual propagators by the application 
of the Duality Theorem to the first loop might also be part of another loop (cf., e.g., 
the “middle” line belonging to 03 in Fig. 4.1). The dual function of the unification 
of several subsets can be expressed in terms of dual and Feynman functions of the 
individual subsets by using Eq. (4.9) (or Eq. (4.10)), and we will use these expressions 
to transform part of the dual propagators into Feynman propagators, in order to apply 
the Duality Theorem to the second loop. Therefore, applying Eq. (4.15) to the loop with 
loop momentum £ 1 , reexpressing the result via Eq. (4.10) in terms of dual and Feynman 
propagators and applying Eq. (4.15) to the second loop with momentum £ 2 , we obtain 
the Duality relation at two loops in the form: 


L^^Hpi,P2, ■ ■ ■ ,Pn) 


(4.16) 



{—GD(ai) 0 ^( 02 ) Goias) + Goiai) Gd{o'2 U as) + Goias) Gd{ 


oi U 02 )] ■ 


This is the dual representation of the two-loop scalar integral as a function of double-cut 
integrals only, since all the terms of the integrand in Eq. (4.16) contain exactly two dual 
functions as defined in Eq. (4.4). The integrand in Eq. (4.16) can then be reinterpreted 
as the sum over tree-level diagrams integrated over a two-body phase-space. 


The integrand in Eq. (4.16), however, contains several dual functions of two different 
loop lines, and hence dual propagators whose dual iO prescription might still depend on 
the integration momenta. This is the case for dual propagators G[){Qi',Qj) where each 
of the momenta qi and Qj belong to different loop lines. If both momenta belong to 
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the same loop line the dependence on the integration momenta in ri{qj — qi) obviously 
cancels, and the complex dual prescription is determined by external momenta only. The 
dual prescription rj(qj — qi) can thus, in some cases, change sign within the integration 
volume, therefore moving up or down the position of the poles in the complex plane. To 
avoid this, we should reexpress the dual representation of the two~loop scalar integral 
in Eq. (4.16) in terms of dual functions of single loop lines. This transformation was 
unnecessary at one-loop because at the lowest order all the internal momenta depend 
on the same integration loop momenta; in other words, there is only a single loop line. 

Inserting Eq. (4.10) in Eq. (4.16) and reordering some terms, we arrive at the following 
representation of the two-loop scalar integral 

L^‘^'>{pi,P2,...,pn) (4.17) 

{Gd{o:i) Gd{oi 2) Gpiofi) + GD{—<y.i) GF{oi2) Goias) + G*{ai) GD{oi2) Goias)} 



il J£2 


where 


G*(ai) = Gi?(a!i) + Gd{o;i) + G£){—ai) 


(4.18) 


In Eq. (4.17), the iO prescription of all the dual propagators depends on external mo¬ 
menta only. Through Eq. (4.18), however, Eq. (4.17) contains also triple cuts, given 
by the contributions with three Goicik)- The triple cuts are such that they split the 
two-loop diagram into two disconnected tree-level diagrams. By definition, however, 
the triple cuts are such that there is no more than one cut per loop line a^- Since there 
is only one loop line at one-loop, it is also clear why we did not generate disconnected 
graphs at this loop order. Eor a higher number of loops, we expect to find at least the 
same number of cuts as the number of loops, and topology-dependent disconnected tree 
diagrams built by cutting up to all the loop lines ak- These results can be generalized 
at three-loops and beyond without any additional effort. The reader is referred to [49] 
for further details. 
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Chapter 5 


Loop—Tree Duality beyond 
Simple Poles 


In the previous chapters, the derivation of the Duality relied on having only simple poles 
in the loop integral, i.e. it was only applicable to Feynman graphs that do not feature 
identical propagators. At one-loop this situation can always be avoided by a convenient 
choice of gauge [48], but for two-loop and higher order corrections this isn’t the case 
anymore. Hence, in this chapter we first want to address the issue of higher order poles. 
We will see that a straightforward application of the Residue Theorem similar to the 
one-loop case is feasible but leads to complex expressions. Therefore a second, more 
elegant, way is shown which involves the use of Integration By Parts (IBP) relations. 
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Figure 5.1; The two-loop A^-point scalar integral with a double pole marked in red. 

5.1 Duality relation for multiple poles 


In the previous chapter we applied the Residue Theorem to one- and two-loop graphs 
that contain only single poles, i.e. no identical propagators. At one-loop this is always 
the case for a suitable choice of gauge [48]. However, at higher loops there exists the 
possibility of identical propagators, i.e. higher order poles [1]. Obviously, we need to 
generalize the Duality Theorem to accommodate for such graphs. The first occurrence 
of higher order poles is at the two-loop level, with the sole double pole generic graph 
shown in Fig. 5.1. The Duality Theorem can be derived for such graphs as well, using 
the Residue Theorem for multiple poles 




1 / 

{k — 1)! \dz^~^ 



(5.1) 


The derivation follows similar steps as with the single pole case and is independent of 
any particular coordinate system. We will derive an expression both in Cartesian and 
light-cone coordinates, to demonstrate this independence. We start with the cartesian 
system. We write the Feynman propagator in a form that makes the poles explicit, i.e. 


Gpiqi) 


1 

(%o - 


(5.2) 


where + mf — iO is the position of the pole. Then, applying the Residue 

Theorem by selecting poles with negative imaginary part, we have 


jjjj <0} Gp{qi ) 


(H 


(+)')3 

iO I 


= - dqo- 


2(<ik¥ 


d+iqf 


(5.3) 
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5.1. Duality relation for multiple poles 


The imaginary component of the new denominator is irrelevant, because it 

is always a positive quantity. We refer the reader to [48] where the calculation for 
the case of simple poles is explained in more detail. Then, we assume the following 
Lorentz-invariant prescription of the residue 

= - J 

where rj^ = (r/o, 0) is a future-like vector, rjQ > 0, in Cartesian coordinates. Contrary to 
the one-loop case, where numerators depending on the loop momentum do not modify 
the Duality prescription, in the two-loop and higher orders cases the derivative in the 
residue calculation introduced by the higher order poles act on every single term in the 
numerator and also on the remaining propagators. Let N{ak) be a function of a set of 
momenta qi, with I £ Ok- Then the residue of a double pole is given by 




YlGD{qi-,qj) 








{H. 


(+)^2 


iO 


--Fiy - ''^{‘^qjo)GD{qi, qj) + 
Qio j^i 


A 

dqo 


N{ak) , 


(5.5) 


which can be written as 


Res{i^q^^<o} I A(^i) f 1 J f 


V 


j¥=i 

1 d 


^ mi 2r]qi dpqi 


N{ak) . 

(5.6) 


In light-cone coordinates we choose our coordinates such that in the plus component 
complex plane the poles with negative imaginary part are located at: 


(+) _ A + mj - iO 

2qi- ’ 


with qi- > 0 . 


(5.7) 


In these light-cone coordinates the Feynman propagator reads: 

GF{q^) = ^ 


2qi-{qi+ - gjj^) 


(5.8) 
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and thus 




(5.9) 


which, at first sight, seems to contradict equation Eq. (5.4). This contradiction can be 
resolved by taking into account the fact that in light cone coordinates, the dual vector 
T] is light-like and therefore rf = 0. Hence equation Eq. (5.4) remains valid. Now, we 
are ready to calculate the residue of a double pole in light cone coordinates: 


Res{img.^<0} <j 9{qi_)G\{qi) | WGpiqj) I N{ak) 

/ 


= 


d{qi-) d 

{2qi-fdq+ 


HGpiqj) I N{ak) 


(ii+=q^iX' 


= I dq+6+{qf - ni^) [ JJ (g*; ) 


^ —Gniq^qj) + 


J_^ 

2riqi drjqi 


N{ak) , 


(5.10) 


where now 77 ^ = ( 77 +, 77 _ = 0, 0_l). Eq. (5.10) has the same functional form as in Eq. (5.6), 
although with a different dual vector rj. Thus we can generalize Eq. (5.6) and Eq. (5.10) 
to an arbitrary coordinate system and combining simple and double poles in a single 
formula we get in Cartesian coordinates: 



U^Fiqj) 


N{ak) = 






^ —GD{qi]qj) + 


_l _^ 

2r]qi drjqi 


+ y^,d{qj)G%{qj]qi) 


Gniq^qk) 




N{ak). 


(5.11) 


Equation (5.11), is the main result of this section. It extends the Duality Theorem 
to integrals with identical propagators or, to put it differently, with double poles in 
the complex plane. For the case of the generic two-loop graph in Fig. 5.1, this re¬ 
sult can be seen as an extension of Eq. (4.10). If we have two groups of momenta, 
ak,a 2 , one of which contains the double propagator, i.e. ak = {qn = (-1 + ^ 2 } and 
02 = {q 2 = h,q 3 = h+Pi,-- ■,qn-i = h +Pi,N-i,q 2 = ^ 2 }, and we denote by a '2 a 
group that contains all the momenta of 02 leading to single poles, namely 
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5.2. Reducing to single poles with IBPs 


“2 = = h,q 3 = h + Pi, ■ ■ ■ ,qn-i = h +Pi,N-i}, then we can write: 


Goiak^ 012) = Hq 2 ) \ GD{q 2 -,qj) 


V 


9( Y1 —GD{q2;qj) 


+ Hqi)Gl{qi-,q2)[\{ Gniq^qj) 


(5.12) 






This result states that for the case of a double pole, one follows the usual procedure of 
cutting every propagator line once, including the double propagator, and transforming 
the rest of the propagators to dual propagators. A similar formula can be derived for the 
case of multiple (triple and higher) poles. The calculation of the residue of a multiple 
pole introduces, however, contributions with powers of dual propagators. In absence of 
a general transformation formula analogous to Eq. (4.9), it is not possible to rewrite 
Eq. (5.11) in terms of dual propagators whose dual +i0 prescription depends on the 
external momenta only. Eor that reason, we will present in the next section a different 
strategy for dealing with higher order poles based on the reduction of the integral using 
Integration By Parts. 


5.2 Reducing to single poles with IBPs 


In this section, we discuss a different approach to the generalization of the Duality 
Theorem to higher order poles. We will use Integration By Parts (IBP) [77, 78] to 
reduce the integrals with multiple poles to ones with simple poles. We emphasize the 
fact the we do not need to reduce the integrals to a particular integral basis. We just 
need to reduce them ’’enough”, so that the higher order poles disappear. 


To give a short introduction to the method and establish our notation, let us consider a 
general m~loop scalar integral in d dimensions, with n denominators Di ,..., Dn raised 
to exponents ai,..., and external momenta pi, ■ ■ ■ ,pn- 



(5.13) 


If we notice that 


d 


tt^ 


= 0 , 


(5.14) 


4 

where s^ = £^,£ 2 ,... ,im, the integrand being a total derivative with respect to the 
loop momenta, we can find relations between scalar integrals with different exponents 
Oj. This will allow us to express integrals with exponents larger than one, in terms of 
simpler ones. In effect, we will be able to write integrals with multiple poles in terms 
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of sums of integrals with simple poles. In the numerator of the integrand of Eq. (5.14) 
we can use ... ,im,Pi, ■ ■ ■ ,p%, to obtain a system of equations that relate the 

various integrals. For simplicity, when referring to an IBP we will use the shorthand 
notation: 

l-t (5,15) 

to denote Eq. (5.14). The differentiation will raise or leave an exponent unchanged, while, 
contractions with the loop and external momenta in the numerator of the integrand, can 
be expressed in terms of the propagators to lower an exponent. Often times though, this 
is not possible, leaving scalar products of momenta, which cannot be expressed in terms 
of denominators. These are called Irreducible Scalar Products (ISP). We will consider 
ISPs as additional denominators, Dij = it ■ pj , with a negative index a*. We use the 
notation: 


F(oia2 ■■■an) = 


>il Ji2 


(5.16) 


to denote a generic two-loop integral with n propagators raised to an arbitrary integer 
power, with Di = Gp^{qi) = qf — mf + iO and qi denotes any combination of external 
and loop momenta. In the following the prescription +z0 for the propagators is un¬ 
derstood. We will use the symbol al*" to denote the raising of the index a, by one i.e. 
l'^E(ai, 02 , • • • an) = F{ai + 1, 02 ,..., an) and the symbol aj” to denote the lowering of 
the index ai by one i.e. 2~F(ai, 02 , • • • o^) = F{ai,a 2 — l,...,an)- A combination of 
the two means that the operators apply at the same time i.e. 1^2~E(ai, 02 , • • • an) = 
F(ai + 1, 02 — 1,..., a„). In the following we will use two automated codes, for the 
reduction, FIRE [79], a MATHEMATICA package for the reduction of integrals and REDUZE 
2 [80] ^ , a package written in C-I-+, using GiNaC [81]. 


5.2.1 The case for two—loop diagrams 

The only generic two-loop scalar graph with Wlegs and a double propagator is shown 
in Fig. 5.1. The simplest case is the two-point function with massless internal lines. 
The denominators are: 

Di=£l , D2 = il , D3 = ii2+p)^ , D^ = {h+i2)^ , D5 = ii^p, 

^ Since the most obvious first approach seems to be to try to express the integrals with multiple poles 
in terms of the same integrals with only single poles, c.f. Eq. (5.18), we used, in addition to the “usual” 
version of REDUZE 2, in some cases a special patch for REDUZE 2 which provides a modification of its 
integral ordering in the final result. This modified version of REDUZE 2 delivered the results for the 
integrals in this desired form stated in the subsequent sections, while we used the normal version of the 
integral ordering for the remaining cases. Note that we also calculated explicitly the relations obtained 
from the modified version, in the easiest cases of the massless two- and three-loop integrals which can 
be built by insertion of the the massless one-loop two-point function, and found agreement. 
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where we have introduced an ISP as an additional denominator For the rest of this 
section the prescription +i0 for the propagators is understood. In our notation, the 
integral we want to reduce is F"( 12110 ) and to this end we use the six total derivatives 

= 1.2. (5.17) 

Applying these IBPs on ^( 0102030405 ) we get a system of recursive equations. Using 
specific values for the exponents o* we can solve this system and obtain U(12110). For 
this particular case, we solve the system explicitly and the reader is referred to the 
Appendix B for details. Finally we arrive at: 

F(12110) = F(llllO), (5.18) 

where s = + i0, a result which contains only single poles and can be treated using the 

Duality Theorem [49] . For the rest of the cases below and in the three-loop case in the 
next section, we have used FIRE and REDUZE 2 to perform the reductions and check our 
results. For three external legs pf = = 0,p| = {pi +^ 2 )^ and massless internal lines, 

we have the denominators: 

Dl=il , 02=^2 , ^3 = (^ 2 +Pl)' , D^ = {i2+Pl+P2f , D5 = ih+i2f , 

De = h • Pi , D 7 = h • P2 , 

where the last two are the ISPs that appear in this case. The integral we want to reduce 
is F(1211100). We use eight IBPs: 

It-w , i,J = l,2. (5.19) 

A similar analysis to the one above, gives: 

F(1211100) = ■^^-^^T(lllllOO) = , (5.20) 

where s = {pi +^ 2 )^ + iO, which, again contains only single poles and can be treated 
with the Duality Theorem. 

The inclusion of masses does not affect the general picture of the reduction. It solely 
introduces numerators in some integrals after the reduction is done. But, as we have 
stressed already, the application of the Duality Theorem is not affected by numerators 
since it only operates on denominators [49]. As an illustrative example, let us consider 
the two-loop graph with two external legs and one massive loop (see Fig. 5.1). For 
^For REDUZE 2 the corresponding propagator is added and used as inpnt instead. 


53 







Chapter 5. Loop-Tree Duality beyond Simple Poles 


the case of the left loop being massive (related to £2); with mass m, the denominators 
involved are 


Di= l\ , D2 = - m'^ , Ds = {£2 + pf' , Di = + £ 2 )^ - , 

with the addition of the irreducible scalar product 


L>5 = £1 • p , 


(5.21) 


needed to perform the reduction. Using the same IBPs of Eq. (5.17), the result of the 
reduction, with FIRE is: 


F(12110) = 


(e - 1) [-es^ + 2m2(9e - 26^ - 3)s + 4m^(-3 + 2e)(-l + 2e)] 


+ 


A 


2 e(2e — l)m^s (4m^ — s) 
(e — 1)^ (2m^ — s) 


2 e(2e — l)m^s (dm^ — s) 

-2^(10110) - 577^^-17^47^(1 “ 1110) 


F(OOllO) 


2 (2e- l)m4s 


, , (e — 1) (4em^ + 2m^ — s) , 

F(OIOIO) - \ ^ V -^F(OlllO) 


+ 


(2e — l)m^s (4m2 — s) 2(2e — l)m^ {Am‘^ — s) 

2 (e — 1) [m? — s) (lOem^ — es — 3m?) 


e(2e — l)m'^s {4m? — sY 


-F(1011 - 1), 


(5.22) 


with s = + iO and 


A = {l-e) [es + 2(3 - 8e)m‘^] s^ + 2(1 - 2e)m'^ [2(3 - 4€)m^ - (6 - 5e)s] . (5.23) 

The reduction generates two integrals with a numerator, namely 

but the double poles have now disappeared. The result with REDUZE 2 reads: 

(e — 1) (4em^ + 2m^ — s) 


F(12110) = - 


-F(OlllO) 


+ 

+ 

+ 


F(1 - 1110) 


2(2e — l)m^ (dm^ — s) 

3(e — 1) (Sem^ — 12em^s + es^ — 4m^ + 4m^s) 

2e(2e — l)m^s (dm^ — s)^ 

- ^ - 2 ^( 10110 ) 

2e(2e — l)m^s (dm^ — s) 

(e — 1) (Se^m^s — 2e^s^ — 16em^ + tSem?s + es^ + 12m'^ — Qm?s) 
2e(2e — l)m^s (dm^ — s)^ 


F(OIOIO) 

(5.24) 
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where A is given by Eq. (5.23). Despite the appearance of different integrals the two 
results are of course equivalent. This is because, the integrals E(OOllO) and E(OIOIO), 
in the result obtained with FIRE, are identical (as can be seen by shifting the loop 
momenta), so the sum of their coefficients gives exactly the coefficient of the result 
obtained with REDUZE 2. The same argument applies for the integrals E(1011 — 1) and 
E(1 — 1110). The appearance of the numerators does not affect the application of the 
Duality Theorem for integrals with single poles as was detailed in [49]. For the case of 
the right loop in Fig. 5.1 being massive (related to ii), we have the denominators: 

, 02=^2 , D3 = i£ 2 +pf , Di = {h+h?-m^ , D 5 = h ■ p . 


Using the IBPs from Eq. (5.17), we get with FIRE: 


, , (32e^m^ + Se^m^s + — 326771“^ — llem^s — + 3m^s) , 

F(12110) = i-iF(lOllO) 


6771*^5^ 

(e-l)(l6 + 4e^s — 20e^777^ — Se'm? — 7es + 3?77^ + 3s) 

6(2e - l)(2e + I)m4s2 

“ ^^^(1011 -1) - (^~ ~ ^) f(oiiio) 


3777^5 2m?s 

(e — 1) (126777^ + es — 3777^) 


6777^5^ 

(e — 1) [6em? + es — 3777^) 
6777 “^s 


F(1 - 1110) 

F(lllOO) , 


F(IOOIO) 


(5.25) 


and with REDUZE 2: 

, , (e — 1) iSem? + es — 2rn?) , 

F(12110) = -iF(l-mO) 


+ 


4777^S^ 

64e^777^ + 16e^777^s + e^s^ — 64e777^ — 22em?s — es^ + 12777'^ + Qm?s 


(e — 1) (l2e777^ + es — 6777 ^) 
6?77^S 


12777^s2 

F(lllOO) 


F(lOllO) 


(e — l)(2e — 3) (l6e^?77^ + 2e^s + 4e777^ + 3es — 2nn? — 2s) 
12(2e- l)(2e + 1 ) 777^52 


F(IOOIO) . 

(5.26) 


For the case of the double pole, two-loop graph, with three external legs and one massive 
loop, we have the denominators: 


Di=l\ , D2 = , D3 = {h+Pif-rri^ , A = ii2+Pi+P2f-m‘^ , D5 = {ii+i2f-m^ , 


Dq=Ii-Pi , D 7 = ii-p2 . 
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Using the IBPs from Eq. (5.19) we get with FIRE: 

(e — 1) {(1 + 4e)s^ — 4em^s(ll — 2e) — 8m^(4e^ — 8e — 1)} 


F(1211100) = 




+ 


2e(2e — l)m®s {Am'^ — s) 
(8e3-12e2 + 4e-l) (e - 1) 
8e(2e — 

(e-1) 


8e(2e — l)m^s — sY 
rF(lOOllOO) 


F(OOOllOO) 


F(OOIOIOO) 


+ 

+ 


2(2e - l)m4 
(e — 1) (2em^ — es — m^) 


[F(OllllOO) - F(lOlllOO)] + 


2(e-l) 


m^s (4m^ — s) 


F(OIOIIOO) 


2e(2e — 

(e — 1)^ {Y>em? — 2es — Qm? + s) 


F(1 - 101100) 


E(OIOOIOO) 


2(2e — l)mPs [Am^ — s) 

2(e — 1) (m^ — s) {—es^ + m^s(6e — 1) + 8m^(2e — 1)} 
e(2e — l)m®s2 {Am‘^ — sY 


F(10011 - 10) 

(5.27) 


where s = (pi + P 2 Y + ^0; and: 


= e(l-e)s^+m2(-l+e)(-l+9e)s''+m''(l-e)(5+6e)s+2m'^(l+2e)(-3+4e) , (5.28) 


„2 , 


while, with REDUZE 2, we get: 


F(1211100) = - 

+ 

+ 


e-1 


rF(OllllOO) 


where 


2(2e- l)m4^ 

F(lOlllOO) 


2(2e- l)m4 
2(e-l) 


-F(OIOIIOO) 


m‘^s {Am‘^ — s) 

3(e — 1) (4em^ — SemYs + es^ + 2m‘^ + rrYs) 
2e(2e — 1 ) 771^5 {Am? — sY 
A 2 


F(1 - 101100) 


+ 


2e(2e — l)m®s (4m^ — s)^ 
^3 

4e(2e — lYm^s {Am? — s) 


F(lOOllOO) 
2 ^( 0100100 ) , 


(5.29) 


A _ i/z?2 6 /?2 4 if\2 22 23 

A 2 = loe m — be m s + 'de m s — e s 

— 4em® + em^s — lOem^s^ + es^ — 6m® + 5m^s + m?s^ , (5.30) 
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and 

A 3 = (e — 1) (l28e^m^ — 64e^m^s + 8e^s^ — 256e^m^ + 112 e^m?s — 12 e^s^ + 192e^m'^ 
— 92e^m^s + 8e^s^ — 40em^ + 26em^s — es^ — 12m^ + Am^s — s^) . (5.31) 

The cases with additional external legs can be treated in a similar manner. It can always 
be rednced to sums of integrals with single propagators at the expense of introducing 
numerators. Although no formal proof exists, in all cases studied so far it has been 
possible to reduce to integrals where only single poles appear. The generality of this 
result seems plausible [82]. 

Our strategy is now clear. For a two-loop calculation, first we reduce all double pole 
graphs using IBPs or any other method. The remaining integrals all contain single poles 
and can be treated using the Duality Theorem at two-loops. The appearance of vector 
or tensor integrals does not spoil this strategy since the Duality Theorem for single poles, 
affects only the denominators of the integrands. 


5.2.2 The case for three—loop diagrams 


For three-loop graphs there exists one topology with a triple propagator and a number 
of topologies with a double propagator. All topologies are shown in Fig. 5.2. The 
arguments for the two-loop case are valid here as well. We first reduce the multiple 
pole integrands by using IBPs until we have integrals with only single poles (possibly 
with numerators) and then we can then apply single-pole Duality Theorem as it was 
described for the three-loop case in Ref. [49]. In the following, we show explicitly the 
reduction of the two-point function for the different topologies and for massless internal 
lines. We use the notation: 


F{aia2 ■■■an) 




(5.32) 


to denote a generic three-loop integral with n propagators raised to an arbitrary integer 
power, with Di = Gp^{qi) = g? — mf -|-i0 or Di equal to any ISP and qi any combination 
of external and loop momenta. We also have s = + iO. The IBPs to be used for the 

reduction are: 

I-'P . >.i = 1 . 2.3 (5.33) 

In the following, we present first the result obtained with REDUZE 2 and then with FIRE. 
For the single triple pole graph (a) in Fig. (5.2), we have the following expressions: 
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Figure 5.2: Master topologies of three-loop scalar integrals with multiple powers of 
internal propagators. Each internal line already dressed with some external leg can be 
dressed with additional external lines. 


(a) 


The denominators used are: 

= £? , T »2 = ^2 , £>3 = , ^4 


(£2-?? , D, 


{il-£2)\ De = i£3-£2f , 


Di — £i ■ p , Ds — £3 ■ p , Dg — £i ■ £3 , 

with the result: 

(5.34) 


For the graphs with doubles poles, (b)-(g), Fig. (5.2), we find: 

(b) The denominators are: 

Di = £1 , D2 = £l , D3 = £l , D^ = {£2-pf , D3 = {£i-£2? , De = (^3-^2+p)' , 
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Dj — ii ■ p , Ds — £i ■ is , Dg — is ■ p , 


with the result: 


F(121211000) 


(1 + 

6(—2 + 3e)(—1 + 3e)(l + 3e)(—3 + 4e)(—1 + 4e) 
e^(l + e)^(—1 + 2e)s^ 


F(IOIOIIOOO) . 


(5.35) 


(c) The denominators are: 


D,=il, D2 = il, Ds=il, Di = (is+pf , Ds = (is-hf , A = {h - hf , 


Dj = ii ■ P , Ds = i 2 ■ p , Dg = ii ■ is 


with the result: 
F(122111000) 


2e(—1 + 4e)(—1 + 3e) 


F(llllllOOO) 


(l + 2e)(l + e)2s2 
2(—2 + 3e)(—1 + 3e)(—3 + 4e)(—1 + 4e) 

e(l + e)2(l + 2e)s4 


F(lOOlllOOO) . 

(5.36) 


(d) The denominators are: 

D,=il, 02=^2, Ds=ih ^4 = (4-p)^ Ds = ii2+i3-h)\ 

D& = h ■ P ^ Dj = h ■ P , Ds = ■ h , Dg = ii ■ is , 

with the result: 


F(112110000) = ^ F(lllllOOOO) = ^ F(llOllOOOO) . (5.37) 


es 


es^ 


(e) The denominators are: 


Di = il , D2 = il , Ds = il , D^ = {is-pf , Ds = {ii+h? , Dg = (£2 + 4 )' , 


D 7 = il ■ p , Ds = il ■ is , Dg = £2 ■ p 


with the result: 


F(112111000) = ^-^^F(llllllOOO) = ^ F(lOOlllOOO) 

(1 ~t~ 2e)s e (1 ~t~ 2e)s 

(5.38) 
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(f) The denominators are; 

D,=il, D 2 = il, = D^ = ih+pf , T »5 = {h+h? , Dq = {li+isf , 
D 7 = {h-i2f , Ds=il-P, Dg=h-P, 

with the result: 

F(121111100) = 

+ 

(g) The denominators are: 

Di=il, 02 = 4 , Ds = el, D^ = (^ 3 -pf , 

D5 = {h —12)^ ■, Dq = (£3 — £2)^ , -D7 = (£3 — ^2 — p4 

= , D^ = {h-pf ( 5 . 40 ) 

with the result: 

The difference between the results of FIRE on the one hand and REDUZE 2 on the other 
is due to the fact that the second is expressed in terms of basis integrals while the first 
is expressed in terms of integrals with single poles of the same type as the multiple pole 
integral (in effect the first result can be further reduced to the second). Since we do not 
seek a particular basis for our reduction, as was stressed earlier, both results are equally 
useful as far as application of the Duality Theorem is concerned. 

5.3 Conclusions 

We have extended the Duality Theorem to two- and three-loop integrals with multiple 
poles. A Lorentz-invariant expression for the residues of double poles has been derived, 
which can be extended straightforwardly to triple and, in general, multiple poles. In the 
absence of a systematic procedure to reexpress dual propagators in terms of Feynman 


---F(lllllllOO) 

(l + e)s 

2(—2 + 3e)(—1 + 3e)(—3 + 4e)(—1 + 4e) 
e2(l + e)(l + 2e)s4 
2(-l + 2e)2(-l+ 4e) -|n-|QQ'> 

———^ F(iomoioo), 


[F(OOllllOOO) +F(100101100)] 


(5.39) 
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propagators (cf. Eq. (4.9) for the case of simple poles) we have explored an alternative 
approach. We use IBP identities to reduce the integrals with identical propagators to 
ones with only single poles. Therefore, the essential features of the Loop-Tree Duality 
now remain intact. We reiterate that our goal is not to reduce everything to some 
set of master integrals. Rather, we reduce the integrals until there are no multiple 
poles left. Then, we can use the Duality Theorem in its original form for single pole 
propagators, to rewrite them as integrals of a tree-level object over a modified phase- 
space. The appearance of additional tensor integrals, due to the reduction, does not 
affect our procedure, since applying the Duality Theorem in its single-pole version, only 
cuts propagators, leaving the numerators of the integrals unaffected. 
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Chapter 6 


On the Cancellation of 
Singnlarities 


In the previous chapters the Loop-Tree Duality was extended to cover Feynman dia¬ 
grams which involve multiple loops as well as higher order poles [1, 49]. 

Keeping in mind that the final aim of the Loop-Tree Duality is to treat virtual and 
real corrections at the same time. In this chapter we analyse the singular behaviour 
of one-loop integrals and scattering amplitudes in the framework of the Loop-Tree Du¬ 
ality method. We begin with a discussion of the cancellation of singularities among 
dual contributions at the integrand level. After that, we present a phase space mapping 
between virtual corrections in the dual representation and the real corrections for the 
local cancellation of infrared divergencies. 
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6.1 The singular behaviour of the loop integrand 


We consider a general one-loop V-leg scalar integral 


L^^\pi,P2, • • ■, 



d'^i 

-r • 


( 6 . 1 ) 




where 



( 6 . 2 ) 


are Feynman propagators that depend on the loop momentum i, which flows anti¬ 
clockwise, and the four-momenta of the external legs pi, i G ai = {1, 2,... iV}, which are 
taken as outgoing and are ordered clockwise. We use dimensional regularization with d 
the number of space-time dimensions. The momenta of the internal lines = {qifi,c[i), 
where Qi^o is the energy (time component) and qj are the spacial components, are defined 
as Qi = i + ki with ki = pi + ... + pi, and /cat = 0 by momentum conservation. We also 
define kji = Qj — qi. 

The loop integrand becomes singular in regions of the loop momentum space in which 
subsets of internal lines go on-shell, although the existence of singular points of the 
integrand is not enough to ensure the emergence in the loop integral of divergences 
in the dimensional regularization parameter. Nevertheless, numerical integration over 
integrable singularities still requires a contour deformation [83-90], namely, to promote 
the loop momentum to the complex plane in order to smoothen the loop matrix elements 
in the singular regions of the loop integrand. Hence, the relevance to identify accurately 
all the integrand singularities. 

In Cartesian coordinates, the Feynman propagator in Eq. (6.2) becomes singular at 
hyperboloids with origin in —ki, where the minimal distance between each hyperboloid 
and its origin is determined by the internal mass m*. This is illustrated in Fig. 6.1, 
where for simplicity we work in d = 2 space-time dimensions. Figure 6.1 (left) shows 
a typical kinematical situation where two momenta, ki and ^ 2 , are separated by a 
time-like distance, fell > 0, and a third momentum, k^, is space-like separated with 
respect to the other two, fe|i < 0 and k‘^2 < 0- The on-shell forward hyperboloids 
(?i,o > 0) are represented in Fig. 6.1 by solid lines, and the backward hyperboloids 
(9i,o < 0) by dashed lines. For the discussion that will follow it is important to stress that 
Feynman propagators become positive inside the respective hyperboloid and negative 
outside. Two or more Feynman propagators become simultaneously singular where their 
respective hyperboloids intersect. In most cases, these singularities, due to normal or 
anomalous thresholds [91, 92] of intermediate states, are integrable. However, if two 
massless propagators are separated by a light-like distance, fe|j = 0, then the overlap 
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of the respective light-cones is tangential, as illustrated in Fig. 6.1 (right), and leads 
to non-integrable collinear singularities. In addition, massless propagators can generate 
soft singularities at = 0. 



Figure 6.1: On-shell hyperboloids for three arbitrary propagators in Cartesian coor¬ 
dinates in the {(.q/z) space (left). Kinematical configuration with infrared singularities 
(right). In the latter case, the on-shell hyperboloids degenerate to light-cones. 


As we have seen in Chapter 3, the dual representation of the scalar one-loop integral in 
Eq. (3.1) is the sum of N dual integrals [48, 49]: 


L^^\pi,P2,---,Pn) = f d{qi) GD{qi\qj) , (6.3) 

i&ai ^ iScfi 

j¥=i 


where 


GD{qi;qj) 


1 

q‘j - m'j - iO r] ■ kji 


(6.4) 


are the so-called dual propagators, as dehned in Ref. [48], with rj a future-like vector, rj^ > 
0, with positive definite energy r/o > 0. The delta function 5{qi) = 27ri 0{qifl) 6{qf — m?) 
sets the internal lines on-shell by selecting the pole of the propagators with positive 
energy qi^Q and negative imaginary part. In the following we take t]^ = (1,0), and thus 
—iOr] ■ kji = —iOkjifi. This is equivalent to performing the loop integration along the 
on-shell forward hyperboloids. Let us mention that in the light-cone coordinates (£+, 
1 _l), where £± = {io ±id-i)/y/2, Feynman propagators vanish at hyperboloids in the 
plane {£+,£-) which are similar to those depicted in Fig. 6.1 but rotated by 45 degrees. 
Consequently, by selecting the forward hyperboloids the integration limits of either £+ 
or £- are restricted and the restrictions are different for each dual contribution. For 
this reason, although Eq. (6.3) is valid for any system of coordinates, we will stick for 
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the rest of the thesis to Cartesian coordinates where all the dual contributions share the 
same integration limits for the loop three-momentum. 


A crucial point of our discussion is the observation that dual propagators can be rewritten 
as 


5{qi) Goiquqj) = i27r 






( + )^ 


(6.5) 


where 




( 6 . 6 ) 


is the loop energy measured along the on-shell hyperboloid with origin at —ki. By 
definition we have Re(gj^Q^) > 0. The factor can become singular for mi = 0, 

but the integral ffh(qi^o — Qjo^)/Qi~o^ is still convergent by two powers in the infrared. 
Soft singularities require two dual propagators, where each of the two dual propagators 
contributes with one power in the infrared. From Eq. (6.5) it is obvious that dual prop¬ 
agators become singular, G~j^{qi] qj) = 0, if one of the following conditions is fulfilled: 


+ ^jifi — 0 , (6.7) 

= 0 • ( 6 - 8 ) 


The first condition, Eq. (6.7), is satisfied if the forward hyperboloid of —ki intersects 
with the backward hyperboloid of —kj. The second condition, Eq. (6.8), is true when 
the two forward hyperboloids intersect each other. 

In the massless case, Eq. (6.7) and Eq. (6.8) are the equations of conic sections in the 
loop three-momentum space; and are the distance to the foci located at —kj 
and —kj, respectively, and the distance between the foci is i/kT. If internal masses 
are non-vanishing, Eq. (6.6) can be reinterpreted as the distance associated to a four¬ 
dimensional space with one “massive” dimension and the foci now located at (— kj, — mj) 
and (—kj, —mj), respectively. Then, the singularity arises at the intersection of the conic 
sections given by Eq. (6.7) or Eq. (6.8) in this generalized space with the zero mass plane. 
This picture is useful to identify the singular regions of the loop integrand in the loop 
three-momentum space. 

The solution to Eq. (6.7) is an ellipsoid and clearly requires kji^ < 0. Moreover, since it 
is the result of the intersection of a forward with a backward hyperboloid the distance 
between the two propagators has to be future-like, kjj^ > 0. Actually, internal masses 
restrict this condition. Bearing in mind the image of the conic sections in the generalized 
massive space so we can deduce intuitively that Eq. (6.7) has solution for 


A:|j — {mj + mi)'^ > 0 , kjifl < 0 , forward with backward hyperboloids . (6.9) 
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The second equation, Eq. (6.8), leads to a hyperboloid in the generalized space, and 
there are solutions for kjifi either positive or negative, namely when either of the two 
momenta are set on-shell. However, by interpreting the result in the generalized space 
it is clear that the intersection with the zero mass plane does not always exist, and 
if it exists, it can be either an ellipsoid or a hyperboloid in the loop three-momentum 
space. Here, the distance between the momenta of the propagators has to be space-like, 
although also time-like configurations can fulfill Eq. (6.8) as far as the time-like distance 
is small or close to light-like. The following condition is necessary: 

kj^ — {mj — mi)^ < 0 , two forward hyperboloids . (6.10) 

In any other configuration, the singularity appears for loop three-momenta with imagi¬ 
nary components. 


6.2 Cancellation of singularities among dual integrands 

In this section we prove one of the main properties of the Loop-Tree Duality method, 
namely the partial cancellation of singularities among different dual integrands. This 
represents a significant advantage with respect to the integration of regular loop integrals 
in the d-dimensional space, where one single integrand cannot obviously lead to such 
cancellation. 

Let us consider first two Eeynman propagators separated by a space-like distance, < 0 
(or more generally fulfilling Eq. (6.10)). In the corresponding dual representation one 
of these propagators is set on-shell and the other becomes dual, and the integration 
occurs along the respective on-shell forward hyperboloids. See again Eig. 6.1 (left) for a 
graphical representation of this set-up. There, the two forward hyperboloids of —ki and 
—ks intersect at a single point. Integrating over along the forward hyperboloid of —ki 
we find that the dual propagator Goigiigs), which is negative below the intersection 
point where the integrand becomes singular, changes sign above this point as we move 
from outside to inside the on-shell hyperboloid of —k^. The opposite occurs if we set 
qs on-shell; GD{q 3 ',qi) is positive below the intersection point, and negative above. The 
change of sign leads to the cancellation of the common singularity. Notice that also the 
dual iO prescription changes sign. In order to prove analytically this cancellation, we 
define x = -|- kjifl. In the limit x —>■ 0: 

lim (5{qi) Gnigi] qj) + {i ^ J)j = 

\x xj 
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and thus the leading singular behaviour cancels among the two dual contributions. The 
cancellation of these singularities is not altered by the presence of other non-vanishing 
dual propagators (neither by numerators) because 

lim GD{qj-,qk) = hm --rr^- - ---- 7 -wr = 1 ™ Gniqu qk) , ( 6 . 12 ) 

where we have used the identity /cfcj.o = kkifl~kjifl. If instead, the separation is time-like 
(in the sense of Eq. (6.9)), we define x = -|- kjifl, and find 

(6{qi) Goiqi] qj) + {i ^ j)) = -^{-kjifl) ^ 5{qi) + {i ^ j) + 0(x°) . (6.13) 

In this case the singularity of the integrand remains because of the Heaviside step func¬ 
tion. 


We should consider also the case in which more than two propagators become simulta¬ 
neously singular. To analyse the intersection of three forward hyperboloids, we define 

^^ - qj'fl + kjifi , Xy = - qi^d + kki,o . (6.14) 

As before, we use the identity kkjfl = kkifl — kjifl, and thns gj — q^dj + kkjfl = \ {y — x). 
In the limit in which the three propagators become simultaneously singular: 


lim 

A-:-0 


:) G'd(%; qj) Gniqi; qk) + perm.) = 


1 

W 


( — + ^ 

\xy x{x-y) 


+ 


1 


y{y-x) 




^5(g,) + 0(A-i) , 


(6.15) 


and again the leading singular behaviour cancels in the sum. Although not shown 
for simplicity in Eq. (6.15), also the 0(A“^) terms cancel in the sum, thus rendering 
the integrand finite in the limit A —)• 0. Eor three propagators there are also more 
possibilities: two forward hyperboloids might intersect simnltaneously with a backward 
hyperboloid, or two backward hyperboloids might intersect with a forward hyperboloid. 
In the former case, we define Ax = -|- q^^d + and Xy = -|-g[^^ -|- k^jfi, with 

kkifl < 0 and kkjfi < 0) and hence g|^^ — g^Q^ -|- kjifi = X{x — y). In the A —)• 0 limit 


lim 

A-:>0 


) Goiqv, qj) Goiqi] qk) + perm. = 


e{ kk^,o)ei kkjfl) ^2 [xiy-x)^ yix-y) 


2 ( 7 ^+) 

^yfl ^%fl 


5{qi) + 0{X-^) . (6.16) 


Notice that the singularity in l/(x — y) cancels in Eq. (6.16) (also at C>(A ^)). In the 
latter case, we set as before Xx = gj^^ -|- q^j^d + kkifl, and define Az = g^-^^ -|--|- kjifl, 
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then 


lim (^S{qi)GDiqi;qj)GD{qi;qk) + perm.^ = -6{-kkifl) 

(^) 2 ^) 

\ / zg^.Q zg^ Q 

Similarly, it is straightforward to prove that four forward hyperboloids do not lead to 
any common singularity and more generally that the remaining multiple singularities 
are only driven by propagators that are time-like connected and less energetic than the 
propagator which is set on-shell. 

Thus, we conclude that singularities of space-like separated propagators occurring in 
the intersection of on-shell forward hyperboloids, are absent in the dual representation 
of the loop integrand. The cancellation of these singularities at the integrand level 
already represents a big advantage of the Loop-Tree Duality with respect to the direct 
integration in the four-dimensional loop space; it makes unnecessary the use of contour 
deformation to deal numerically with the integrable singularities of these configurations. 
This conclusion is also valid for loop scattering amplitudes. Moreover, this property can 
be extended in a straightforward manner to prove the partial cancellation of infrared 
singularities. 

Collinear singularities occur when two massless propagators are separated by a light¬ 
like distance, k^- = 0. In that case, the corresponding light-cones overlap tangentially 
along an inhnite interval. Assuming ki^ > kj^, however, the collinear singularity for 
Iq > —kjfi appears at the intersection of the two forward light-cones, with the forward 
light-cone of —kj located inside the forward light-cone of —ki, or equivalently, with 
the forward light-cone of —ki located outside the forward light-cone of —kj, Thus, the 
singular behaviour of the two dual components cancel against each other, following the 
same qualitative arguments given before. For —kifi < io < -kj^, instead, it is the 
forward light-cone of —ki that intersects tangentially with the backward light-cone of 
—kj according to Eq. (6.7). The collinear divergences survive in this energy strip, which 
indeed also limits the range of the loop three-momentum where infrared divergences 
can arise. If there are several reference momenta separated by light-like distances the 
infrared strip is limited by the minimal and maximal energies of the external momenta. 
The soft singularity of the integrand at = 0 leads to soft divergences only if two 
other propagators, each one contributing with one power in the infrared, are light-like 
separated from —ki. In Fig. 6.1 (right) this condition is fulfilled only at = 0, but 
not at q!f\^ = 0 neither at q^~\^ = 0. 

^ Including light-like and time-like conhgurations such that Eq. (6.10) is fulhlled. 
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In summary, both threshold and infrared singularities are constrained in the dual rep¬ 
resentation of the loop integrand to a finite region where the loop three-momentum is 
of the order of the external momenta. Singularities outside this region, occurring in the 
intersection of on-shell forward hyperboloids or light-cones, cancel in the sum of all the 
dual contributions. 


6.3 Cancellation of infrared singularities with real correc¬ 
tions 


Having constrained the loop singularities to a finite region of the loop momentum space, 
we discuss now how to map this region into the finite-size phase-space of the real cor¬ 
rections for the cancellation of the remaining infrared singularities. The use of collinear 
factorization and splitting matrices, encoding the collinear singular behaviour of scat¬ 
tering amplitudes as introduced in Ref. [93, 94], is suitable for this discussion. 



Figure 6.2; Factorization of the dual one-loop and tree-level squared amplitudes in 
the collinear limit. The dashed line represents the momentum conservation cut. 


We consider the interference of the one-loop scattering amplitude with the corre¬ 
sponding Wparton tree-level scattering amplitude which is integrated with the 

appropriate phase-space factor 


j d4>7v(pi;p2, • • • ,PAr) 


N 


N 






i=l 


(6.18) 


where we assume that only the external momentum pi is incoming (pi,o < 0). Then, we 
select the corresponding dual contribution with the internal massless line Qi on-shell 



X 


2Re j d^N{pi-,P 2 ,-■ ■ ,Pn) 

{M^^\pi,.. .,Pn)\-Q hQi^Pi+i^ ■■■)) > 


(6.19) 


where the loop energy in Eq. (6.19) is restricted by the energy of the adjacent external 
massless particle pip to select the infrared sector, according to the discussion of the 
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previous sections. We also consider the N + 1-parton tree-level scattering amplitude 




( 0 ), 2 r / 


; Pir 


Pi+Pr,...)) 


( 6 . 20 ) 


where an extra particle is radiated from parton i, with = p[ + p'^, and the com¬ 
plementary scattering amplitude that contains all the tree-level contributions 

with the exception of those already included in ■ The corresponding interference, 

integrated over the phase-space of the final-state particles, is 

/®=2Re J (6-21) 


For the simplicity of the presentation, we do not consider explicitly in this section the 
square of , which is related with a self-energy insertion in an external leg and 

whose infrared divergences are removed by wave-function renormalization [48] . The final- 
state external momenta of the loop and tree amplitudes in Eq. (6.19) and Eq. (6.21), 
although labelled with the same indices, are constrained by different phase-space mo¬ 
mentum conservation delta functions. A mapping between the primed (real amplitudes) 
and unprimed (virtual amplitudes) momenta is necessary to show the cancellation of 
collinear divergences. 

In the limit where p* and q* become collinear the dual one-loop matrix element 
in Eq. (6.19) factorizes as 


\-M^N+ 2 i- ■ -^Ph-Qi^Qi, •••))= Sp^^\pi, -qp, -Qi-i) \Mn+i{- ■ ■,-Qi-i,qi, ■ ■ , 

( 6 . 22 ) 

where the reduced matrix element is obtained by replacing the two collinear 

partons of ^ single parent parton with light-like momentum 


~a u 

C-1 = Qi-i 


2nqi-i 


(6.23) 


with a light-like vector, = 0. Similarly, in the limit where p( and p), become 
collinear the tree-level matrix element factorizes as 


{M%\"{{pi,P 2 , ... ,p'n+i)\ = ■ ■ ,p'i-i,Pir,Pi+i, ■ ■ 01 Sp^'^'>'^{Pi,Pr;Pi, 


/T^(0) 


where and 


^11 111 
Pir = Pir 




2 np'- 


(6.24) 


(6.25) 


is the light-like momentum of the parent parton. A graphical representation of the 
collinear limit of both virtual and real corrections is illustrated in Fig. 6.2. This graph 
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suggests that in the collinear limit the mapping between the four-momenta of the virtual 
and real matrix elements should be such that pi = pj = p'j{j 7 ^ i), —qt-i = p' and 
Qi = p'^ in the collinear limit. Notice that p'^ is restricted by momentum conservation but 
qi is not. However, the relevant infrared region is bound by < pi^ in Eq. (6.19). This 
restriction allows to map qi to p(,. The mapping, nevertheless, is not as obvious as can be 
induced from Fig. 6.2 as the propagators that become singular in the collinear limit in 
the virtual and real matrix elements are different. Reconsidering p( as the parent parton 
momentum of the collinear splitting, we find the following relation between splitting 
matrices entering the real matrix elements 




(6.26) 


where (p(^ ~ Pr)‘^lsir — show now that the factor —np(/np'^ is compen¬ 

sated by the phase-space. By introducing the following identity in the phase-space of 
the real corrections 

1 = J § 0 ) (p'^ - p' - p(.) , (6.27) 

and performing the integration over the three-momentum p( and the energy component 
of p'^, the real phase-space becomes 


d^>Ar+i(pi;P2, 


d$Ar(pi;...,p'^,...) / 5(p(. 

■'p'r 


El 

E[ 


(6.28) 


where the factor {np[/npD{El/E'f) equals unity in the collinear limit. Inserting Eq. (6.22) 
in Eq. (6.19), and Eq. (6.24), Eq. (6.26) and Eq. (6.28) in Eq. (6.21) the loop and tree 
contributions show to have a very similar structure with opposite sign and match each 
other at the integrand level in the collinear limit. Correspondingly, soft singularities at 
p(, —>■ 0 can be treated consistently as the endpoint limit of the collinear mapping. 


6.4 Conclusions and outlook 

The Loop-Tree Duality method exhibits attractive theoretical aspects and nice proper¬ 
ties which are manifested by a direct physical interpretation of the singular behaviour 
of the loop integrand. Integrand singularities occurring in the intersection of on-shell 
forward hyperboloids or light-cones cancel among dual contributions. The remaining 
singularities, excluding UV divergences, are found in the intersection of forward with 
backward on-shell hyperboloids or light-cones and are produced by dual propagators 
that are light-like or time-like separated and less energetic than the internal propagator 
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that is set on-shell. Therefore, these singularities can be interpreted in terms of causal¬ 
ity and are restricted to a finite region of the loop three-momentum space, which is of 
the size of the external momenta. As a result, a local mapping at the integrand level 
is possible between one-loop and tree-level matrix elements to cancel soft and collinear 
divergences. One can anticipate that a similar analysis at higher orders of the Loop-Tree 
Duality relation is expected to provide equally interesting results. We leave this analysis 
for future work. 
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Chapter 7 


Multi-leg Scalar Integrals 


After all these theoretical foundations haven been established in the previous chapters, 
let us turn to an actual numerical implementation. We have written a computer program 
in C++ to give a practical proof-of-concept of the theoretical framework, and to be used 
for the calculation of one-loop diagrams. The program at its current state will serve as 
the basis of any further developments. The specifics of the program are given in Section 
7.7.1. For the moment let us add here that from the point of view of a dual integral, 
calculating a phase-space point in the absence of threshold singularities is comparatively 
easy. Things become more challenging for points with complex results, because these 
feature singularities that do not cancel among dual contributions which means they have 
to be dealt with by contour deformation. Since the energy component is already fixed 
due the prior application of the Residue Theorem, a dual integral and consequently the 
contour deformation lives only in three spatial dimensions. 
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7.1 Preparation of the Numerical Implementation 


The goal is to calculate a one-loop scalar integral in its dual representation: 


i&ai ^ j&ai 


(7.1) 


Therefore, the program relies on a set of fundamental equations. These shall be presented 
here. Since the program shall work for generic masses, we need, as a first ingredient, 
the massive dual propagator 


GD{qi\qj) 


1 

q‘j - mj - iO rj kji 


(7.2) 


However, in practice we will make use of its rewritten form as in Eq. (6.5) 


o{qi)GD{qi]qj) = ^ 27 ^--^ 


where 


ql^^ = -\-mj -iO and kji = kj - ki . 


(7.3) 


(7.4) 


This means actually, that the integration to perform is a three-dimensional one, because 
q* = £ -|- kj where I is the loop three-momentum and kj = X]}=i Pi 1^^® 
external three-momenta up to the i-th. The zero-component has already been integrated 
out as explained in Chapter 3 and as a result the gj^o-component is fixed according to 
Eq. (7.4). This is indicated by the presence of the delta function 5{qi). 

The next ingredient is given by 


k'^i — {mj -|- ruj)^ > 0, kji^ < 0 ellipsoid singularity (7-6) 

kji — {mj — mi)^ < 0 hyperboloid singularity (7-6) 


These two equations are used to identify the type of singularity that we are dealing with. 
In Section 7.3 we will see how they come into play. 


7.2 The Mapping 

We use Cuhre and VEGAS from the Cuba library [71] as numerical integrators. Cuhre 
is a deterministic integrator that uses cubature rules and performs well as long as the 
dimensionality of the integral is not too high. As stated before, dual integrals are only 
three-dimensional making Cuhre an excellent choice for the task. VEGAS on the other 
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hand, is a Monte Carlo-integrator that uses importance sampling for variance reduction. 
While it is much slower, it is useful to crosscheck results and most probably will be used 
in later stages when real and virtual corrections are combined. 

Both integrators have in common that they assume the integration region to be the unit 
cube. Hence, in order to perform an integral over the entire phase-space, we must find 
a mapping (— 00 , 00 )^ —)• [0,1]^. 

The first and easiest option that comes to mind is simply rescaling every dimension. 
One possible mapping to do so is 

= tan ^TT ^ , f = l,2,3 (7.7) 

When Xi —)• 0, the argument of the tangent goes to —7r/2 and the tangent of it to — 00 , 
when Xj —)■ 1, the argument goes to -|-7r/2 and the tangent of it to -|-oo. 

However, with regard to program stability it is advantageous to choose a mapping that 
relies on the use of spherical coordinates. In fact, it is a two-step process: In the first 
step the Cartesian coordinates lx, ^y, (-z running from —00 to -|-oo get mapped to r, cos(0) 
and (/>, with r being the radius, (j) azimuth and 6 the polar angle, respectively. 

^x = r cos((/))\/W-'cos2^ 

^y = rsm.{4>)\Jl — cos^(0) 

£z = rcos(0) (7.8) 

In the second step these get mapped to x,y,z G [0,1]. 

X 

r = - 

X — 1 

cos{0) = 1 — 2y 

(/> = 2Trz (7.9) 

Despite seeming more complex at first glance, it is actually slightly faster. Of course, 
every change of variables demands the inclusion of the corresponding Jacobian. 
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7.3 Singular behaviour of dual contributions 

Dual integrals feature certain types of singularities. This has already been thoroughly 
discussed in Section 6.2, thus, in this section, we want to quickly recapitulate the parts 
that are going to be relevant to a numerical implementation. For generic masses, the 
loop integrand becomes singular at on-shell hyperboloids with -|- m? — iO 

(forward-hyperboloids, solid lines of Fig. 6.1) and mf — iO (backward- 

hyperboloids, dashed lines of Fig. 6.1). The origins of the corresponding on-shell hyper¬ 
boloids are at —kj. There are two main types of singularities to distinguish: 

• Forward-foward intersection. In Fig. 6.1, these are the singularities that corre¬ 
spond to the intersection of two solid lines. They cancel among dual contribu¬ 
tions. In short, the reason is the following (for more details, see Section 6.2): 
Propagators are positive inside and negative outside. When integrating along the 
forward-hyperboloids, every singularity is passed twice. One time going from the 
inside to the outside (or vice versa) and the second time from the outside to the 
inside (or vice versa). The crucial point is that therefore, the contributions coming 
from the two integrations have opposite sign and thus cancel out. 

The situation for this type of singularity is drawn in Fig. 7.1, once in loop three- 
momentum space (7.1b) and for easier legibility in two dimensions (7.1a). The 
surface (curve) represents all points for which the dual propagator in considera¬ 
tion has a pole. The dots indicate the positions of the foci (= —k*) as introduced 
in Chapter 6. 

• Forward-backward intersection. These singularities originate from the intersection 
of a solid with a dashed line in Fig. 6.1. They remain and require to be dealt with 
by contour deformation. In Section 7.5 we will see how this is done. 

In Fig. 7.2, the shapes of the surfaces at which propagators with such a singularity 
become infinite are illustrated. Again we show the actual loop three-momentum 
space plot (Figure 7.2b) alongside a simpler two-dimensional version (Fig. 7.2a). 

Because Fig. 6.1 is only a two-dimensional plot these singularities show up as mere 
points. In reality, i.e. in 1 -|- 3 dimensions (with the expression 1 -|- n, n G N, dimensions, 
we refer to a Minkowski space with 1 time and n spacial dimensions), these singularities 
have the shape of hyperboloids or ellipsoids, respectively, as has been shown in Fig. 7.1b 
and 7.2b. Therefore, we will call them hyperboloid or ellipsoid singularities from now 
on. 
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4 

(a) Two-dimensional hyperbola 




(b) Three-dimensional hyperboloid 


Figure 7.1: Example of a hyperboloid singularity. In 1+2 dimensions the hyperboloid 
degenerates to a hyperbola, as displayed on the left. 



h 

(a) Two-dimensional ellipse 


(b) Three-dimensional ellipsoid 



Figure 7.2: Example of an ellipsoid singularity. In 1-1-2 dimensions the ellipsoid 
degenerates to an ellipse, as displayed on the left. 


7.4 Interplay of dual contributions 


Loop-Tree Duality transforms one-loop integrals into so-called dual integrals: 

r _ r ^ 

L^^\pi,P2, ... ,pn) =n^“ X] y n ^Diqiqj) (7.10) 




JSOI 
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with a set of internal momenta belonging to the same loop. 

Symbolically speaking, we could express Duality with the following scheme: 


Gf-Gf---Gf 


Duality 


<5 

Gd 

Gd 

••• Gd 

Gd 

5 

Gd 

••• Gd 

Gd 

Gd 

5 

••• Gd 

Gd 

Gd 

Gd 

... s 


Figure 7.3: Symbolical representation of the Loop-Tree Duality. Gp = Feynman 
propagator, Go = dual propagator, 5 = 5{qi). 


Each line in the matrix-like structure on the right side of the arrow of Fig. 7.3 represents 
a dual contribution whereas the columns give the number of the corresponding external 
leg. This scheme can now be used to indicate the positions of the different singularities in 
a dual integral. Consider the following small example, in which zero means no singularity, 
H stands for hyperboloid singularity and E for ellipsoid singularity: 


0 

H 

0 

H 

H 

0 

H 

E 

E 

H 

0 

E 

H 

0 

0 

0 


The correct interpretation would be: There is a one-loop box integral, hence we have 
four dual contributions. On the main diagonal we have only zeros because the dual 
delta functions cannot produce hyperboloid or ellipsoid singularities. The first dual 
contribution has two hyperboloid singularities, namely at positions two and four, the 
second dual contribution has hyperboloid singularities at positions one and three and 
an ellipsoid singularity at position four. The third dual contribution has two ellipsoid 
singularities at positions one and four and one hyperboloid singularity at position two 
and the fourth dual contribution has only one hyperboloid singularity at position one. 
This way of denoting singularities will become useful momentarily. 

There is one further observation to make: Apparently the hyperboloid singularities are 
distributed symmetrically around the main diagonal. This is not by accident. Going 
back to Section 6.1 we have Eq. (6.8) which is the defining equation for hyperboloid 
singularities. Due to its symmetry under the exchange of i {i counts dual contributions) 
and j {j counts leg positions) the hyperboloid singularities always appear in pairs and 
are distributed symmetrically around the main diagonal of (7.11). Inspecting Eq. (6.7), 
which is the defining equation for ellipsoid singularities, we see that this equation is not 
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symmetric under the exchange of indices. Thus for every ellipsoid singularity in 7.11 we 
have a zero as its counterpart. 

In Section 7.3, and with more detail in Section 6.2, we established that hyperboloid 
singularities cancel among dual contributions and therefore do not need to be treated 
via contour deformation. Yet they impact the way we have to deform. In order to 
preserve the cancellation of hyperboloid singularities, dual contributions featuring the 
same hyperboloid singularity (pair) must receive the same deformation. To further 
illustrate this point, let us look at the following pentagon example: 


0 

H 

0 

0 



Contributions are coupled: 

H 

0 

H 

E 



Every contribution receives all deformations 

E 

H 

0 

E 

0 

1 

that occur within the coupling. 

E 

E 

E 

0 

E -4 

Deform with ellipsoids that itself contains. 

0 

0 

0 

0 

0 ^ 

No deformation needed here. 


Figure 7.4: Pentagon with dual contributions coupled by hyperboloid singularities. 


In Fig. 7.4, contributions one, two and three are coupled via their common hyperboloid 
singularities. Thus, they need to receive the very same deformation that accounts for all 
ellipsoid singularities occurring within those contributions. These are found at position 
four of the second contribution and positions one and four of the third contribution. The 
fourth dual contribution is not coupled to any other contribution and can be deformed 
as standalone. The fifth contribution does not require any treatment. 

As a general strategy, one organizes the dual contributions into groups. A group is a 
set of pairwisely coupled contributions. Each of the groups is deformed independently 
from the others. Within a group every contribution receives the same deformation that 
accounts for all the ellipsoids of the group. 

Turning back to the example of Fig. 7.4, we would have three groups: the first group 
involves contributions one to three, the second group constitutes of contribution four 
and the third group constitutes of contribution five. 


7.5 Contour deformation 

The ellipsoid singularities (forward-backward type) lie on the real axis. To avoid them, 
we need to deform the integration path into the imaginary space. 
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7.5.1 A one-dimensional example 


To motivate the basic concept of contour deformation, let us have a look at the following 
simple example. The function 


/( 4 ) 


1 

+ iO 


(7.12) 


has poles at Ix± = -ii{E — it)). Simply integrating along the real axis would lead to 
infinities. Therefore, an integration path that goes around the singular points is needed. 
One possible way to achieve this would be: 


— 4 + i^(-x 


exp 


4-^2 

2E'^ 


(7.13) 


The parameter A serves to scale the deformation along the imaginary axis. At the 
position of the pole, the exponent becomes 0 and thus the exponential function hits its 
maximum, which is 1. Far away from the poles the exponent is a large negative number, 
hence exponentiating it suppresses the deformation. 


Shape of the contour deformation 



Re(lx') 

Figure 7.5: Contour deformation as in Eq. (7.13) for A = 5 and E = 2. 


7.5.2 Deformation in three dimensions 

Every valid deformation must satisfy a certain set of requirements [87] : 
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1. The deformation has to respect the iO-prescription of the propagator: 

In general, a (d-dimensional) contour deformation has the form: 

= £ + (7.14) 

where k usually is a function of the loop momentum £. In our case, we want to per¬ 
form the integration over a product of dual propagators. Plugging the deformation 
of Eq. (7.14) into the on-shell energy (Eq. (7.4)), we obtain 

2iK • q* -h qf -h m? - iO . (7.15) 

The Eeynman prescription —iO tells us in which direction to deform when coming 
close to a singularity. Hence, any valid deformation must match this prescription. 
Consequently, we need to have 


k: • q* < 0 . 


(7.16) 


2. The deformation should vanish at infinity: 

We are looking for a deformation that does not change the actual value of the 
integral. Therefore, we do not want |ft:| to grow for \i\ —?■ oo. An easy way to 
satisfy this condition is to choose k, such that |«:| —)• 0 as \i\ —)• oo.^ 


With these conditions in mind, we construct the deformation in the following way: 

As explained in Section 7.4, we first organize the dual contributions into groups. Eor 
every ellipsoid singularity of the group we include a factor: 




^ ) 


(7.17) 


with qi = i + 'k.i and i the loop three-momentum. It is made up of two main components: 
The vector-part {f^i/\J~^ + ^.j/ is designed to always point to the outside of the 
singularity ellipsoid, see Eig. 7.6. q* is a vector that points from focus z to a particular 
point in loop three-momentum space. This is indicated by the dashed lines. Keep in 
mind that, despite its name, a ‘focus’ is not a focal point of the ellipse/ellipsoid, but 
rather refers to the foci introduced in Section 6.1, i.e. focus i = —kj. Dividing by 
rescales the vector to length 1. The same procedure is repeated with q^. Adding the 
two together results in the angle between q^ and the resulting vector part always being 


^ Strictly speaking, there is a third condition: 

The deformation must vanish at the position of soft or collinear singularities: 

This point is of importance for the matching of soft and collinear singularities between real and virtual 
corrections. If the deformation shifts those singularities, alongside everything else, the cancellation will 
be spoilt. However, in the scope of this thesis, we are only dealing with finite diagrams. 
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< 7r/2. By choosing all the scaling parameters Ajj < 0 for all possible combinations {ij} 
we satisfy the first condition. 



Figure 7.6: Two-dimensional slice of the singularity ellipsoid of dual contribution i 
at position j. The resulting vector gives the orientation of the vector part. 


Inside the singularity ellipsoid the two vectors qi/y qf and q^/y q^ cancel totally along 
the major axis of the ellipsoid. This cancellation in the vicinity of the center of the 
singularity ellipsoid is shown in Fig. 7.7. In order to be able to plot the situation, the 
graphs are restricted to 1-1-2 dimensions. Keep in mind that the imaginary part of the 
plotted deformation is a two component vector where each component is a function of 
ix and iy. Therefore, the graphs 7.7a and 7.7b indicate how far the deformation goes in 
the direction of the imaginary part of the ix- and I’y-axis. 



(a) X- component 



(b) y-component 


Figure 7.7: Imaginary part of the deformation in 1-1-2 dimensions. The vector part 
causes the deformation to flatten inside the singularity ellipsoid. 
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7.5. Contour deformation 


The suppression factor ex.p{—Gf)^{qj;qi)/Aij) fulfills the second condition. At the po¬ 
sition of the singularity, Gf^{qj]qi) is 0 and thus the suppression factor reaches its 
maximum. Far away from the singularity, —Gfj^{qj',qi) is a large negative value and 
thus the exponential is close to zero. This behaviour is demonstrated in Fig. 7.8. 



Figure 7.8: Volcano-shaped suppression factor in 1-1-2 dimensions. The crater line is 

exactly the singularity ellipse. 



Xij is a scaling factor analogous to A in Section 7.5.1. Aij is the width of the deformation. 
The indices ij in Xij and Aij indicate that those parameters can be chosen individually 
for each contribution to the deformation for optimization purposes. Then we sum over 
the entire group and arrive at: 


in 


i ^ Xij 
z jEgroup 



(7.18) 


where, of course, i is the imaginary unit. Finally, we add the imaginary contribution to 
the loop momentum to make the deformation complete. 


i 


i' 


i + iK 


(7.19) 


The corresponding Jacobian can be calculated analytically. 
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7.6 Choosing the parameters 

As mentioned before, the parameters and Aij are the scaling and the width of 
the deformation. Typically, \ij is a dimensionless negative number and Aij a positive 
number of dimension q^] setting all Ajj = 0 completely switches off the deformation. 
Although they can all be chosen differently, selecting Ajj = A = —0.5 and Aij = A = 10® 
produces very good results for most of the examples that we are presenting within this 
thesis. As a rule of thumb, Aij ~ lO^m^ where is a measure for the energy^ of the 
process in consideration. 

7.6.1 The scaling parameter A 

The parameter Ajj determines how far an individual contribution to the deformation 
goes around a pole. This is illustrated in Fig. 7.9, in which for the sake of clarity we 
restricted ourselves to 1 + 1 dimensions and a deformation that consists of only one 
contribution. 



Figure 7.9; The deformation at different values of A. The dashed vertical lines mark 

the positions of the singularities. 


From the dotted to the dashed and finally to the solid line, Ajj = A has been doubled 
every time while keeping the width constant. Therefore, the peak of the solid curve is 
four times higher (lower) than the dotted one. Between lx = 0 and lx = ^ there is a 
region where the deformation vanishes. This is the cancelling effect of the vector part 
mentioned in Section 7.5.2. In 1+1 dimensions this leads to a total cancellation. It is 
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equivalent to slicing Fig. 7.7 along the axis that connects the two foci of the singularity 
ellipse. 

There is one further aspect to it which is inspired by [89]. On the one hand, we would 
like to choose A as large as possible to stay away from the singularity. One the other 
hand, we would like to make A as small as possible in order to ensure that we do not 
enclose other (unwanted) singularities. To escape from this dilemma, consider the on 
shell energy with the deformation I + zAk plugged in: 

+ 2i\K • qi + q2 + - iO . (7.20) 

Here we have introduced an overall factor A. It scales the entire deformation as opposed 
to the \ij which scale individual contributions to the deformation (see Eq. (7.18)). 
Since > 0, this puts a constraint on A. We set the square root to zero and solve the 
quadratic equation in A: 

A± = iXi ± (7.21) 


with 




qz • 1 ^ 

9 ’ 


Yi = 


qf +mf 

k2 


(7.22) 


where W < 0 and 17 > 0. If Xf > 1), the poles of lie on the imaginary A-axis and 
the real value of A may take any value. If however, X'f < Yi and Xj —?■ 0 we have a pole 
at A = XYi. Thus we define A in the following way: 


A 


'\XYi if Xl<Yil2, 

< ^x1 - Y if y,/2 < xl < y,, 
.1 if Y,<Xl 


(7.23) 


After performing this check for every dual contribution in the deformation group, we end 
up with a set of “safety As” out of which we pick the smallest. Using this procedure, we 
make sure that our deformation does not come too close to other singularities originating 
from the dual delta function. 


7.6.2 The width of the deformation 

The parameter Aij determines how broad or narrow the peaks of the deformation con¬ 
tributions are. This is illustrated in Fig. 7.10. Again, to improve clarity, the plot is 
restricted to 1 + 1 dimensions and the deformation used has only a single contribution 
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to it. This time the width has been increased by a factor of five from the dotted to the 



Figure 7.10: The deformation at different values of A. The dashed vertical lines mark 

the positions of the singularities. 


dashed curve and another factor five from the dashed to the solid curve. The scaling 
parameter A has been kept constant. 

As mentioned, setting Aij = A = 10® works fine for the majority of momentum config¬ 
urations that we tested. Nonetheless, there are some cases where a bigger (or smaller) 
choice of width is advised to get good precision. To understand why, let us have a look 
at the dual propagator in its rewritten form: 


GniquQj) 


+ kjiflY 



(7.24) 


We see that the inverse of the dual propagator is proportional to momentum G~j^ ~ q. 
In the suppression factor the inverse of the dual propagator gets squared. Hence the 
exponent of the suppression factor is ~ q'^ jA. In order to avoid the deformation getting 
to wide or too narrow, it is necessary choose the width proportional to the physical 
energy scale of the external momenta. 


7.7 Numerical Results 

Now we are in possession of all the necessary tools to make the numerical integration 
work. We have a code which runs on a desktop machine with an Intel i7 (3.4GHz) 
processor with 8 cores and 16 GB of RAM. Memory consumption is negligible. 












7 . 7 . Numerical Results 


7.7.1 Details on the implementation 

The program is written entirely in C++ and uses the Cuba library [71] as a numerical 
integrator. This means it can be run on any machine on which C++ and the Cuba 
library are available. As input the user has to specify the number of external legs, the 
external momenta themselves, the internal masses and, if necessary, the parameters of 
the contour deformation. The momenta and masses can be read in from a text file. 
With regard to the numerical integrator the user can choose between Cuhre [95, 96] and 
VEGAS and give the desired number of evaluations. At this stage, the program is ready 
to be compiled and executed. Schematically speaking, it performs the following steps: 

1. Read in momenta and masses. 

2. Check where ellipsoid singularities occur. 

3. Check where hyperboloid singularities occur and group the dual contributions 
accordingly. 

4. Call the integrator, taking into account the number of external legs as well as the 
results from the previous steps. 

We use MATHEMATICA 8.0 [97] to generate randomized momenta and masses and LoopTools 
2.10 [98] to produce reference values for comparison. The goal is to scan as much of 
the phase-space as possible to make sure the program works properly in all regions. The 
momenta of all the example points and scans of the following sections are collected in 
Appendix C. 

We see two main paths along which the complexity of the calculation grows. 

1. Increasing the number of external legs. In general, going from n —n + 1 legs 
means for the Duality: 

• One extra dual contribution. 

• Since each dual contribution does consist of one more dual propagator, there 
are much more possibilities for ellipsoid singularities to occur. Thus the 
deformation picks up more contributions. 

Therefore, we investigate multi-leg scalar integrals up to the pentagon level in this 
chapter. 

2. The presence of tensor numerators. Having tensors of increasing rank in the numer¬ 
ator will render the function to be integrated more complex. We address tensor 
integrals in Chapter 8 although it does not affect the singular behaviour of the 
integral. 
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Our setup is such that we run points asking for a fixed number of evaluations. Then we 
modify the parameters of the deformation for best results.^ 

7.7.2 Scalar Triangle 

We consider first infrared finite triangle (three external legs) scalar integrals. Momentnm 
configurations that do not need deformation (i.e. whose loop integral is pnrely real) are 
integrated in about 0.15 seconds with a precision of at least 4 digits. 5 • 10^ evalnations 
are snfficient to achieve this result. 



Point 1 

Point 1 Error 

Point 2 

Point 2 Error 

LoopTools 

-5.85694E-5 

0 

-3.39656E-7 

0 

Loop-Tree Duality 

-5.85685E-5 

2.4E-9 

-3.39688E-7 

5.3E-11 


Table 7.1: Example for two non-deformation phase-space points. Since there is no 
deformation, all values are purely real. 


Point 1 of Table 7.1 has all internal masses equal while Point 2 has three different in¬ 
ternal masses. Momenta and masses were chosen randomly between —100 and -1-100. 
This even allows for nnphysical momentum configurations, but at this stage we want to 
test stability and precision regardless of whether we are in the physical region or not. 
For example, althongh Triangle 2 represents a sitnation that is not realized in nature, 
it is computed without any problems. From the program’s point of view there is no 
difference; the internal masses are mere parameters. 

For momentnm configurations that reqnire deformation the function to integrate is more 
complicated, hence we have to evaluate it more often. With 10® iterations, the calcu¬ 
lation time increases to around 2.5 seconds. Four digits of precision are achievable by 
optimizing the parameters of the deformation. 



Real Part 

Real Error 

Imaginary Part 

Imaginary Error 

LoopTools P.3 

5.37305E-4 

0 

-6.68103E-4 

0 

Loop-Tree Duality P.3 

5.37307E-4 

8.6E-9 

-6.68103E-4 

8.6E-9 

LoopTools P.4 

-5.61370E-7 

0 

-1.01665E-6 

0 

Loop-Tree Duality P.4 

-5.61371E-7 

7.2E-10 

-1.01666E-6 

7.2E-10 


Table 7.2: Example for two phase-space point that need deformation. With a defor¬ 
mation applied, the results are complex. 

^Part of the results which we are going to show in the following will be published in a forthcoming 
paper [3]. 
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7 . 7 . Numerical Results 


Similar to the non-deformation points before, Point 3 of Table 7.2 has all internal masses 
equal whereas in Point 4 all three of them have different values. Regardless of kinemat¬ 
ics, the Loop-Tree Duality is capable of producing accurate results in both situations. 



0 H E 
HOE 

0 0 0 


Figure 7.12: Singularities of a triangle graph in loop-momentum space and singularity 

scheme for comparison. 

Fig. 7.12 gives an additional view on a triangle example phase-space point: Each of the 
cubes in the first line represents one dual contribution. Drawn are the ellipsoid (violet 
surfaces) and hyperboloid (orange surfaces) singularities in loop-momentum space. The 
blue dots are the positions of the foci of the on-shell hyperboloids, i.e. —k,, i £ {1, 2, 3}. 
In the second line we have the singularity scheme. This means each line of the singular¬ 
ity scheme corresponds to one box. 

Another important check is the scan around threshold. This region is usually numeri¬ 
cally unstable. The Loop-Tree Duality excels here, because the algorithm does not have 
to do extra work. In Fig. 7.13 all internal masses are equal, i.e. rm = m, i £ {1,2,3}, 
and the center-of-mass energy s was kept constant while the mass m was varied. As you 
can see from the plots, there is no drop in precision around threshold. Calculation time 
remains constant, as well. 
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Figure 7.13: Scan of the region around threshold. The red curve is LoopTools, the 
blue points are the Loop-Tree Duality. 
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7.7.3 Scalar Boxes 

We repeated the analysis for the box (four external legs) case. The main difference is 
that now there are four dual contributions with three dual propagators each. 

To get good precision, we chose to evaluate boxes that need deformation for 4 • 10® 
times, while the number of evaluations for non-deformation phase-space points was kept 
at 5 • 10“^ evaluations, the same as in the triangle case. This is reflected in the program 
runtimes. Deformation points take 16 seconds, non-deformation points take about 0.25 
seconds, respectively. 

While it is practically guaranteed to get the non-deformation points with good precision, 
the quality of results depends on the proper choice of parameters for points with de¬ 
formation. Therefore, we focus our attention towards such points by giving many more 
examples for them. 

A very good method to do so is to perform a mass-scan as we already did in the tri¬ 
angle case. From Fig. 7.16 to Fig. 7.18 we present scans for the following momentum 
configurations: 

• All internal masses equal. The center-of-mass energy s is kept constant while the 
mass is varied. This is shown in Fig. 7.16. 

• Two adjacent internal masses equal and the two opposing masses equal, i.e. mi = 
m 2 and m 3 = m 4 . mi and therefore also m 2 is varied, the other momenta and 
masses are kept constant. This is shown in Fig. 7.17. 

• Scan in which the Mandelstam variable t is varied. This is realized by varying 

while keeping P 3 constant. Of course, due to momentum conservation, this involves 
7*4 = {pi + P 2 not being constant. This is shown in Fig. 7.18. 

Concrete values for each of those situations are found in Table 7.3. 

The scans show how the program performs in certain slices of the phase space. This is 
important because we want the program to be generically applicable to any momentum 
configuration. From the plots and tables, you can see that the program deals well with 
all kinds of boxes, even when many different kinematical scales are involved. 

Point 5 and 7 of Table 7.3 correspond to a momentum configuration, in which all four 
internal masses are equal. In Point 6 and 8 all masses are different. In Point 9 two 
adjacent internal lines have equal masses as well as the two opposing ones. Point 10 
represents a situation in which opposite lines have equal masses. 
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Point 5 

Point 5 Error 

Point 6 

Point 6 Error 

LoopTools 

2.15339E-13 

0 

1.39199E-11 

0 

Loop-Tree Duality 

2.15319E-13 

5.2E-17 

1.39199E-11 

6.3E-16 


(a) Two non-deformation phase-space points. 



Real Part 

Real Error 

Imaginary Part 

Imaginary Error 

LoopTools P.7 

-2.38766E-10 

0 

-3.03080E-10 

0 

Loop-Tree Duality P.7 

-2.38798E-10 

8.2E-13 

-3.03084E-10 

8.2E-13 

LoopTools P.8 

-4.27118E-11 

0 

4.49304E-11 

0 

Loop-Tree Duality P.8 

-4.27127E-11 

5.3E-14 

4.49301E-11 

5.3E-14 

LoopTools P.9 

6.43041E-11 

0 

1.61607E-10 

0 

Loop-Tree Duality P.9 

6.43045E-11 

8.4E-15 

1.61607E-10 

8.4E-15 

LoopTools P.IO 

-4.34528E-11 

0 

3.99020E-11 

0 

Loop-Tree Duality P.10 

-4.34526E-11 

3.5E-14 

3.99014E-11 

3.5E-14 


(b) Four phase-space points with different kinematics that need deformation. 


Table 7.3: Examples for boxes with and without deformation. 


Similar to the triangle case, Fig. 7.15 shows the singularities in loop three-momentum 
space. It becomes immediatly obvious that, only by going to box graphs, the level of 
complexity has grown substantially. Because the hyperboloid singularities couple dual 
contributions together, an individual contribution receives a much more complicated 
deformation. 
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0 

H 

H 

E 

H 

0 

H 

E 

H 

H 

0 

E 

0 

0 

0 

0 


Figure 7.15; Singularities of a box graph in loop-momentum space and singularity 

scheme for comparison. 


In Fig. 7.15, the dots give the locations of the foci —kj, the surfaces are ellipsoid and 
hyperboloid singularities. Note how the hyperboloids always appear pairwisely across 
the dual contributions. 
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0 


-2.xl0“‘“ 


-4.xl0“‘“ 
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Figure 7.16: Mass-scan of a box integral. The red curve is LoopTools, the blue points 

are the Loop-Tree Duality. 




In Fig. 7.16, two thresholds are passed at 2m/y^ = 0.65 and 1. From right to left, the 
number of ellipsoid singularities grows by one after each threshold from one to three. 
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-2.x 10 


-4. X 10 


-6.x 10 


-8.x 10 


-l.xlO 






V7 


Figure 7.17: Mass-scan of a box integral. The red curve is LoopTools, the blue points 

are the Loop-Tree Duality. 


In Fig. 7.17, one threshold is crossed at 2m\l\fs ~ 0.65. From right to left the number 
of ellipsoid singularities goes from one to two. 
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t 



t 


Figure 7.18: t-scan of a box integral. The red curve is LoopTools, the blue points are 

the Loop-Tree Duality. 


In Fig. 7.18, the Mandelstam variable t is being varied. Two thresholds are passed, 
one at t ~ 10^ and one at t ~ 150 • 10^. From left to right, the number of ellipsoid 
singularities increases by one after each crossing of a threshold from zero to two. 
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7.7.4 Scalar Pentagons 


After the evaluation of triangles and boxes, the next step is to check pentagon (five 
external legs) graphs. Complexity grows once more, which means: 


• A dual integral now consists of five dual contributions. 


• Each dual contribution contains one additional (= four in total) dual propagator. 


Therefore, we had to increase the number of evaluations: Non-deformation points are 
evaluated 10® times which takes approximately 0.5 seconds. Points with deformation 
demand 5 • 10® evaluations to maintain the level of precision of the triangles and boxes. 
This results in an average calculation time of 28 seconds. This enhanced complexity 
is especially well-illustrated by the following Figure in which the singularities of the 
different dual contributions are plotted in loop-momentum space: 







0 
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0 
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Figure 7.20: Singularities of a box graph in loop-momentum space and singularity 

scheme for comparison. 


Again, dots are foci of the on-shell hyperboloids, surfaces are ellipsoid and hyperboloid 
singularities. From the singularity scheme we can read off that all five dual contributions 
are coupled together due to the interaction of the hyperboloids. 
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Point 11 

Point 11 Error 

Point 12 

Point 12 Error 

LoopTools 

Loop-Tree Duality 

-1.24025E-13 

-1.24027E-13 

0 

1.6E-17 

-1.48356E-14 

-1.48345E-14 

0 

l.lE-17 


(a) Two non-deformation phase-space points. 



Real Part 

Real Error 

Imaginary Part 

Imaginary Error 

LoopTools P.13 

1.02350E-11 

0 

1.40382E-11 

0 

Loop-Tree Duality P.13 

1.02353E-11 

l.OE-16 

1.40385E-11 

l.OE-16 

LoopTools P.14 

7.46345E-15 

0 

-9.13484E-15 

0 

Loop-Tree Duality P.14 

7.46309E-15 

6.1E-18 

-9.13444E-15 

6.1E-18 

LoopTools P.15 

6.89836E-15 

0 

2.14893E-15 

0 

Loop-Tree Duality P.15 

6.89848E-15 

6.5E-18 

2.14894E-15 

6.5E-18 


(b) Three phase-space points that need deformation. 


Table 7.4: Examples for pentagons with and without deformation. 


Table 7.4 displays a collection of pentagon sample results for different kinematical con¬ 
figurations. In Points 11 and 13 all internal masses are equal; in Point 14 they are all 
distinct from each other and in Point 15 we have mi = m 2 = m 3 7 ^ m^ = ms. Again, 
the Loop-Tree Duality shows its robustness by producing accurate results regardless 
of the kinematical situation. This statement is further supported by various scans we 
performed: 

• All internal masses equal. The center-of-mass energy s is kept constant while the 
mass gets varied. This is shown in Fig. 7.21. 

• All five internal masses different. Mass mi gets varied, the rest is kept constant. 
This is shown in Fig. 7.22. 

• Scan in which the center-of-mass energy s is varied. This is realized by varying pi. 

Of course, due to momentum conservation, this involves pi = {pi P 2 not 

being constant. This is shown in Fig. 7.23. 

From the plots, you can see that the Loop-Tree Duality is able to achieve good precision 
for pentagons as well. Even situations in which many ellipsoid singularities are involved, 
are handled well by the program. For example, the number of ellipsoid singularities in 
Figure 7.21 increases from two to five when going from right to left. 
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Figure 7.21: Mass-scan of a pentagon integral. The red curve is LoopTools, the blue 

points are the Loop-Tree Duality. 




Three thresholds are passed in this scan. From right to left, we start at two ellipsoid 
singularities and arrive at five. 
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V7 

Figure 7.22: Mass-scan of a pentagon integral. The red curve is LoopTools, the blue 

points are the Loop-Tree Duality. 


One threshold is crossed at 2mi/\/s ~ 0.2. From right to left, the number of ellipsoid 
singularities increases from five to six. 
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Figure 7.23: Energy-scan of a pentagon integral. The red curve is LoopTools, the 
blue points are the Loop-Tree Duality. 


In this scan, we pass three thresholds at s ~ —8.5 • 10^,—13.5 • 10^ and —21 • 10^ 
which divide the graph into four zones. From right to left, we start with zero ellipsoid 
singularities in the first zone, then we have one in the second zone, two in the third zone 
and finally one in the last zone. 
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Chapter 8 


Tensor Integrals 


Up to this point, we have shown that the Loop-Tree Duality is capable of dealing with 
all kinds of kinematical configurations. The only constraint that we have made was 
to limit ourselves to scalar integrals. In Chapter 3, we argued that the inclusion of 
numerators different from 1 should not be an issue for the Loop-Tree Duality because 
it acts only on the denominator of the loop integral. Hence, in this chapter, we are 
going to relieve this restriction and repeat some of the analysis of the previous chapter 
with tensor integrals. Additionally, we are going one step further by calculating even 
hexagons, which are usually tough to compute, but from the perspective of LTD are a 
straightforward generalization. 
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Chapter 8. Tensor Integrals 


8.1 The influence of numerators 


To start off, we restate the Duality Theorem for scalar integrals (cf. Eq. (6.3)): 


L^^\pi,P2,---,Pn) = / S{qi) Goiqiqj) , 


( 8 . 1 ) 




j&ai 


where 

Gd{ qt , Q]) ^2 _ ^2 _ ^0 r]{qj -qi) ' 

If our one-loop integral features a non-trivial numerator AA(t', {pj}). 
Tree Duality Theorem takes the form 


( 8 . 2 ) 


Then, the Loop- 


L^^\pi,P2 ,---,Pn) = j n GD{qi]qj)M{(.,{pi}) . (8.3) 

i£ai ^ jSal 


While the numerator is formally left unchanged, there actually is a practical impact. 
The presence of the dual delta function demands g-which is equivalent 
to having the energy component of the loop-momentum forced to 


4 


~ki,0 + 


q2 + m2 


(8.4) 


In other words, whenever we perform a single cut of a Feynman graph, the numerator 
has to be evaluated at the position of the cut which is fixed by the dual delta function. 
As a direct consequence, the numerator takes a different form in each dual contribution. 

Another important aspect to take into consideration is the cancellation of singularities 
among dual contributions. In Section 6.2, we argued that numerators do not spoil the 
cancellation of the hyperboloid singularities. Here, we would like to make explicitly show 
why. A typical numerator is a product of scalar products of the form “loop-momentum 
contracted with external momentum”; i- pk- Let us see what happens to a single factor 
when it hits the singularity. Note first, that the hyperboloid singularity is given by 
Eq. (6.8) which we rewrite in the more suitable form 

■ (8-5) 

Using Eq. (8.4), the loop-momentum i contracted with some external momentum is: 

^ * Pk |i-th cut — ( ki -\- q )P/c,0 ^ ‘ Pfc 

— ( kj ~\~ Qj^Q )Pk,0 ^ * P/c — ^ * Pk Ij-th cut (^'6) 
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8.1. The influence of numerators 


where we have used Eq. (8.5) for the first equality on the second line of Eq. (8.6). It 
means that the numerators of two dual contributions i and j take the same value at their 
common pole, thus leaving the cancellation of hyperboloid singularities intact. This is 
an important property to take advantage of, because it allows us to straightforwardly 
apply the Loop-Tree-Duality to such diagrams without any additional effort. 


8.1.1 An illustrative example 


To make things even more explicit, consider the following simple example. We take the 
triangle integral of Section 3.5.2 and give it the numerator Af{i,pi,P 2 ) = ^ • Pi- 

r 3 

1= i-pi nGF(%) (8.7) 

T=i 

Pi,P 2 and p 3 = —pi —p 2 are the external momenta. Since we are dealing with a triangle, 
obviously N = 3. Consequently, the Loop-Tree Duality translates the Feynman integral 
into three dual contributions. 


1= h{qi)GD{qi]q2)GD{qr,qfl)^-pi 
Je 

+ / GD{q2]qi)^{q2)GD{q2]qii)Gpi 

Je 

+ / GDiq3-,qi)GD{q3;q2)Siq3)i-pi 


( 8 . 8 ) 


This time, we investigate how the numerator affects the calculus. The dual delta func¬ 
tions of the different contributions fix the energy component of the numerators to dif¬ 
ferent values. Hence the scalar product becomes 


i-pi^ i-pi,o + \/{£ + Pi)^ + ml) pifl -i -pi 
i-pi^ i-pi,o - P2,o + \l{(- + ^1+ P 2 Y + ml) pifl 
£-pi^ + mlpifl -£ pi 


first contribution 
£. ■ Pi second contribution 

third contribution (8.9) 


This means, that finally we will arrive at with three dual contributions similar to the 
ones of Eq. (3.70). 

In the next sections, we present the results for multi-leg tensor integrals. We only 
consider IR- and UV-finite diagrams. The explicit momenta of all the points and scans 
presented within this chapter are found in Appendix C. 
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Chapter 8. Tensor Integrals 


8.2 Tensor Triangles 

In Chapter 7, we convinced ourselves that the Loop-Tree Duality can produce accurate 
results for scalar one-loop integrals independent of the kinematical configuration. Now, 
we focus on checking various different numerator functions within a similar analysis. 
We start with infrared-finite one-loop triangle graphs. The simplest non-trivial nu¬ 
merator function possible is the loop-momentum i contracted with one of the external 
momenta pi or p 2 , i.e. i ■ pi oi i ■ p 2 - Since we are limiting our discussion to UV-finite 
graphs, these are the only possible numerators at the triangle level. 

The number of evaluations per phase-space point is the same as in the scalar case; this 
means 5 • 10^ evaluations for non-deformation points and 10® evaluations for phase-space 
points with deformation. Due to the presence of numerators, the integrand function is 
a bit more complex. Hence the calculation time, compared to scalar triangles increases 
from 2.5 seconds to 3.5 seconds for deformation-points. For non-deformation points 
there is hardly any measurable difference, i.e. they stay at around 0.15 seconds 
A compilation examples is shown in Table 8.1. Points 16 and 18 correspond to kinemat¬ 
ical situations in which all internal masses are equal and the numerator is i-p 2 , whereas 
Points 17 and 19 have all internal masses are different from each other and numerator 
i-pi. 



Point 16 

Point 16 Error 

Point 17 

Point 17 Error 

LoopTools 

-1.07284E-2 

0 

-1.59964E-3 

0 

Loop-Tree Duality 

-1.07281E-2 

5.4E-6 

-1.59985E-3 

6.6E-7 


(a) Two non-deformation phase-space points. 



Real Part 

Real Error 

Imaginary Part 

Imaginary Error 

LoopTools P.18 

-2.64773E-3 

0 

1.37469E-2 

0 

Loop-Tree Duality P.18 

-2.64726E-3 

5.5E-6 

1.37448E-2 

5.5E-6 

LoopTools P.19 

-1.19501E-2 

0 

1.35834E-3 

0 

Loop-Tree Duality P.19 

-1.19511E-2 

3.6E-5 

1.35859E-3 

3.6E-5 


(b) Two phase-space points that need deformation. 

Table 8.1: Examples for triangles involving numerators, with and without deforma¬ 
tion. 


In the triangle mass-scan of Fig. 8.1, we varied the mass in the same way as we did in 
the previous chapter. Three thresholds are passed at 2m/\/s Ri 0.15, 0.45 and 1. 
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Imaginary part Real part 


8.2. Tensor Triangles 





V7 

Figure 8.1: Mass-scan of a triangle with numerator l-p^. The red curve is LoopTools, 
the blue points are the Loop-Tree Duality. 
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Chapter 8. Tensor Integrals 


8.3 Tensor Boxes 

Again, we repeat the analysis of the triangles at the box-level. We only want to deal 
with UV-finite diagrams, thus tensors up to rank three are allowed in the numerator. 
Since having a rank-one tensor is quite similar to the triangle case, we will focus on rank 
two and rank three situations in order to have a more challenging test for the program. 
We use the same number of evaluations as we did for scalars, this means 4-10® evaluations 
for deformation points and 5 • 10^ evaluations for points that do not need deformation. 
As far as deformation points are concerned, calculation times differ according to the rank 
of the tensor function in the numerator. Their general order of magnitude is around 20 
seconds, the span between a rank zero and a rank three tensor amounts to 8 seconds. 
Non-deformation points take about 0.25 seconds, almost independently of the rank of 
the numerator function. 



Point 20 

Point 20 Error 

Point 21 

Point 21 Error 

LoopTools 

-3.42913E-4 

0 

9.64909E-b0 

0 

Loop-Tree Duality 

-3.42905E-4 

5.9E-8 

9.64924E-b0 

4.6E-3 


(a) Two non-deformation phase-space points. 



Real Part 

Real Error 

Imaginary Part 

Imaginary Error 

LoopTools P.22 

1.33032E-3 

0 

-1.07780E-3 

0 

Loop-Tree Duality P.22 

1.33033E-3 

5.6E-7 

-1.07779E-3 

5.6E-7 

LoopTools P.23 

-2.15448E-b2 

0 

-1.10792E-b2 

0 

Loop-Tree Duality P.23 

-2.15451E-b2 

8.7E-2 

-1.10789E-b2 

8.7E-2 


(b) Two phase-space points that need deformation. 

Table 8.2; Examples for boxes involving numerators of tensor rank two and three, 

with and without deformation. 


In Table 8.2, Points 20 and 22 are boxes with all internal masses equal and numerator 
i ■ Pi X i ■ p 2 . Points 21 and 23 represent situations in which all internal masses are 
different and the numerator has the form i ■ pi x i ■ p^ x i ■ p 4 ^. 

We also did a scan of the Mandelstam variable t, see Fig. 8.2, for which we took a 
box with three scalar products in the numerator, i ■ pi x i ■ p 2 x i ■ p^, and varied p^ 
while keeping constant. Since p^ also appears in the numerator, both numerator and 
denominator are affected by this scan. 
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8 . 3 . Tensor Boxes 




Figure 8.2; Scan of the region around threshold. The red curve is LoopTools the blue 

points are the Loop-Tree Duality. 
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Chapter 8. Tensor Integrals 


8.4 Tensor Pentagons 

Next, we investigate pentagon tensor integrals at the one-loop level with numerators 
up to rank three. The number of evaluations is chosen to be the same as in the scalar 
case, i.e. 10^ times for non-deformation and 5 • 10® times for phase-space points that 
require deformation. This results in calculation times of 0.7 seconds and 33 seconds, 
respectively. The runtime difference between a point with a numerator of rank zero and 
rank three is negligible for non-deformation points and about 10 seconds for deformation 
points. 

Table 8.3 shows a selection of sample points. Point 24 and 26 feature the rank two 
numerator x 7 • p 4 while Point 25 and 27 have the numerator i-p^x i-p^x i-p^. In 
all points all internal masses are equal. At first glance, Point 27 seems to be weaker but 
it actually contains six ellipsoid singularities whereas the other points usually have two 
to three. We include this point to demonstrate that the program holds well together 
even under such challenging circumstances. 



Point 24 

Point 24 Error 

Point 25 

Point 25 Error 

LoopTools 

-1.86472E-8 

0 

1.74828E-3 

0 

Loop-Tree Duality 

-1.86462E-8 

2.6E-12 

1.74808E-3 

2.8E-7 


(a) Two non-deformation phase-space points. 



Real Part 

Real Error 

Imaginary Part 

Imaginary Error 

LoopTools P. 26 

-1.68298E-6 

0 

1.98303E-6 

0 

Loop-Tree Duality P. 26 

-1.68298E-6 

7.4E-10 

1.98299E-6 

7.4E-10 

LoopTools P. 27 

-8.34718E-2 

0 

1.10217E-2 

0 

Loop-Tree Duality P. 27 

-8.34829E-2 

7.5E-5 

1.10119E-2 

7.5E-5 


(b) Two phase-space points which need deformation. 

Table 8.3: Examples for pentagons involving numerators of tensor rank two and three, 

with and without deformation. 


We also performed several scans; a sample is presented in Figure 8.3. In that scan, 
similarly to what we have done with scalar pentagons, we varied pi and thus the center- 
of-mass energy s = {pi+P 2 )- The corresponding numerator function is i-pi x£-p 2 xi-p^, 
which means that both numerator and denominator take part in the scan. From the 
plot, we can see that the Loop-Tree Duality is able to pass this challenging test. 
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8 . 4 . Tensor Pentagons 




Figure 8.3: Scan of the region around threshold. The red curve is LoopTools, the 
blue points are the Loop-Tree Duality. 
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Chapter 8. Tensor Integrals 


8.5 Tensor Hexagons 

We also analysed hexagon (six external legs) tensor integrals, for which we raised the 
number of evaluations for non-deformation points to 10® and for deformation points to 
8-10®. The corresponding calculation times are 8 and 75 seconds, respectively. Since we 
evaluate non-deformation points ten times more compared to pentagons, the different 
numerators have a measurable impact on the calculation time. The difference between 
rank zero and rank three is 6 seconds and 15 seconds for deformation points. 

LoopTools can provide reference values only up to pentagons. Instead, we used the 
program SecDec version 3 [38] for comparison. 

Again, we present a selection of sample points in Table 8.4. Point 28 and 30 feature the 
rank-one numerator i-pi, in the former all internal masses are different and in the latter 
they are all equal. Point 29 has six distinct internal masses and the numerator function 
£ • P 2 X ^ • P 4 X ^ • P6) Point 31 possesses the numerator i ■ P 2 x i ■ p^ and six different 
masses, as well. Finally, in Point 32, which exhibits the numerator i ■ p/^ x i ■ p^ x i ■ pQ, 
all momenta are distinct form each other. 



Point 28 

Point 28 Error 

Point 29 

Point 29 Error 

SecDec 

-1.21585E-15 

1.2E-19 

4.46117E-9 

3.7E-13 

Loop-Tree Duality 

-1.21552E-15 

3.5E-18 

4.46136E-9 

2.6E-15 


(a) Two non-deformation phase-space points. 



Real Part 

Real Error 

Imaginary Part 

Imaginary Error 

SecDec P. 30 

1.01359E-15 

2.3E-19 

2.68657E-15 

2.6E-19 

Loop-Tree Duality P. 30 

1.01345E-15 

1.3E-18 

2.68633E-15 

I.3E-I8 

SecDec P. 31 

2.45315E-12 

2.4E-16 

-2.06087E-12 

2.0E-16 

Loop-Tree Duality P. 31 

2.45273E-12 

7.3E-15 

-2.06202E-12 

7.3E-15 

SecDec P. 32 

-2.07531E-6 

1.9E-10 

6.97I58E-7 

5.6E-1I 

Loop-Tree Duality P. 32 

-2.07526E-6 

7.5E-13 

6.97192E-7 

7.5E-13 


(b) Two phase-space points that need deformation. 

Table 8.4: Examples for pentagons involving numerators of tensor rank two and three, 

with and without deformation. 
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Chapter 9 


Summary and Future Work 


9.1 Summary 

One of the main difficulties in calculating scattering amplitudes at NLO and higher is 
the cancellation of infrared singularities among real and virtual corrections. The Loop- 
Tree Duality aims to solve that problem by recasting the virtual corrections in a form 
which closely resembles the real ones. This appealing property motivates the idea of 
combining the two, thus treating them simultaneously in a common Monte Carlo event 
generator. 

Initially the LTD was introduced for scalar one-loop integrals [48] and in later works 
its applicability has been expanded by a large margin via a systematic procedure has 
been established to calculate Feynman graphs with an arbitrary number of loops [49] . In 
Chapter 5 and in [I] we have shown how to deal with loop integrals involving higher order 
poles. Instead of applying the Residue Theorem for higher order poles, using Integration 
By Parts techniques reduced the integrals featuring such poles to integrals with simple 
poles. After that we could go ahead and straightforwardly apply the known formulae. 
At that point, the major theoretical questions were answered and the important was to 
check the efficiency of the method. To that direction, as a preparatory step (Chapter 6 
and [2]), we investigated the singularities of the dual contributions and found a partial 
cancellation of singularities among contributions happening at the integrand level. In 
the case of massless diagrams, we were able to show that the infrared singularities are 
conhned to a hnite region of loop-momentum space which can then be mapped to the 
real corrections. With this information at hand, we did a numerical implementation of 
the Loop-Tree Duality in the form of a computer program written in C+-I-. Being a 
hrst step towards the combined treatment, the program calculated one-loop integrals. 
Due to the partial cancellation of singularities, we only had to deal with singularities 
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of the ellipsoid type. We dealt with the ellipsoid singularities by contour deformation. 
Still, the hyperboloid singularities affected the way we deformed: In order to keep the 
cancellations intact, we grouped dual contributions featuring the same hyperboloid sin¬ 
gularities together. Then, we employed a contour deformation that accounted for all 
the ellipsoid singularities of the entire group. This contour deformation was constructed 
to respect the iO-prescription of the propagators and fall rapidly to 0 where it was not 
needed so that the value of the overall integral would not be altered. Since dual inte¬ 
grals are three-dimensional, we only had to deform in loop three-momentum space. As 
integrators we used Cuhre and VEGAS from the Cuba library. This setup has proven to 
be successful. We were able to calculate (UV- and IR-finite) scalar with good precision. 
This has been illustrated by showing individual points as well as several scans through 
the phase space, in which we picked one mass or momentum to be varied while the other 
ones were kept constant. In Chapter 8 we also demonstrated that the Loop-Tree Duality 
deals equally well with tensor integrals (numerators different from 1) up to rank three. 
The numerical implementation is still in an early state. Neverthelss, the wide range of 
checks we have performed proves the potential of the method. In particular, the imple¬ 
mentation works well independently from the number of kinematical scales (bilinears, 
trilinears, different masses) involved or whether a phase space point is in the vicinity of 
a threshold or not. 


9.2 Future Work 

Having shown that the Loop-Tree Duality method works and is efficient, the next step 
will be to address any issue the code might run into. Furthermore we would like to 
provide a systematical way to choose the parameters A and the width for optimal re¬ 
sults. We also aim to try the LTD for a real, physical process (i.e. not a toy model), 
for example the six- or eight-photon amplitude. To that end, we will move towards 
treating a case where we combine the real and virtual corrections. The goal is to have 
a highly automated program which calculates the entire cross-section at the one-loop 
level. Finally, the next step in our list is to implement the LTD for two-loop amplitudes. 
This is quite an extensive list, but with the presented work we have laid a good basis. 
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Appendix A 


Duality Nomenclature 


• Loop-Tree Duality: 

The Loop-Tree Duality (in short: Duality) is a Feynman integral transformation 
that maps loop integrals to a sum of phase-space (tree-level) integrals. Within the 
conventions of Chapter 3, the Loop-Tree Duality at one-loop level looks like: 



(AT) 




This transformation is achieved by performing the integration over the energy 
component of the loop integral. The resulting integration runs only over the loop 
three-momentum and is very similar to the real radiation corrections. Thus, it 
encourages the idea of combining the two, treating them simultaneously in a com¬ 
mon Monte Carlo event generator. 

Each summand is called dual contribution and is constructed according to a special 
pattern: One of the internal lines gets on-shell, i.e. it is replaced by a dual delta 
function, while all the other (non-cut) Feynman propagators are promoted to dual 
propagators. This procedure is repeated for every internal line once and the results 
are added together to yield the dual integral. 

• Dual delta function: 

The dual delta function S{qi) is the delta function that appears after applying the 
Loop-Tree Duality: 


6{qi) = 2TTie{qifl)6{qf - m-) = 2-71 i S+{qf - mf) 


It serves several purposes. Formally it is there (cf. Eq. (A.l)) to indicate that the 
energy integration has already been carried out, leaving an integral over the loop 
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three-momentum only. It fixes the £o-component to a certain value determined by 
the argument of the delta function. The theta function, which is included in the 
dual delta function, ensures that the positive energy solution gets picked from the 
argument of the dual delta function. 


• Dual propagator: 

Applying the Loop-Tree Duality involves converting all Feynman propagators to 
dual propagators, except the cut one. There are two representations of the dual 
propagator. 


S{Qi)GD{qi;qj) 


Sjqi) 

qj -mj -iOr]- kji 


and Gniqi^qj) 


where 






’ 

(A.2) 


(A.3) 


The main difference between a dual propagator and a Feynman propagator is the 
fO-prescription. The Feynman propagator depends on the loop four-momentum as 
opposed to the dual propagator which depends on the loop three-momentum. 

The iO-prescription indicates which way to go around the singularities when inte¬ 
grating. Since the integration of the energy component has already been carried 
out, the dual propagator can only depend on the loop three-momentum and the 
iO-prescription gets modified (for details on how and why this happens, see Sec¬ 
tion 3.3). As a direct consequence, the new dual prescription is merely a way of 
bookkeeping on how to go around singularities in loop three-momentum space. 
The form on the left of Eq. (A.2) is a natural result of the derivation of the 
Loop-Tree Duality, the form on the right is a little more suited for numerical 
implement at ions. 


• Dual prescription: 

The modified iO-prescription of the dual propagator is called dual prescription. 
The dual prescription of the dual propagator Goigi', qj) has the form 


-iOrj ■ kji , (A.4) 

where 77 is a future-like four-vector, rf > 0 ,770 > 0- Typically, 77 is chosen to be 
(1, 0, 0, 0). In general, r/ depends on the choice of the coordinate frame. Nonethe¬ 
less, once all dual contributions have been added together, 77 cancels out. In other 
words dual integrals are independent from the choice of 77 whereas individual dual 
contributions do depend on it. At one-loop, the vector kji = qj — qi does not 
depend on the loop momentum and determines the sign of the dual prescription. 
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• Dual contribution: 

A dual contribution is one individual summand of the sum in Eq. (A.l), 

r ~ ^ 

/ HQi) Y{GD{qi;qj) ■ (A.5) 

j=i 

It represents one single cut of a one-loop diagram. The cut line gets replaced 
by a dual delta function, the Feynman propagators of all other internal lines are 
promoted to dual propagators. This is demonstrated explicitly in Section 3.5.2. 
Loop-Tree Duality converts a one-loop integral with N external legs into a sum of 
N dual contributions. While the original one-loop integral runs over the loop four- 
momentum, the integrations of the dual contributions are only three-dimensional 
thus simplifying the calculation. 

Individual dual contributions can feature singularities (hyperboloid singularities, 
cf. Section 7.3) that partially disappear once all dual contributions are combined. 

• Dual integral: 

The dual integral is the end result obtained by applying the Loop-Tree Duality to 
a one-loop integral: 

N . N 

i=i j=i 

It features a couple of appealing properties: Since the ^o-integration has already 
been performed, the resulting dual integral only involves a three-dimensional in¬ 
tegration. In fact, it resembles the form of the real corrections. Therefore, the 
Loop-Tree Duality aims towards a combined treatment of loop and tree contribu¬ 
tions. 
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Proof of the Reduction of Eq. 
( 5 . 18 ) 


Here we solve the system of equations, explicitly, to arrive at Eq. (5.18). We note that 
we are not aiming for a full reduction to a set of master integrals but rather to reduce 
the multiple poles to single poles. Therefore, any integral which has single propagators 
is to be considered known. 

Using the IBPs, Eq. (5.17), on the generic integral ^( 0102030405 ), we get the system of 
equations: 

d — 2ai — 04 — 05 — a 44 '^l“ + 044 ^ 2 ^ = 0 , (h.l) 

01 - 04 + ^so 55 + + 044 +( 1 - - 2-) + ail+(2- - 4’) + ^ 055 +( 2 ^ _ 3 -) = 0 , (B.2) 

044 +(s + 2“ - 3~ - 2 5“) - SO 55 + - 2ail+5“ = 0 , (B.3) 

02 - 04 + 022+(1- - 4 -) + 033+(1- + 2- - 4 - - 2 5 -) + 044+(2- - I”) = 0 , (B.4) 

d - 2 o 2 - 03 - 04 + 033 +(s - 2~) + 044 +( 1 ^ - 2”) = 0 , (B.5) 

a 2 -a 3 + s(a22+-033+ + a44+) + 044+(2"-3"-25“)+033+2"-a22+3" = 0 , (B. 6 ) 

where s = + iO. The appearance of the operator 5“ signals that we have the ISP ii -p 

in the numerator of an integral. As long as these integrals possess single propagators, 
we will not reduce them further but consider them known. We also note that a lot of the 
integrals that appear after setting particular values to the parameters Oj in this system, 
are zero in dimensional regularization (in the massless case). Let us start by setting 
02 = 2, oi = 03 = 04 = 1, 05 = —1 in (B.2). We get: 

F(2111 - 1) - F(1112 - 1) - ^sF(12110) - ^E(llllO) = 0 . (B.7) 
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Taking the sum of (B.5) and (B. 6 ) and setting ai = 1, 02 = 03 = 04 = 1, 05 = 0 we get: 


{d - 4)F(11110) + sF(12110) + sF(11120) - 2F(1112 - 1) = 0 . (B. 8 ) 


Taking the difference between (B.5) and (B. 6 ) and setting ai = 2, 02 = 03 = 04 = 1, 05 = 
0 we get: 


2F(2111 - 1) - sT(12110) + sF(11210) - F(10210) = 0 . (B.9) 


Finally, setting ai = 2, 02 = 03 = 04 = 1, 05 = —1 in (B.l), we get: 


{d - 4)F(2111 - 1) - F(1112 - 1) + F(2012 - 1) = 0 . 


(B.IO) 


The integrals T(10120) and F(11210) and F(2012 — 1), in this system of equations, can 
be computed simply by taking further, appropriate combinations of Eqs. (B.1)-(B.6). 
Setting oi = 02 = 03 = 04 = 1, 05 = 0 in (B.l) we get F(10120) = (3 — d)T(llllO). It 
also holds that F’( 10120 ) = F( 10210 ), as can be seen by making the shifts in the loop 
momenta £2 —> —£i — £2 — P- From Eq. (B.5), by setting oi = 02 = 03 = 04 = 1, 05 = 0 
and using the value of F(10120), we get: sF(11210) = (10 — 3(i)F(11110). Einally, 
adding Eqs. (B.l) and (B.2) and setting 04 = 02 = 1,03 = 2,04 = 1,05 = —1 we get: 
F(2012 — 1) = {d — 3){d — 4)E(11110). The rest of the system of equations (B.7)-(B.10) 
can now be solved sequencially, arriving at: 



(B.ll) 


a result which, after putting d = 4 — 2e, agrees with Eq. (5.18). 
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Momenta of the Example Points 


Here we give the momenta and masses of the different phase-space points and scans 
shown in Chapter 7 and 8. Due to momentum conservation pj\f = — Pi, therefore 

it is sufficient to give only the momenta pi to pn-i- 


C.l Individual points 

Point 1 Pi = {5.23923,-4.18858,0.74966,-3.05669} 

P2 = {6.99881, -2.93659, 5.03338, 3.87619} 
mi = m 2 = m 3 = 7.73358 

Point 2 pi = {13.42254,58.79478,-73.11858,-91.95015} 
P2 = {81.65928, -68.52173, 8.75578, -95.05353} 
mi = 49.97454, m 2 = 86.92490, m 3 = 80.22567 

Point 3 Pi = {10.51284,6.89159,-7.40660,-2.85795} 

P2 = {6.45709, 2.46635, 5.84093,1.22257} 
mi = m 2 = m 3 = 0.52559 

Point 4 pi = {95.77004,31.32025,-34.08106,-9.38565} 
P2 = {94.54738, -53.84229, 67.11107, 45.56763} 
mi = 83.02643, m 2 = 76.12873, m 3 = 55.00359 

Figure 7.12 pi = {44.38942,17.84418,12.70440,-23.67441} 
P2 = {11.62982, -35.11756, -9.52573,1.27635} 
mi = m 2 = m 3 = 7.89824 
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Point 5 


Point 6 


Point 7 


Point 8 


Point 9 


Point 10 


Figure 7.15 


Figure 7.20 


Pi = {31.54872,-322.40325,300.53015,-385.58013} 

P2 = {103.90430, 202.00974, -451.27794, -435.12848} 

P3 = {294.76653,252.88958,447.09194,311.71630} 
mi = m 2 = m 3 = m 4 = 4.68481 

p^ = {50.85428, -55.74613,11.69987, 94.92591} 

P2 = {0.69914, 67.19262, -5.78627, 91.52776} 

P3 = {52.35768,76.32258,43.82222,13.05874} 

mi = 54.29650, m 2 = 53.54058, m 3 = 55.96814, m 4 = 51.74438 

Pi = {62.80274, -49.71968, -5.53340, -79.44048} 

P2 = {48.59375,-1.65847,34.91140,71.89564} 

P3 = {76.75934, -19.14334, -17.10279,30.22959} 
mi = m 2 = m 3 = m 4 = 9.82998 

Pi = {98.04093,77.37405,30.53434,-81.88155} 

P2 = {73.67657, -53.78754,13.69987,14.20439} 

P3 = {68.14197, -36.48119, 59.89499, -81.79030} 

mi = 81.44869, m 2 = 94.39003, m 3 = 57.53145, m 4 = 0.40190 

Pi = {90.15393, -60.44028, -18.19041,42.34210} 

P2 = {75.27949, 86.12082,19.15087, -95.80345} 

P3 = {14.34134,2.00088,87.56698,39.80553} 
mi = m 2 = 21.23407, m 3 = m 4 = 81.40164 


Pi = {56.88939, 87.04163, -34.62173, -42.86104} 
P2 = {92.86718, -91.88334, 59.75945, 38.70047} 
P3 = {55.98527,-35.20008,9.02722,82.97219} 
mi = m 3 = 67.88777, m 2 = = 40.77317 


Pi = {95.95213,65.25140,-40.62468,30.93648} 
P2 = {68.47023, -60.09584,18.23998, 84.29507} 
P3 = {12.99839,12.08603, -99.08246, -34.58997} 
mi = m 2 = m 3 = m 4 = 11.50163 


pi = {62.72107,-90.39721,46.19119,67.90895} 

P2 = {57.27070,84.69403,21.86256,93.85136} 

P3 = {97.74567,11.18811, 68.32347, -31.80762} 

Pa = {51.13803, -98.84980, -16.22424, -84.44255} 
mi = m 2 = m 3 = m 4 = ms = 5.12753 
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C.l. Individual points 


Point 11 Pi = {33.74515,45.72730,31.15254,-7.47943} 

P 2 = {31.36435, -41.50734,46.47897,2.04203} 

P 3 = {4.59005,17.07010,32.65403,41.93628} 

Pa = {29.51054, -28.25963,46.17333, -35.08918} 
mi = m 2 = m 3 = 777,4 = ■'77-5 = 5.01213 

Point 12 Pi = {33.76482,45.44063,-10.68084,16.41925} 

P 2 = {72.93498, 67.49170, -11.81485, -36.28455} 
P3 = {8.01673, -49.40112, -66.09200, -0.11414} 

P4 = {-86.54188,-97.01228,68.12494,32.94875} 
mi = 98.42704,7772 = 28.89059,7773 = 40.51436 
7774 = 75.45643,7775 = 11.08327 

Point 13 Pi = {1.58374,6.86200,-15.06805,-10.63574} 

P 2 = {7.54800, -3.36539, 34.57385, 27.52676} 

P 3 = {43.36396, -49.27646, -25.35062, -17.68709} 
Pa = {22.58103, 38.31530, -14.67581, -3.08209} 

7771 = 'fn 2 = 'ITT’S = ITIA = ''^75 = 2.76340 

Point 14 pi = {-93.06712,-36.37997,-27.71460,38.42206} 
P2 = {-46.33465,-11.90909,32.33395,46.42742} 
P 3 = {8.41724, -83.92296,56.21715,34.04937} 

= {-15.23696, 71.33931, 48.68306, -53.67870} 
7771 = 59.10425,7772 = 60.25099,7773 = 76.79109 
7774 = 65.27606 ,7775 = 5.99925 

Point 15 pi = {-32.14401,-64.50445,46.04455,-75.56462} 
P2 = {-96.90340, -27.60002, -71.50486, 86.25541} 
P3 = {-37.95135,46.18586,25.67520,-71.38501} 
P4 = {-87.67870,66.66463, -36.20151, -27.37362} 
7771 = 7772 = m 3 = 79.63229,7774 = 7775 = 51.70237 

Point 16 Pi = {4.38045,-34.94383,6.34750,-24.71187} 

P 2 = {37.32326, -11.93504, -29.29012,31.65489} 
7771 = 7772 = 7773 = 7.27906 

Point 17 Pi = {-42.91300,62.46030,42.97296,-18.69577} 
P2 = {0.56662, 67.82929,16.29586, -8.45274} 

7771 = 85.86724,7772 = 81.52322,7773 = 76.12179 

Point 18 Pi = {31.22641,-11.44458,11.90191,5.51688} 

P2 = {39.14661,48.24103, -22.72685,18.32313} 

7771 = 7772 = 7773 = 6.89867 
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Appendix C. Momenta of the Example Points 


Point 19 Pi = {91.31175,59.18576,-1.99192,-24.84234} 

P2 = {71.50353, -58.97396, -27.79064, -4.80365} 
mi = 2.21844, m 2 = 36.02116, m 3 = 19.04975 

Point 20 Pi = {-2.85715,23.32734,73.44393,-19.93952} 

P2 = {-78.71545,23.81905,-94.38603,7.22745} 

P3 = {75.48933, -0.49095, -37.15293, -43.56893} 
mi = m 2 = m 3 = m 4 = 99.24288 

Point 21 pi = {96.38513,72.23517,25.81581,23.30879} 

P 2 = {-7.41341, -96.03161, -79.52827, -28.56493} 
P 3 = {-5.92534,98.48338,97.62395,-72.38643} 
mi = 32.69503, m 2 = 95.90407, m 3 = 53.60699, m 4 = 

Point 22 pi = {18.47170,12.68290,4.34693,-11.26827} 

P 2 = {-69.33772,-87.02423,59.44602,4.30660} 

P3 = {-91.88519,10.31462, 20.61165, -46.63546} 
mi = m 2 = m 3 = m 4 = 9.94045 

Point 23 pi = {-72.75921,-40.60802,-46.57923,73.25257} 
P 2 = {-92.33823, -18.31403, -5.67422, -86.71053} 
P 3 = {-76.47950,44.67795, -51.72040,57.95624} 
mi = 44.76393, m 2 = 29.71647, m 3 = 48.95514, m 4 = 

Point 24 Pi = {69.70234,62.68042,25.44429,-97.78603} 

P2 = {-65.98494,-85.19920,98.05702,-70.89141} 
P 3 = {-26.75642, -30.42288, -26.84633,14.81944} 
PA = {-69.44800,56.74842, -32.23649,96.45829} 
mi = m 2 = m 3 = m 4 = ms = 87.00572 

Point 25 Pi = {-45.80756,95.63842,-55.04954,44.01174} 

P2 = {36.09562, 52.66752, -11.22354, -87.48918} 

P3 = {-4.90798,41.11273,14.29379,2.15944} 

Pa = {49.48233,40.26756, -23.16581, -96.89362} 
mi = m 2 = m 3 = m 4 = ms = 56.97318 

Point 26 Pi = {-18.90057,-97.14671,44.69176,-16.67528} 
P 2 = {-70.86315, -81.27489, -3.71628,18.79403} 

P 3 = {-89.53092,50.02356,33.39784,-51.66031} 

= {-96.59097,-34.80215,-83.24353,44.73888} 
mi = m 2 = m 3 = m 4 = ms = 43.87459 


42.59971 


72.76312 
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C.l. Individual points 


Point 27 Pi = {-88.70322,37.98826,62.19352,-35.86433} 
P2 = {-58.60617,-58.60074,-83.75298,61.78210} 
= {-83.73607,46.98912,67.44602,78.40612} 

P 4 = {-96.41508, 71.69925, -14.47818, -61.82390} 


mi = m 2 = m 3 = m 4 = ms = 16.73899 


Point 28 pi = {-3.43584,4.73492,17.31242,61.53467} 

P2 = {12.12233,32.23256,87.57836,-58.25073} 

P 3 = {-38.67209,-54.27020,21.15570,79.15640} 

Pi = {-90.90573, -79.70266, -88.26463, -66.00973} 
= {-34.40043,-88.73043,84.41781,-4.21221} 
mi = 54.36459, m 2 = 30.96600, m 3 = 51.03652, 


m 4 = 16.03115, ms = 2.25657, me = 59.45020 


Point 29 


= {-9.85384,15.70678,80.94234,-84.96387} 

P2 = {90.11707, -74.59469, -70.73997, 54.32748} 

P 3 = {-55.84212, -34.47531, -87.20597, -27.73882} 
Pi = {16.72808, 64.83574, -31.16733, 63.94189} 

= {-42.62943,49.91058,-46.12974,59.76096} 
mi = 42.61768, m 2 = 22.13590, m 3 = 34.87263, 
m 4 = 54.00634, ms = 79.54844, me = 87.50131 


Point 30 Pi = {35.27512,36.08798,-89.66662,18.22907} 

P2 = {-32.58939,14.45447, 86.93898, -47.20827} 
P 3 = {-76.40210,-62.22587,-63.59955,41.03465} 
Pi = {-2.30248,0.45058, -76.74256, -64.19292} 

= {-88.80252,18.06504, -6.53891,49.34535} 


mi = m 2 = m 3 = mi = ms = me = 82.87370 


Point 31 Pi = {-99.20747,-68.16217,95.24772,68.87644} 
P2 = {-95.09224,78.51258,-82.38270,20.36899} 

= {-56.04092,22.93681,-72.82681,96.81954} 
Pi = {78.53840, -86.40143, -82.49674, -57.42855} 
P5 = {13.70265,77.87278,99.79126,8.31677} 
mi = 63.23680, m 2 = 86.48449, m 3 = 44.51361, 


mi = 79.73599, ms = 74.43246, me = 70.11421 
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Appendix C. Momenta of the Example Points 


Point 32 Pi = {-70.26380,96.72681,21.66556,-37.40054} 

P2 = {-13.45985,2.12040,3.20198,91.44246} 

P 3 = {-62.59164, -29.93690, -22.16595, -58.38466} 
= {-67.60797, -83.23480,18.49429,8.94427} 

= {-34.70936, -62.59326, -60.71318, 2.77450} 
mi = 94.53242, m 2 = 64.45092, m 3 = 74.74299, 


m 4 = 10.63129, ms = 31.77881, me = 23.93819 


C.2 Scans 


Figure 7.13 pi = {27.95884,25.55639,-29.88288,-2.17433} 

P 2 = {27.45521,-7.81292,3.19651,6.05088} 

6.05088 < mi = m 2 = m 3 < 31.53414 

Figure 7.16 pi = {67.40483,49.44993,-20.67085,48.63654} 

P 2 = {54.64295, -58.23071, 9.55042, -16.59411} 

P 3 = {41.37620,11.75178, -40.77655, -8.25014} 

2.33822 < mi = m 2 = m 3 = < 70.14658 

Figure 7.17 pi = {56.91015,-1.53764,58.47812,-64.49872} 

P 2 = {74.07761, 77.46365, 97.18347, -63.75152} 

P3 = {32.85182, -13.87036, -68.91291,87.62289} 

3.29978 < mi = m 2 < 98.99340 
m 3 = m4 = 22.70133 

Figure 7.18 pi = {-44.68568,-76.98829,54.10542,-72.67533} 

P2 = {30.54785, -3.09363, -50.03807, 88.12743} 

P 3 ,min = {71.46086,28.68398,47.82601, -93.53563} 

P 3 ,max = {3265.93916,860.51940,1434.78028, -2806.06893} 
mi = 68.14502, m 2 = 85.10360, m 3 = 75.43127, m 4 = 50.79371 

Figure 7.21 pi = {86.88039,29.81650,28.58637,88.04167} 

P 2 = {79.07341,18.41517, -37.32140, -80.98985} 

P 3 = {48.59046, -32.90460, -14.18821,85.89247} 

P 4 = {26.61931, -21.93050, 21.22092, 26.23208} 

3.16786 < mi = m 2 = m 3 = m 4 = ms < 95.03580 
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C.2. Scans 


Figure 7.22 pi = {-97.64085,75.52538,15.99594,35.63304} 

P2 = {-99.78509, -71.86820, -55.31804, -58.19254} 

P 3 = {-87.92245, -85.76444,19.12306, -59.11598} 

P 4 = {-87.70477,83.72569, -0.57040,35.79677} 

1.39548 <mi< 41.86441, 

m2 = 35.09895, m 3 = 74.42515, m 4 = 15.47168, ms = 41.24757 

Figure 7.23 = {-15.22437,-26.74156, 6.65483, 29.13661} 

Pi,max = {-83.38733, -146.46961,36.45001,159.58778} 

P2 = {-91.22611,-63.97875,55.07507,-52.90153} 

P3 = {0.95105, 75.90791, -10.13814, -88.40860} 

Pa = {43.04908, 77.11321, -50.69469, -7.60198} 
mi = 49.12560, m 2 = 57.87487, m 3 = 26.47098, 
m 4 = 0.42094, ms = 62.31320 

Figure 8.1 pi = {49.84064,5.20440,40.73789,-27.57262} 

P2 = {61.25771, 32.42640, -14.71631, -3.67014} 

1.92598 < mi = m 2 = m 3 < 57.77950 

Figure 8.2 pi = {-64.20450,59.44455,7.63891,57.71454} 

P2 = {1.91777,-70.66999,4.99626,48.86490} 

P3,min = {-95.66288, -30.14530,14.13910, -68.75167} 

P3,max = {-2292.42604,-904.35922,424.17308,-2062.55008} 
mi = 67.86064, m 2 = 88.88038, m 3 = 74.26926, m 4 = 34.85441 

Figure 8.3 pi^min = {-51.76504, -81.75539, -46.42422, -40.15540} 

Pi,max = {-283.52879, -447.79274, -254.27593, -219.94019} 
P2 = {-63.76533, -2.53015,16.27485, 69.16770} 

P 3 = {-78.50262,46.32052,13.19246, -54.00166} 

Pa = {25.40582,81.48058,39.11105,93.24648} 
mi = 78.45208, m 2 = 42.71315, m 3 = 91.94256, 
m 4 = 61.59730, ms = 16.75672 
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